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ABSTRACT 


A  previous  reliability  research  effort  by  Battelle  generated  an  analytical  technique 
called  the  Moment  method  for  predicting  drift-type  failures  in  components  and  systems. 

This  technique  is  based  on  the  Propagation  of  Variance  formula,  which  combines  data 
derived  from  the  system  or  system  model,  and  parameter  tolerance  data,  to  formulate 
an  estimate  of  the  probability  that  the  system  will  perform  within  specifications  at  the 
initial  or  some  future  time  of  operation. 

This  report  discusses  the  application  of  the  reliability-prediction  technique  to  a 
momentum- exchange-type  attitude  control  and  stabilisation  subsystem  for  a  satellite,  a 
liquid-metal  servo  actuating  subsystem,  and  an  electrohydraulic  servo  valve. 

A  mathenaatical  model  was  completed  for  the  momentum-exchange-type  attitude 
control  (fluid  flywheel),  significant  system  parameters  identified,* the  analysis  of  the 
expected  system  degradation  over  time  was  made,  and  a  prediction  of  efficiency  and  I 

reliabilitv  is  presented  in  the  report.  I 

i 

Performance  characteristics  used  to  define  system  performance  for  the  fluid-  I 

flywheel  system  were:  | 

(1)  Settling  time  for  an  input  step 

(2)  Open-loop  peak  efficiency.  j 

I 

I 

Since  parameter  variability  data  as  a  function  of  time  were  essentially  nonexistent,  it 
was  estimated  that  system  gain  and  fluid-friction  factor  would  be  the  system  param¬ 
eters  most  likely  to  drift  over  time.  (A  test  plan  designed  to  obtain  the  desired  param¬ 
eter  variability  data  is  appendixed  to  the  report.)  Settling  time  was  observed  to  be  a 
discontinuous  function  of  both  system  gain  and  friction  factor  due  to  the  definition  of  i 

settling  time  rather  than  a  physical  phenomenon.  Open- loop  peak  efficiency  for  the 
flywheel  is  defined  in  the  same  way  as  the  efficiency  for  an  a-c  or  d-c  servo  motor. 

The  results  of  the  investigation  give  the  parameters  to  which  their  performance  criteria 
are  most  sensitive  and  the  percentage  contribution  of  each  parameter  to  the  performance 
criteria  variances. 

The  pump  and  actuator  portion  of  the  liquid-metal  servo  actuating  subsystem  were 
similarly  modeled,  sensitive  system  parameters  identified,  and  critical  points  of 
degradation  estimated  from  analysis  of  the  Propagation  of  Variance  study.  Study  of  the 
ar.alog  model  for  the  pump  and  actuator  yielded  the  partial  derivatives  of  system  and 
component  performance  characteristics  with  respect  to  those  parameters  that  were 
identified  as  being  susceptible  to  degradation.  The  results  indicated  that  wear  in  the 
actuator  piston  rings  is  probably  the  most  critical  area  in  the  system  with  respect  to 
drift  reliability.  This  is  based  on  the  assumption  that  the  rod  seal  and  bearing  in  the 
actuator  perform  satisfactorily.  The  servo  control  valve  was  not  analysed  because  the 
design  for  this  component  is  incomplete. 

! 

The  servovalve  was  modeled  on  the  analog  computer,  and  the  sensitivity  of  per-  j 

formance  characteristics  with  respect  to  variations  in  the  physical  parameters  was  | 

investigated.  Design  optimisation  was  achieved  through  reduction  of  those  partials  of 
the  Propagation  of  Variance  Formula  that  had  large  values  and/or  were  related  to 
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parameters  whose  time  rate  of  chu^nge  wa<«  expected  to  be  significant.  The  second  step 
in  this  study  was  the  development  of  a  quantitative  prediction  of  the  relative  reliability 
improvement  expected  irom  a  new  va.'’e  design  with  respect  to  the  characteristic  of 
peak  quiescent  flow.  The  third  activity  on  this  task  was  the  performance  of  life  tests 
on  two  modifled  and  two  unmodified  servo  valves.  The  modified  valves  were  those  con¬ 
structed  in  accordance  with  the  findings  of  the  analysis  of  the  analog  model.  The  pre¬ 
dicted  improvements  in  quiescent  leakage  and  frequency  response  were  verified  by  the 
laboratory  teste.  In  addition,  a  significant  improvement  in  feedback-wire-tip  wear 
occurred  in  the  redesigned  valves. 
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INTRODUCTION 


Moat  present  reliability  analysis  and  estimating  techniques  are  based  on  statistical 
techniques  and  catastrophic  failures.  This  approach  is  not  appropriate  ior  electro¬ 
mechanical  devices  and  small-quantity  production  items,  such  as  advanced  vehicles, 
which  are  designed  to  be  subjected  to  extreme  environments.  In  addition  to  small- 
quantity  production,  the  development  and  evaluation  of  new  concepts  and  designs  for 
advanced  missions  and  environments  are  constrained  by  increasing  reliability  require¬ 
ments,  decreasing  calendar  time  for  the  evolution  of  a  design,  more  cfnnplex  functions, 
and  smaller  packaging.  New  methods  must  be  developed  to  predict  the  reliability  of 
new  and  critical  flight- control  components  having  unique  design  and  application,  small** 
production  quantities,  and  relatively  short  design  lead  tVne. 

Battelle  has  been  conducting  research  on  the  development  of  the  Moment  method 
of  analysis  of  drift  failures  in  electromechanical  devices  and,  under  the  present  con¬ 
tract,  has  undertaken  the  task  of  extending  and  modifying  the  techniques  developed  and 
applying  these  techniques  to  three  specified  systems:  a  momentum-exchange-type 
attitude -control  and  stabilisation  subsystem,  a  liquid-metal  servo-actuating  subsystem, 
and  an  electrohydraulic  servo  valve.  These  systems  and  components  have  been 
simulated  mathematically  on  the  analog  computer.  The  results  show  the  relationships 
between  system  performance  characteristics  and  individual  parameters,  thus  praviding 
information  for  design  optimization  and  reliability  predictions. 

Although  the  technique  may  be  considered  as  still  in  its  basic  stage  of  develop¬ 
ment  with  some  simplifying  assumptions  required  to  implement  it,  the  applications 
demonstrated  in  this  project  indicate  its  feasibility,  as  well  as  providing  guidelines 
for  similar  applications  to  other  systems.  The  technique  described  is  of  proven  value 
for  design  analysis  and  optimization,  and  reliability  prediction  of  drift  parameter 
degradation.  Battelle  has  been  conducting  research  and  development  of  the  technique 
for  3  years  under  Air  Force  contract.  The  initial  study,  conducted  under  Contract 
No.  AF  33(657)-7978,  to  develop  the  analytical  techniques  was  reported  in  Technical 
Documentary  Report  No.  ASD-TDR-63'287  dated  January,  1963.  The  research  was 
continued  through  the  application  and  further  development  of  the  techniques,  using  as 
an  example  an  electrohydraulic  servo  valve  and  an  associated  flight-type  linear 
hydraulic  actuator.  This  work,  conducted  under  the  same  contract,  is  reported  in 
Technical  Documentary  Report  No.  FDL-TDR-64-50  dated  August,  1964. 

This  report  is  divided  into  three  distinct  sections,  each  representing  a  complete 
md  independent  report  for  its  respective  task.  JLists  of  symbols  are  applieable  only  to 
the  particular  section  of  the  report  with  which  they  are  associated.  Brief  abstracts, 
stimmaries,  and  conclusions  are  contained  within  each  of  the  three  report  sections. 
Supplemental  material  containing  information  which  may  be  considered  proprietary  by 
other  contractors  has  been  delivered  to  the  program  Technical  Monitor.  Included  in 
this  supplemental  material  are: 

(1)  Equations  for  Mathematical  Model  of  Centrifugal  Pumps 

(2)  Derivation  of  a  Method  for  Determining  the  Nonsteady-State  Power 
Requirements  for  a  Centrifugal  Pump 

(3)  Approximate  Method  for  Determining  the  Pressure  Differential  Across 
the  Faces  of  a  Centrifugal  Impeller  Blade 
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(4)  AbaEc^  Computer  Wiring  Oiagram  and  Potentlomater  Setting*  for 
Pump  Simulation  Program 

(§)  Ralatipnahip  Between  Volume  of  Eroded  Material  ant  Increase  in 
Iap  S^e«, 

SUMMARY 


Tbia  raseateb  *uiie  been  concerned  primarily  with  the  application  of  the  Moment 
metitod  el  vdeiabitity  analysis  as  a  means  of  reliability  prediction  applied  to  three  speci¬ 
fic  fli|^t-contrel  ^etems.  Emidiasis  was  on  the  application  of  this  previously  developed 
teeJmiiin*  to  provide  a  basis  for  ethers  to  utilixe  it  with  their  own  systems.  The  analysis 
is  ecmcemed  with  the  gradual  shift  in  one  or  more  performance  characteriatics  as  a  re¬ 
sult  of  internal  parameter  variations  produced  by  environmental  and  time -dependent  ef- 
Cects,  thus  causing  the  performance  to  shift  outside  of  specification  limits.  The  goal 
then  is  to  predict  the  probability  that  a  device  will  exhibit  performance  characteristics 
within  ths  r*<piir*d  limits  for  a  specified  period  of  time  in  the  specified  environments 
The  Moment  method  is  based  oh  the  Propagation  of  Variance  formula,  equating  total 
variability  for  each  performance  characteristic  to  the  sum  of  the  contributions  of  varia¬ 
bility  from  all  the  related  internal  and  input  parameters.  The  results  of  the  analyses 
idealify  the  relationships  between  the  input  parameters  and  the  output  characteristics, 
and  provide  a  weighting  of  these  relationships,  thus  enaUing  identification  of  those 
parameter*  most  likely  to  cause  severe  degradation  of  the  system  and  permitting  a  pre¬ 
diction  of  system  reliability  ovor  time  if  the  input  parameter  data  are  known. 

The  three  systems  studied  on  this  project  were:  (1)  momentum-exchange-type 
attitude  contrpl  and  stabilisation  subsystem,  (Z)  liquid  metal  (NaK)  servo- actuating  sys¬ 
tem,  and  (3)  an  electrohydraulic  seirvovalve.  The  research  was  conducted  in  three 
separate  tasks,  and  this  report  has  been  organized  accordingly. 

Task  1,  Momentum- Exchange- Type  Attitude  Control  and 
Stabilisation  Subsystem 

The  investigatloii  into  ths  drift  reliability  of  the  fluid-flywheel  system  considered 
the  system  performance  specified  by  the  performance  characteristics:  (1)  settling  time 
for  an  input  step  and  (2)  open-loop  peak  efficiency.  Settling  time  was  given  an  upper 
•Pacification  limit  of  14.  5  seconds,  and  peak  efficiency  was  given  a  lower  specification 
limit  of  3  percent. 

Parameter-variability  data  as  a  function  of  time  wer*  essentially  nonexistent. 
Battelle  generated  a  test  plan  designed  to  obtain  the  parameter- variability  data.  The 
test  plan  is  submitted  as  an  appendix  to  this  report. 

From  muaufaebarers'  information  and  past  Battelle  esqperience,  it  was  decided  that 
system  gain  and  the  fhiid-friction  factor  would  be  the  system  parameters  most  likely  to 
drift  over  time.  The  system  gain  was  assumed  to  decrease  linearly  over  time,  and  the 
friction  factor  was  considered  to  have  an  exponential  increase  over  time.  The  result*  of 
the  investigation  give  the  parameter*  to  which  the  performance  criteria  are  most  sensi¬ 
tive  and  the  percentage  contribution  of  each  parameter  to  the  performance  criteria 
variances. 


Settling  time  wa«  obaerved  to  be  a  dlecontlnuoua  function  of  both  ay  stem  gain  and 
friction  factor.  A  aharp  decreaae  in  aettling  time  occurred  after  1/2  year  of  operation. 
The  diacontinuoua  relationahip  la  a  reault  of  the  atandard  definition  for  settling  time  and 
not  the  reault  of  phyaical  phenomena,  which  indicates  the  need  to  consider  carefully  what 
performance  characteristics  are  to  be  used  to  specify  system  performance.  Initiailly, 
the  nominal  settling  time  was  13.25  seconds,  with  a  tolerance  of  dfel.OS  seconds  and  a 
reliability  estimate  of  0.9999.  After  5  years  of  operation,  the  nominal  settling  time  was 
11.75  seconds,  with  a  tolerance  of  ±2.56  seconds  and  a  reliability  estimate  of  0.9999. 
The  high  reliability  estimate  at^the  end  of  5  years  is  due  to  the  decrease  in  nominal  value 
that  resulted  from  the  discontinuous  relationship. 

Open-loop  peak  efficiency  for  the  flywheel  has  the  same  definition  us  that  for  an 
a-c  or  d-c  servo  motor  and  is  a  continuous  function  of  the  parameters.  The  nominal 
value  of  peak  efficiency  decreased  steadily  over  time,  and  its  tolerance  value  increased 
steadily  over  time.  Initially,  the  nominal  peak  efficiency  was  4.27  percent,  with  a 
tolerance  of  ±0.  576  percent  and  a  reliability  estimate  of  1.000.  After  5  yearv'  of  opera¬ 
tion,  the  nominal  peak  efficiency  was  3.  25  percent,  with  a  tolerance  of  ±1.  31  percent 
and  a  reliability  estimate  of  0.7164. 

The  Moment  method  for  predicting  drift-type  failures  requires  quantitative  data  on 
the  degradation  of  parameters  as  a  function  of  time.  Since  these  data  were  lacking  for 
this  system,  it  was  necessary  to  rely  on  qualitative  information,  which  proved  to  be  suf¬ 
ficient  to  make  reasonable  assumptions  for  the  analysis.  It  was  concluded  that  the  sys¬ 
tem  has  a  high  reliability  with  respect  to  the  characteristics  chosen  and  the  arbitrarily 
assigned  specification  limits.  It  has  been  demonstrated  that  the  Moment  method  has  the 
ability  to  provide  insight  to  a  designer  on  the  behavior  of  the  system  characteristics  be¬ 
ing  considered  as  a  means  of  performance  evaluation. 


Task  II.  Liquid-Metal  (NaK-77)  Servo- Actuating  Subsystem 


A  satisfactory  model  of  a  centrifugal  pump  and  linear  actuator  was  developed,  and 
manipulation  of  the  analog  model  yielded  the  necessary  partial  derivatives  of  system  and 
component  performance  characteristics  with  respect  to  those  parameters  identified  as 
being  susceptible  to  degradation. 

The  analysis  of  the  liquid-metal  servo-actuating  subsystem  has  shown  that  the 
variance  of  the  system  performance  characteristics,  in  this  example  the  actuator  fre¬ 
quency  response,  is  due  primarily  to  tolerances  in  the  effective  clearance  around  the 
actuator  piston  rings.  Tolerances  in  pump  parameters  contribute  only  15  to  25  percent 
of  the  total  variance  in  system  performance.  System  performance  is,  however,  sensi¬ 
tive  to  degradations  in  the  tip  radius  of  the  pump  impeller  and  in  the  suigle  of  the  blade 
tip,  and  these  parameters  produce  up  to  twice  the  change  in  system  performance  as  an 
equal  percentage  increase  in  the  piston- ring  clearance.  Leakage  between  stages  of  the 
centrifugal  pump  has  only  a  slight  influence  on  the  pump  and  system  performance. 

Study  of  the  pump  and  actuator  configurations  leads  to  the  conclusion  that  changec 
in  blade-tip  radiufe  and  angle  due  to  liquid-metal  erosion  would  be  very  slight.  Also, 
since  interstage  leakage  does  not  significantly  affect  frequency  response,  it  could  be 


•xpect«4  that  wear  in  th«  piston  rings  of  the  linear  actuator  would  be  the  most  critical 
wear  area  in  the  syetonn  (assuming  satisfactory  performance  of  the  actuator  rod  seal  and 
bearing  over  the  required  lifo  of  the  syetem).  If  this  is  not  the  case,  the  bearing  and 
seal  may  bo  the  tnost  critical  because  degradation  can  lead  to  a  '‘catastrophic"  type  of 
failure. 

.  The  results  of  this  study  indicate  the  ability  of  the  Moment  method,  of  variability 
analysis,  in  conjunction  with  mathematical  and  aiutlog  models,  to  predict  the  variation 
that  can  be  expected  in  the  performance  of  a  servo  control  system  and  in  the  individual 
system  components.  The  reliability  observations  provide  useful  design  guidance  during 
development  stages  of  new  equipment  by  establishing  a  basis  upon  which  to  reduce  or  in¬ 
crease  manufacturing  tolerances  on  parts  affecting  various  system  parameters. 

Further,  study  of  the  normalized  partial  derivatives  indicating  the  sensitivity  of  per¬ 
formance  to  changes  in  system  parameters  assists  in  pinpointing  critical  wear  areas, 
thus  identifying  subjects  for  reliability  improvement. 


Task  111.  Redesign  of  an  Electrohydraulic  Servo  Valve 


The  first'part  of  the  analytical  study  for  the  electrohydraulic  servo  valve  employed 
an  analog  computer  simulation  of  a  two- stage  servo  valve  to  obtain  partial  derivatives  of 
various  performance  characteristics  with  respect  to  valve  parameters  that  are  suscepti¬ 
ble  to  degradation.  The  partial  derivatives  were  calculated  and  plotted  for  different 
values  by  key  design  parameters.  These  plots  were  evaluated  and  their  interrelation¬ 
ships  noted  in  order  to  identify  design  changes  that  would  reduce  the  over-all  sensitivity 
of  the  valve  to  parameter  degradation  and  to  the  effects  of  parameter  variations  intro¬ 
duced  in  the  manufacturing  process. 

The  study  showed  that  the  existing  design  was  already  close  to  optimum  as  far  as 
the  critical  areas  of  degradation  were  concerned,  and  it  was  not  expected  that  design 
changes  would  bring  about  significant  improvement  in  the  gross  reliability  of  the  valve. 
The  critical  areas  were  identified  in  previous  programs  as  those  parts  of  the  spool  and 
bushing  that  influence  the  quiescent  leakage  at  null.  In  the  reliability- improvement 
derivation  for  quiescent  leakage,  it  was  shown  that  in  the  statistical  sense,  i.e.,  for  the 
averaged  histories  of  a  large  number  of  valves,  a  definite  improvement  would  occur  in 
time  to  failure  with  regard  to  allowable  leakage. 

The  experimental  life-testing  program  carried  out  on  a  sample  of  two  modified 
and  two  unmodified  valves  essentially  confirmed  the  above  conclusions.  Because  the 
spread  in  initial  characteristics  of  the  valves  was  relatively  great,  in  both  the  inter¬ 
sample  and  the  intrasample  sense,  and  because  of  the  small  sample  size,  the  observed 
improvement  cannot  be  confirmed  on  a  statistical  basis. 

The  design  changes  that  were  ultimately  derived  did  reduce  slightly  the  sensitivity 
of  quiescent  leakage  to  degradation  effects,  and  also  reduced  the  initial  value  of  quies¬ 
cent  leakage  and  improved  the  frequency  response  of  the  valve.  In  addition,  these  modi¬ 
fications  reduced  significantly  the  sensitivity  of  several  performance  characteristics  to 
degradations  in  the  first- stage  orifice  and  nozzle-flow  characteristics. 
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An  interesting  unexpected  difference  between  the  designs  was  a  noticeable  im¬ 
provement  in  the  wear  of  the  ball  on  the  end  of  the  feedback  wire  in  the  new  design.  It 
is  felt  that  this  improvement  would  lead  to  significantly  better  performance  of  the  "aged" 
valve  in  the  critical  areas  very  close  to  null. 

Since  the  existing  valve  was  already  at  a  near  optimum  design,  the  reliability  im¬ 
provement  anticipated  from  the  design  modifications  was  relatively  minor.  However, 
the  results  of  the  study  are  significant  because  the  successful  application  of  the  Moment 
method  of  variability  analysis  to  the  redesign  and  improvement  of  an  existing  control- 
system  component  was  successful,  and  both  the  feasibility  of  application  and  the  pre¬ 
dicted  results  were  confirmed  by  practice. 


REVIEW  OF  THE  MOMENT  METHOD  ANALYSIS 


There  are  two  ways  in  which  a  eyetem  can  fail.  The  firet  ie  caueetl  by  drifts  in 
system  parameter  values  occurring  over  a  period  of  time,  which  results  hi  overall 
system  performance  outside  of  specification  limits.  Ehcamples  of  tbW  type  of  failure 
would  be  changes  in  transistor  characteristics,  chemical  reaction  bev-een  a  fluid  and 
its  container  that  would  be  sufficient  to  change  the  characteristics  of  the  fluid  but  not 
weaken  the  container,  and  wear  between  moving  surfaces.  The  second  cause  of  failure 
is  some  radical  or  catastrophic  occurrence,  such  as  a  break  between  a  power  supply 
and  its  load,  a  rupture  in  a  hydraulic  line,  or  a  break  in  a  mechanical  linkage.  This 
reliability-prediction  technique  applies  only  to  the  drift  type  of  failure. 


A  detailed  discussion  of  the  technique  has  been  reported  in  previous  reports  and 
paperst  ^#2,3)*.  therefore,  only  a  general  description  will  be  given  here.  The  Propaga¬ 
tion  of  Variance  formula  (Figure  1)  states  that  the  variability  of  a  performance  char¬ 
acteristic,  Aj,  represented  by  its  standard  deviation  depends  on  the  sensitivity  of 

the  p':jrformance  characteristic  to  parameter  changes  represented  by  the  partials 
(dA^/dP^)p. ,  the  variability  of  the  parameters  Pj  represented  by  Op-,  and  any  cross 

coupling  between  parameters  represented  by  Pj.g. 


Information  Necessary  to  Perform  Analysis 


In  order  to  perform  the  necessary  computation  to  obtain  the  variability  of  a 
performance  characteristic,  a  mathematical  model  of  the  system  must  be  established 
so  that  the  partials  may  be  obtained.  Data  must  also  be  available  that  give  the 
variability  of  the  parameters.  To  obtain  a  reliability  estimate,  a  set  of  specification 
limits  that  defines  failure  for  critical  performance  characteristics  must  be  known  and 
the  shape  of  the  distribution  curve  for  the  performance  characteristic  must  be  deter¬ 
mined.  In  order  to  make  any  prediction  on  future  performance,  the  variability  of  each 
parameter  with  time  must  be  known. 

The  partials  can  be  determined  analytically  from  the  system  model  or  empirically 
by  programming  the  model  on  a  digital  or  analog  computer  and  varying  the  parameters 
one  at  a  time  to  show  the  resulting  changes  in  the  performance  characteristic.  If  the 
systv^m  has  any  nonlinearities,  the  use  of  a  computer  is  practically  a  necessity. 
Knowledge  of  parameter  variability  might  be  obtained  from  component  manufacturers, 
independent  tests,  or  knowledge  of  the  present  state  cf  the  art  of  component  variability. 
For  computational  purpoties,  the  tolerance  values  on  the  parameters  arc  considered 
equivalent  to  3  Op^  for  each  parameter.  The  tark  or  mission  for  a  system  often  pro¬ 
vides  the  specifications  for  that  system. 

The  shape  of  the  distribution  curve  is  not  obtained  from  the  Propiagation  of 
Variance  formula.  However,  it  is  known  that,  as  the  number  of  parameters  affecting 
a  performance  characteristic  increases,  the  distribution  curve  for  that  characteristic 
approaches  a  Gaussian  or  normal  distribution.  In  applying  the  reliability  technique, 
the  distribution  curves  of  the  performance  characteristics  are  considered  normal. 


*r.eferences  ire  on  page  129. 
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“J,  =  I  (S^).  z  z  (Si). 

^  .1=1  ^  Pj  ^  ’*r  **• 

2 

o.  Is  the  varl'^nce*  of  component  perfoimanee  characteristic  A. 
''i  ^ 

N  ;.s  the  number  of  contrlbuttn  parameters 

2 

a„  is  the  variance  of  contributing  parameter  P. 

5  ^ 

Pj  is  the  mean  value**  of  parameter  Pj 

(^] 

y j  is  the  first  partial  derivative  of  component  performance 
j  o  characteristic  A.  with  respect  to  parameter  Pj  evaluated 

atPj  '  ^ 

p„  is  the  correlation  coefficient  that  relates  paramet-rs  P 
and  P,  ' 

2 

*  or  =  Variance  is  defined  as  the  second  moment  of  a  distribution  about 
its  mean  value. 

<T  =  Standard  deviation  is  defined  as  the  square  root  of  the  variance. 

**  Mean  value  is  defined  as  the  first  moment  of  a  distribution  about  the 
origin. 


figure  1.  PROPAGATION  OF  VARIANCE  FORMULA 


The  variability  of  each  parameter  with  time  and  environment  is  very  difficult  to 
obtain  in  present  practice.  To  determine  these  values,  strong  reliance  must  be  placed 
on  previous  life  tests  for  components  and  component  parts.  On  occasion,  the  time- 
dependent  variability  values  (or  functions)  would  be  generated  from  knowledge  or  insight 
gained  from  experience  in  working  with  a  component  or  a  part.  When  the  partial 
derivatives  are  known,  worst-case  operation  can  be  obtained  and  compared  with  the 
performance  specifications.  (Worst-case  operation  is  system  operation  at  the  extreme 
values  of  each  parameter  that  maximizes  and/cr  minimizes  a  performance  characteris¬ 
tic,  depending  on  whether  the  specification  limits  considered  are  upper  and/or  lower 
bounds. ) 


Procedure  to  Obtain  Reliability  Estimate 


When  the  mean  value  of  the  performance  characteristic  Aj  has  been  determined 
from  the  mean  values  of  the  parameters  Pj  and  the  variability  of  the  performance 
characteristic  a^.  has  been  determined  by  the  Propagation  of  Variance  formula,  the 

following  steps  will  give  the  reliability  estimate: 


(1)  Assxime  that  the  distribution  of  the  performance  characteristic  is 
normal.  This  has  the  effect  of  defining  the  shape  of  the  distribution 
in  terms  of  known  values,  i.  e. ,  the  mean  (nominal)  value  and  the 
variance. 


(2)  Determine  two  values  nj__  and  ny  so  that  n^  *  mean  value  of  the  per¬ 
formance  (Aj)  minus  the  lower  specification  limit  (SLl),  the  difference 
divided  by  the  standard  deviation,  Oy^.,  which  is  defined  as  the  square 

'  2  ^  r  ^  -  SLjl 

j.  Therefore,  n^  = 


root  of  the  variance 


and 


"A. 

1 


SLy  -Xi 


°Ai 


by  a  similar  reasoning  for  the  upper  specification 


limit.  If  there  is  either  no  upper  specification  limit  (SL^)  or  no 
lower  specification  limit  (SLy),  then  let  ny  or  ny  respectively 
equal  infinity. 


(3)  From  a  statistical  table  for  the  normal  distribution,  determine  the 
portion  of  the  distribution  corresponding  to  ny  cind  ny.  The  sum  of 
the  two  is  the  probability  that  the  performance  characteristic  values 
will  fall  within  the  tolerance  limits.  This  probability  is  the 
reliability  estimate. 


When  the  partial  derivatives  and  the  variemces  are  normalized,  the  critical  sys¬ 
tem  parameters  and  the  proportion  each  parameter  contributes  to  the  variance  of  the 
performance  characteristic  are  easily  distinguished.  Critical  parameters  are  those 
to  which  the  system  is  most  sensitive. 
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Limitation  of  the  Technique 


The  use  of  the  Propagation  of  Variance  formula  to  determine  reliability  presents 
three  possible  areas  for  error  that  would  affect  the  results  obtained  by  this  technique: 

(1)  The  first  two  moments  of  a  parameter  distribution,  mean  value  and 
variance,  used  in  the  Propagation  of  Variance  formula  are  not  sufficient 
to  completely  describe  a  distribution  that  is  not  normal.  The  Propaga¬ 
tion  of  Variance  equation  can  consider  higher  moments;  however,  from 
model  and  data  accuracy  consideration,  only  the  first  two  moments  were 
used.  The  characteristics  of  non-normal  distribution  are  not  reflected 
in  the  Propagation  of  Variance  formula. 

(2)  The  partial  derivatives  used  in  the  Propagation  of  Variance  formula  are 
considered  constant  over  the  tolerance  value  of  the  parameter  (e3opj). 

In  actuality,  the  partial  for  any  parameter  may  be  a  function  of  the  param¬ 
eter  value  over  the  tolerance  interval.  A  plot  of  tho  partial  over  the 
tolerance  interval  may  give  a  carved  line  rather  than  the  straight  line 
signified  by  the  constant  partial  value.  The  degree  of  curvature  will 
atfect  the  results  of  the  technique,  since  some  straight-line  approxima¬ 
tion  must  be  made  to  calculate  the  performance  characteristic  variance, 

2 


(3)  Any  prediction  of  future  reliability  will  only  be  as  good  as  the  variation 
estimates  of  the  parameter  values  and  variances  with  time  and 
environment. 

Complete  details  of  the  technique  development  and  procedures  are  contained  in 
Technical  Documentary  Report  No.  ASD-TDR-63-287  (AD  412496),  dated  January,  1963, 
and  Technical  Documentary  Report  No.  FDl-TDR-64-50,  dated  August,  1964. 
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TASK  L  MOMEWTUM-EXCHaNGE-TYPE  ATTITUDE-CONTROL 
AND  STABILIZATION  SUBSYSTEM 


Introduction 


A  previous  rollabillty  rosMirch  effort  by  Battelle  generated  an  analytical  technique 
called  ^e  Moment  method  for  predicting  drift-type  failures  in  components  and  sys¬ 
tems.  This  technique  is  based  on  the  Propagation  of  Variance  formula,  which  com¬ 

bines  data  derived  from  the  system  or  system  model  and  parameter  variability  data  to 
fornnilate  an  estimate  of  the  probability  that  the  system  will  perform  within  specifica¬ 
tions  at  initial  or  future  times  of  operation.  The  probability  value  is  interpreted  as  the 
reliability  estimate. 

This  report  discusses  the  application  of  the  reliability-prediction  technique  to  a 
momentum-exchange-type  attitude-control  and  stabilisation  subsystem  for  a  satellite. 
Several  types  of  controllers  would  fit  into  the  momentum-exchange  category.  This  re¬ 
port  considers  the  type  called  a  fluid  flywheel.  The  principle  of  the  fluid  flywheel  is  to 
accelerate  a  fluid  within  a  closed  loop  in  the  satellite.  This  action  generates  a  torque 
about  the  axis  of  the  loop  to  which  the  satellite  responds  and  consequently  is 
repos  itiemed. 


Engineering  Activity 


Description  of  Fluid-Flywheel  System 

The  fluid  flywheel  considered  in  this  report  Is  one  designed  and  built  by  the 
General  Electric  Company.  General  Electric's  goal  was  to  prove  the  feasibility  of  a 
fluid  flywheel  as  an  a'.titude-control  device.  All  numerical  values  for  the  system 
parameters,  most  of  the  parameter  tolerance  values,  and  most  of  the  system  equa¬ 
tions  are  presented  in  the  General  Electric  report  on  the  fluid  flywheel.  There  is 
effectively  one  flywheel  for  each  of  the  three  axes  of  a  satellite,  all  powered  from  a 
single  electrical-energy  source. 


Functional  Block  Diagram  and  General  Discussion.  The  functional  block  diagram 
for  one  of  the  three  axer  of  the  fluid-flywheel  system  is  shown  in  Figure  2.  The  opera¬ 
tion  of  each  of  the  three  axes  of  the  flywheel  is  identical,  with  the  exception  ol  the  loadinj 
imposed  by  the  inertia  of  the  satellite  about  each  axis. 

General  ^'.lectric  tested  for  crcss-coupling  between  the  axes  and  found  that  no 
appreciable  amount  exists.*  Therefore,  Battelle's  investigation  considered  the  opera¬ 
tion  about  a  single  axis  with  no  cross-coupling. 

The  torque  to  maintain  a  given  attitude  under  disturbances  and  to  obtain  a  desired 
attitude,  given  a  position  error  indication,  is  generated  by  accelerating  a  fluid  in  a 


a,.  General  Electric  Cosn^'ny,  Jotnian  City.  New  Terfc,  tdephone  converutinn.  December  19Ct, 


closed  loop  (loops)  about  the  appropriate  axis  (axes).  The  action  is  similar  to  a  servo- 
motor-driven  inertia  wheel  for  attitude  control.  The  torque  available  for  cmitrol  is 
proportional  to  the  product  of  the  fluid  inertia  value  and  its  acceleration. 


FIGURE  2.  FUNCTIONAL.  BLOCK  DIAGRAM  OF  FLUID-FLTWKBEL  ATTITUDE- 
CONTROL  SYSTEM 


The  heart  of  the  fluid  flywheel  is  Ute  pump,  which  accelerates  the  fluid,  (ht  this 
flywheel  the  fluid  was  mercury. )  The  pump  is  an  electromagnetic  d-c  conduction  pump 
(E.  M.  pump),  which  operates  on  the  same  basic  principles  as  an  electric  d-C  motor. 
For  the  operation  of  this  type  of  pump,  the  fluid  must  be  a  liquid  metaL 


Brief  Discussion  of  Electromagnetic  Pump.  The  electromagnetic  pump  is  not  a 
new  device.  It  has  been  used  since  the  early  19S0's  for  pumping  liquid  metals,  partic¬ 
ularly  sodium  and  NaK  (mixtvre),  as  coolants  of  nuclear  reactors.  Oflier  liquid  metals 
considered  as  coolants  have  been  bismuth,  lithium,  and  magnesium.  The  basic  con¬ 
struction  of  the  pump  is  also  applicable  to  metering  liquid-metal  flow  rates.  A  liquid 
metal  that  flows  within  a  constant  magnetic  field  will  generate  a  current  proportional  to 
the  flow  rate.  Electromagnetic  pumps  vary  in  sins  from  those  pumping  less  Uian 
1  gpm  and  requiring  several  amperes  to  those  pumping  30  gpm  and  requiring  several 
thousand  amperes.  The  development  of  electromagnetic  pun^  has  received  good 
coverage  in  the  literature.  ^ 

The  pumping  of  the  mercury  is  achieved  by  passing  an  electric  current  through 
the  mercury  under  the  influence  oi  a  magnetic  field.  The  mercury  flew  path,  the 
cxurrent  path,  and  the  magnetic  field  are  mutually  orfliogonal.  The  interaction  of  the 
current  and  magnetic  field  applies  a  force  to  the  mercury,  causing  the  mercury  to 
move.  The  back  electromotive  force  of  the  pump  and  the  fluid  losses  in  the  loop  pro¬ 
vide  sufficient  damping  so  that  no  rate  feedback  or  other  compensatioa  is  required. 


System  Electronics.  The  force  that  accelerates  the  mercury  is  proportional  to 
the  current  flow  geuerated  by  the  error  signal  from  the  attitude  error  sensor.  The 
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•igMkl  i*  by  a  Toltege-to~curr«nt  Amplifier.  The  load  presented  to  the 

■WipWftar  by  tlta  merttury  Is  In  tike  mUliohm  range.  Therefore,  a  prime  design  goal  is  to 
bbey  tite  outyot  tnafedence  of  tiie  amplifier  as  small  as  possible. 

The  amfiiflev  tased  eras  designed  by  General  Electric,  and  Is  based  on  the  use  of 
parallel  transistors  to  reduce  the  output  impedance  seen  by  the  pump.  The  amplifier 
echematie  diagram  Is  shoem  in  Figure  3.  Detailed  discussions  of  the  voltage>to~current 
ana^lifier  and  Its  devel^nent  have  been  recorded  previously  by  General  Electric 
authors.  (■*» 

Efficiency  of  the  fluid-flywheel  system  is  quite  low  due  to  the  fluid  losses  in  the 
loop  and  electrical  losses  at  tiie  amplifier-pump  interface.  The  amplifier  gave  approxi¬ 
mately  a  4  to  I  ratio  between  an^lifier  output  impedance  and  mercury  load  resistance. 
Improved  high-current,  low-voAtage  characteristics  of  new  transistors  give  promise  of 
amplifier  output  impedance  to  mercury  resistance  ratios  of  0. 4  to  1. 

System  Model 

A  model  of  the  fluid  flywheel  was  reported  by  Nichol.  This  model  considers  the 
eeneor-amplifier  combination  as  an  ideal  current  amplifier  with  a  gain,  a  small  dead 
band,  and  a  saturation  level  that  limits  the  current  to  60  amperes.  In  most  respects, 
Nichol's  model  parallels  the  models  described  in  previous  reports.  9, 13) 

A  sketch  of  the  electromagnetic  d-c  conduction  p*imp  is  shown  in  Figure  4.  The 
basic  equation  for  eyetem  operaticn  is  the  vector  force  equation 

F  a  i  w  X  B 
~  ”P  ” 

for  the  force  exerted  on  the  fluid  to  accelerate  it, 
where 

1^  is  the  current  in  the  pumping  channel 
w  is  the  width*  of  channel 
B  is  the  magnetic  flux. 

With  the  current  path  at  a  right  angle  to  the  flux  field  and  using  the  pound-feet  system  of 
units: 

F=i  wxB  =  6.  87  X  10-6  Bi  w  lb,  (1) 

—  —p  —  p 

with  B  in  gauss. 


The  pressure  developed  in  the  pump  is 


_  .  m  _ 

p  «  JL  =  6.  87  X  lO'®  Ib/ft*  , 

t  wb  h  ' 


and  can  be  written  in  terms  of  the  voltage  generated  by  the  amplifier  by  considering  the 
circuit  shown  in  Figure  5. 


>4cComb.  R. ,  Coer*)  electric  Comp«ny,  Schenccudy,  New  York.  B*tuU«-G.E.  conference,  October  6,  1904. 
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FIGURE  5.  EQUIVALENT  CIRCUIT  OF  AMPLIFIER- PUMP  INTERFACE  (BATTELLE) 


Figure  5  represents  the  interface  between  the  amplifier  and  the  pump.  From 
Figure  5  and  standard  circuit-analysis  procedures,  it  can  be  shown  that 


and 


K  0  R,  “  E.  (R  +  R*) 
V  «DB  f  «»  *  I 


R^Rf  Rj^Rp  +  RpR£ 


(3) 


6.87  X  10'^  Bip 


6.  87  X  10*^  B 
h 


Rf  -  Eb  (Ra  +  Rf) 
R^R£  -t-  R|^Rp  *  RpRf 


(4) 


An  expression  for  Eb  is  needed  in  terms  of  the  system  parameters  and  the  fluid  angular 
velocity,  <ii£,  and  is  obtained  from  the  vector  equation 

Eb  »  back  emf  »  ^  (V  x  B)  •  dL  .  (5) 

The  fluid  velocity,  the  flux  density,  B,  and  the  length,  L,  of  the  current-carrying 
medium  (which  in  this  case  is  the  pump  channel  width)  are  mutually  orthogonal,  and 
therefore 


Eb  -  9.  3  X  10-6  BwV. 


(6) 


The  fluid  velocity  through  the  pump  channel  can  be  esqkressed  in  terms  of  the  fluid 

2 

D. 

(7) 


angular  velocity  in  the  loop,  uy,  by 

V 


TTDjDjWy 


8hw 


giving 


Eb 


9.3  X  10-4  7rBD|P,Mf 
8K 


Nichol  defines  an  overall  diffuser  ratio.  A,  as 


A  ■ 


4hw 


(8) 
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Th«  diffuse r  ratio  should  be  greater  than  one  to  keep  the  magnetic  gap  small  and  to  ob¬ 
tain  a  large  magnetic  flux  with  small- sised  magnets.  Substituting  A  into  Equation  (8) 

fives 


9. 3  X  10"®  BAw  Di4>»f 
2 


The  expression  for  ip  is  now 


R.  -  4.65  X  10-6  bAwD.  w£(R  +  Rf) 

i  .  DB  _ _  .  {V 

^  Ra®f  ^  ^^^f 

Niehol  does  not  consider  the  amplifier- pump  interface  to  be  represented  by  the  circuit 
in  Figure  5  but  rather  that  shown  in  Figure  6. 


FIGURE  6.  EQUIVALENT  CIRCUIT  OF  AMPLIFIER- PUMP  INTERFACE  (NICHOL) 
From  Figure  5  and  standard  circuit-analysis  procedures,  it  can  be  shown  that 

Kj  R-  .  E. 
a  *DB  ^  ° 

I  s  ■  —I  I  '  , 

P  Rp  +  Rf 


and  substituting  for  Ef^, 


K-Q,  R,  -  4.  65  X  10"®  BAwDiW, 
* 

Rp+Rf 


The  expressions  for  ip  give  it  as  a  ftuiction  of  the  angular  error  after  passing 
through  the  sensor  dead  band,  of  the  angular  velocity,  hf.  However,  from 

Figure  5,  represents  a  voltage  source,  Es,  and  from  Figure  6,  Karens  repre- 

UB  iJn 

sents  a  curront  source,  I5.  The  question  arises,  "Can  the  interface  representation  in 

Figure  5  be  substituted  for  that  in  Figure  6  without  any  change  in  the  pressure  developed 

(Pt)?"  Niehol  gives 


6  87  X  10“6  R^B 
h(Rf  +  Rp) 


6.43  X  10"^^  ttB^dIdi 
8h^(Rf  +  Rp) 
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(13) 


■C~^  a*  *Vg  1.  -  f 

and  Battelle  gives 

6.87xlO"^RfB  ^  6.43  X  10'“  rrB*D|Dj(Ra  +  Rjf> 

‘  h(R^R£+R,.Rp+B^Rf)  8h^(R|iRf4R«Rp-t-RpRf) 

»  KgEs  -  Kb  wf . 

4t  )Qt 

The  question  arises,  "Does  K.«E5  -  K),  uf  «  Kgl  -  Kb  <i)f?" 

From  Figure  5, 

[R  R,  1  r9.3  X  10-6irBD,D,R.'|  ^ 

L  8MRf+V)  p  \  Z 

Substituting  the  expression  for  Eg  mto  Equation  (13)  does  give,  with  some  manipulation, 
Nichol's  expression  for  P  and  the  equality  hold?; 


(14) 


i^uua vu«  cj(p.r5BB&uD  inbu  jcL4u«iiioii  \a 

Nichol's  expression  for  and  the  equality  hold?; 

Vs'<“f’V”V: 


(15) 


It  is  necessary  to  determine  the  value  for  shown  in  Figure  5.  Nichol  states  that 
the  sensor-amplifier  combination  effectively  saturates. at  im  input  error  of  two  degrees, 
since  the  current  I  is  limited  to  60  amperes.  For  a  maximum  current  of  60  amperes, 
the  fluid  velocity  at  steady  state  is  1. 18  radians /second.  From  Nichol's  characteristic 
loss  curves  in  Figure  7,  it  can  be  determined  that 

H  1.1  ft  lb  sec  . 

w,  =  Y"  =  1 - S-’  =  1. '  8  radians /second. 

‘^f  0.  935  u  lb  sec 


Realising  that  during  transient  bpeifation'there  is  a  wide  range  of  possible  combinations 
for  the  input  error,  source  current,  and  fluid  angular  velocity,  Battelle  decided  to  de¬ 
termine  the  values  of  for  steady-state  operation  when  the  input  error,  0  ,  was  two 
degrees,  I  was  60  amperes,  and  was  1.18  radians /second.  This  can  be  accomplished 
by  using  Equation  (14),  where  Eg  is  equal  to  1^^^  »  substituting  in  the  known 

DB 

values  for  9  and  all  the  terms  on  the  right-hand  side  of  the  equation.  This  gives 
*DB 


K^.O.  355  volts/degrec. 


The  pressure 
flywheel  loop.  The 
as  follows: 

(1)  Pressure 


developed  by  the  pump,  P^,  is  opposed  by  the  fluid  drag  in  the  fluid - 
fluid  drag  is  the  sum  of  three  pressure  drops,  which  are  expressed 

2  2  2 

drop  in  the  pump  P^  = 


17 


(2)  Pressure  drop  in  the  diffuser  = 


,,  .  ,2^2  2 
(1  -  "f 


8g 


2 

Wf7F®l"f 

(3)  Pressure  drop  in  the  loop  Pj^  =  - - 


The  expressions  for  the  pressure  drops  are  correct  for  turbulent  flow  only. 
Battelle  has  checked  Nichol's  claim  that  laminar  flow  considerations  are  not  necessary 
and  found  this  claim  to  be  valid.  In  Battelle' s  analog  simulation  of  the  system,  'provi¬ 
sion  was  made  for  reading  out  the  Reynolds  number  as  a  function  of  time  during  a 
dynamic  response  test  of  the  flywheel.  A  plot  of  the  Reynolds  number  readout  is  shown 
in  Figure  8  for  the  system  excited  by  a  step  function  of  two  degrees.  During  the  re¬ 
sponse,  the  Reynolds  number  was  in  the  range  of  turbulent  flow  (R^  >  2,  000)  for  the 
critical  portions  of  the  response  (large  errors),  and  remained  in  the  laminar  flow  range 
for  very  small  values  of  error  toward  the  end  of  the  response. 


FIGURE  8.  REYNOLDS  NUMBER  VERSUS  TIME  FOR  A  SYSTEM  INPUT  STEP 
OF  TWO  DEGREES 

The  Fanning  friction  factor,  f,  is  given  as 

f^ 

(R^)0.16  * 

with  the  friction-iactor  constant  f  initially  equal  to  0.  031.  The  f  appears  in  the 
and  loop  drops.  The  Reynolds  numbers  for  the  pump  and  loop  are  given(d)  as 

R  (pump)  =  -  w, 

*  p  (h  +  w) 


pump 

(16) 


and 


i 
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(17) 


>  D  D 

\  (loop)  =  — Wf  . 

Substituting  th«  above  values  foe  the  Reynolds  numbers  in  the  appropriate  equations 
gives  the  following  relationohips: 

_  1,84 

- - • 


Defining  as 


K  S 
P 


0.34  _  1.84,.  ,1.16  0.16 

fg7p  Li(A  Dj )  (h  +  w)  jLi 


4g(hw) 


ITT?" 


(18) 


P  3  K  u 
P  P 


1.84 


(19) 


^  9.84_  2.84  0,16 

P  -  c^F  1  ^  1.84 

L"  ,,*1.16  _  1.16  “f  . 

(2)  g  D 


Defining  ** 


(,)i.i6^ 


(20) 


1.84 


(21) 


The  maximum  theoretical  diffuser  efficiency,  given  by  Nichol  to  be 


t]d  -  r“  ^  ^ ' 


(22) 


where 


X  = 


0.  88  [  G(R  )]  (A^  -  C^)  +  2  (C  -  J ) 

XT';  “ 


and  G(R  )  »  1  for  R  >  10' 
6  6 


G{R^)  =  0.  325  log  -  0.  625  for  2  x  10  R^  i  lO'*  , 


and 


C  ■  1  for  R  >  10" 
e 
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" ' « 


C  =  - 1 -  for2xlO^<R  <10^. 

0.088{C(R^))  ~ 

The  diffuser  efficiency  depends  upon  the  transition  from  the  rectangular  pump  aection  to 
the  cylindrical  tubing.  General  Electric  performed  numerous  tests  to  determine  the 
best  geometric  configuration  to  maximize  7)^.  (1^) 

The  expression  given  by  Nichol  for  r)^  shows  that  should  be  a  function  of  the 
Reynolds  number  and  would  consequently  vary  over  any  dynamic  response.  From  the 
analog  mechanization  shown  in  Figure  13,  the  variable  value  of  t]^  was  computed  over  a 
dynamic  response  and  inserted  in  the  expression  for  diffuser  pressure  drop.  The  re¬ 
sponse  of  the  fluid-flywheel  system  to  a  step  in  input  position  signal  when  was  varia¬ 
ble  was  similar  to  the  response  of  the  system  when  the  constant  value  for  7)^^  of  65  per¬ 
cent  suggested  by  Nichol  was  used.  Further  investigation  showed  that,  for  a  major 
portion  of  the  dynamic  response,  the  diffuser  efficiency  was  effectively  a  constant  value, 
but  equivalent  to  79.  5  percent,  not  65  percent.  Figure  9  shows  a  plot  of  system  error 
in  the  yaw  a.xis  for  a  step  input  to  the  flywheel  system,  with  equal  to  65  percent, 

T)j^  equal  to  79.5  percent,  and  tj  calculated  over  the  dynamic  re8i>onse.  The  plot  for 
^  of  65  percent  is  distinguishaole  from  the  calculated  7)^  and  the  7)^  of  79.  5  percent. 
Tne  plots  for  the  calculated  7^^  during  a  response  and  the  constant  of  79.  5  percent 
overlap  throughout  the  response  range.  As  a  result,  a  constant  Reynolds  number  of 
39,  000,  corresponding  to  a  nominal  value  of  7)^^  of  79.  5  percent,  was  used  in  this  inves¬ 
tigation.  It  was  concluded-  that  the  system  response  was  relatively  insensitive  to 
changes  in  Reynolds  number.  For  the  reliability  analysis,  the  diffuser  efficiency  value 
changed  for  different  dynamic  responses  due  to  changes  in  parameter  values  that  appear 
in  the  defining  equation  for  diffuser  efficiency. 


FIGURE  9.  COMPARISON  OF  DIFFUSER  EFFICIENCY  VALUES, 

The  change  in  the  nominal  value  does  not  change  the  diffuser  pressure-drop 
expression  in  its  literal  form: 
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'UN  - 


•1 


'23) 


'  - 3i -  “£  *  •'d“f  • 


The  Mt  pressure  to  seeelerste  the  mercury  at  say  time  is 


P  -  P*  -  P  -  P.  -  P, 

B  t  p  d  L 


e  e 

P^  »  K  _  E.  -  - 
n 


1.84 


(24) 


where  K  sad  K,  are  as  defiaed  in  Equations  (18)  and  (20). 

y  ^ 

The  pressure,  P^,  develops  a  torque  with  the  relationship 

T  «  P  (cross>sectional  area)  (moment  arm) 
a  a 


[-  (^)1  (^) 


TD^D 

»  —A  p 

8  n 

«d|d, 


T  »  K*  E-  -K.  w-  -  (K  +  K, 
n  a  S  b  f  P  Id 


,  ,  ,  1.84  •  2 


(25) 


'2('' 


The  last  two  terms  in  the  expression  for  'T  may  cause  some  cozifusion  since,  in  a 
dynamic  run,  the  mercury  flow  can  change  direction  indicating  a  positive  and  negative 
value  for  w^.  To  account  for  this  possibility,  the  expression  for  the  torque  should  be 
rewritten  as 

'  •'i  X  “f  -  <*'/  "L*  hf  ■  ”  -i”  “f  -  Kj  |-,1  ^  .g.  .  (27) 

Tn  ■  Tg  -  Tb  -  ITp  +  Tl)  -  T^. 

The  torque  accelerates  the  mercury  and  the  satellite  about  the  axis  of  the  mercury  loop. 
The  loading  effects  of  the  mercury  and  the  satellite  on  the  system  are  their  respective 
inertias.  Figure  10  shows  the  closed-loop  block  diagram  of  the  system  mode^  used. 

Nichol  considers  a  torque -momentum  expression  (as  shown  in  Figure  6),  rather 
than  the  torque -angular  velocity  expression  used  here.  The  expressions  are  equivalent 
in  that  momentum,  H,  is 


H  =  Jj  wj  ,  '28? 

the  product  of  the  mercury  inertia  and  its  angular  velocity,  and  H  /J^  is  sub'^tltute''*  for 
u)j  in  the  torque  expression. 
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BattfiU*  performed  dyimmic  roue  uelag  the  model  in  Figure  10  for  eomperieon  with 
tho  actual  ayetem  respoaeea  reported  by  Nlchol.  Figure  11a  ahows  the  cloaed-loop  re- 
apflaMoa  of  the  three  ay^em  axea  aa  reported  by  Nicho!,  and  Figure  11b  givee  the  model 
rocpoiiaea  aa  obtained  by  Battelle.  In  all  three  eaaes,  the  agreement  between  the  curvea 
ia  ancal^bMili* 

'  ■ 

Thwaamarical  paramatar  valuea  uaed  in  the  model  simulation  are: 
yp  »  iso  (lb)/ft2 
H  «  l.OdS  «  lO'^  Ib/ft  eee 
B  •  18,800  gauaa 
Dj  «  4  feat 

■  3  X  10‘^  ft 

h  s  6  X  10  ^  ft 
-2 

w  »  4  X  10  ft 

L  »  2.  08  X  lO"*  ft 

Rg  a  lO”^  ohms 

R^  •  2  X  lO”^  ohme 

R  /R,  a  3  X  lO’^ 

P  f 

(R  1-R,)  3  6.  33  X  lO  ^  ohma 
P  * 

■  79.  5  percent 

g  «  32.  2  ft/sec^ 

»  9.  35  X  lO'^  ft/lb/aec^ 

Jy  3  53.0  ft/lb/aec^ 

Dg  3  10  ^  degrees 

Vg  »  7.  04  X  10*'  volts 

K^a  3.  55  X  10  ^  volts/degree. 

The  flywheel-system  model  was  mechanixed  on  the  analog  computer  as  shown  in 
Figure  12.  The  computer  circuitry  for  calculating  diffuser  efficiency  is  shown  in  Fig¬ 
ure  13.  Figure  14  shows  the  computer  mechanisation  for  generation  and  readout  of 
power  efficiency,  which  is  discussed  in  the  following  sections  on  performance  criteria 
and  results. 
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FIGURE  11.  SYSTEM  RESPONSE  CURVES  FOR  STEP  INPUTS 
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FIGURE  14.  MECHANIZATION  FOR  GENERATION  AND  READOUT  OF  POWER  EFFICIENCY 


Performance  Criteria 


Other  than  the  ability  to  operate  in  a  space  environment,  no  specifications  were 
imposed  and  no  mission  was  given  for  the  fluid  flyu'heel.  General  Electric  did  operate 
the  flywheel  system  in  the  temperature  range  of  -19  to  90  F  and  noticed  no  change  in 
the  system-response  characteristics.  General  Electric's  temperature  range  of  opera¬ 
tion  does  not  signify  the  meeting  of  any  specification,  but  it  does  indicate  what  tempera¬ 
ture  variations  could  be  accepted  easily. 

The  lack  of  specifications  makes  the  choice  of  performance  criteria  somewhat 
arbitrary.  However,  in  light  of  the  previous  comments  on  damping  and  efficiency,  the 
following  two  performance  criteria  were  chosen: 

(1)  Settling  time  for  an  input  step  in  position.  This  criterion  is  for  a 
closed-loop  response  and  provides  a  measure  of  the  damping  of 
the  system. 

(2)  Open-loop  peak  efficiency.  This  criterion  presents  an  efficiency 
measure  similar  to  that  used  to  evaluate  a  servemotor. 


Settling  Time.  Settling  time  is  defined  arbitrarily  as  the  amount  of  time  it  takes 
the  satellite  position  to  settle  to  within  five  times  the  nominal  dead-band  value  when  ex¬ 
cited  by  an  input  step  in  an  angular  position  of  2.  5  degrees.  The  2.  5-degree  input  is 
large  enough  to  force  the  system  into  saturation,  and  small  enough  not  to  bury  the  in¬ 
tended  parameter  variations  in  a  saturated  response.  Battelle  considered  only  one  axis 
of  the  General  Electric  flywheel,  namely  the  yaw  axis,  with  moment-of-inertia  of 
53  foot  pound  seconds^.  The  dead  band  of  the  sensor  was  0.  01  degree;  this  determined 
the  settling  time  to  be  the  time  required  to  decrease  to  and  stay  within  the  error  bounds 
of  ±0.  05  degree. 


Open-Loop  Peak  Efficiency.  General  Electric  considered  efficiency  as  a  perfor¬ 
mance  criterion  of  the  system  and  mentioned  that  an  efficiency  of  4  to  5  percent  should 
be  expected  when  efficiency  was  defined  as 

.  power  out  (mechanical) 

Efficiency  =  * - : — j—, — — — rr -  . 

power  in  (electrical) 

Discussions  with  General  Electric  disclosed  that  no  standard  exists  with  which  to 
reference  the  characteristics  of  a  fluid  flywheel.  General  Electric  further  pointed  out 
thatf  while  the  definition  given  for  efficiency  is  good,  one  must  be  careful  in  talking 
about  flywheel  efficiency.  It  is  important  to  know  exactly  how  the  power  terms  are  ob¬ 
tained,  and  what  the  ground  rules  are,  in  any  discussion  of  efficiency. 

General  Electric  obtained  their  efficiency  values  by  a  method  used  for  determining 
the  efficiency  of  servomotor-driven  flywheels  (disks).  The  procedure  for  obtaining  effi¬ 
ciency  w?s  to  excite  the  flywheel,  at  rest  and  under  open-loop  conditions,  with  a  signal 
sufficient  to  obtain  maximum  momentum.  When  the  momentum  reached  one-t'r.lrd  of  its 
maximum  value,  the  power-in  and  power -out  measurements  were  made.  The  ratio  of 
these  two  values  was  determined,  giving  the  efficiency  value.  In  performing  the 
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measurements  necessary  to  obtain  efficiency,  no  standard-load  was  used,  which  permits 
a  possible  variation  in  efficiency  for  a  particular  flywheel  with  different  loads. 

Efficiency  and  momentum  curves  are  noxilinear,  imd  the  one-third  maximum 
momentum  as  the  reference  value  for  the  efficiency  measurement  is  somewhat  arbitrary. 
Battelle  decided  or  approach  for  determining  efficiency  using  the  General  Electric 
definition,  but  eliminating  a  particular  momentum  value  as  a  measurement  reference. 
Battelle  also  decided  to  use  open-loop  peak  efficiency  as  a  performance  criterion. 
Open-loop  peak  efficiency  is  obtained  in  precisely  the  same  manner  as  that  for  a  servo¬ 
motor.  it  Is  the  maximum  value  of  the  power -out  to  power -in  ratio  from  the  standard 
speed-torque  characteristic  curves,  independent  of  the  fraction  of  maxim\un  momentixm 
at  which  it  occurs. 

Battelle  obtained  its  power -out  term  directly  from  the  analog  simulation  by  taking 
the  product  of  the  mercury  angular  velocity  and  the  developed  torque,  and  multiplying 
by  the  conversion  factor  1. 356  to  obtain  the  power-out  in  watts: 

Power  out  =  P^,  =  1.  356  T  w,  watts.  (29) 

O  n  f 

The  power -in  term  cannot  be  obtained  directly  from  the  analog  simulation.  The 
problem  arises  from  the  fact  that  the  electronic  amplifier  was  simulated  as  a  Thevenin 
equivalent  circuit.  A  Thevenin  circuit  is  sufficient  for  giving  voltage -transfer  relation¬ 
ships,  but  a  power  value  obtained  from  the  product  of  the  Thevenin  circuit  voltage  and 
its  current  output  does  not  give,  in  general,  the  power  of  the  device  being  represented 
by  the  Thevenin  circuit.  Briefly  stated,  calculations  involving  the  load  on  an  equivalent 
circuit  (Thevenin  and  Norton)  are  valid,  but  calculations  concerning  the  source  are  not. 

Battelle' s  investigation  of  the  flywheel  electroxucs  shown  in  Figure  3  disclosed 
that  the  correct  operation  of  the  electronic  amplifier  was  not  as  a  variable-voltage 
source  proportional  to  angular  error  with  a. fixed  source  impedance,  but  rather  as  a 
fixed-voltage  source  with  a  source  impedance  that  varied  as  a  function  of  the  angular 
error.  The  equivalent  circuit  describing  the  operation  of  the  electronic  amplifier  as 
simulated  on  the  analog  computer  is  shown  in  Figure  15. 


Rq 


FIGURE  15.  EQUIVALENT  CIRCUIT  OF  AMPLIFIER-PUMP  INTERFACE  SIMULATED 
ON  THE  ANALOG  COMPUTER 

However,  Battelle' s  investigation  showed  that  the  equivalent  circuit  describing  the 
operation  of  the  electronic  amplifier  is  actually  as  shown  in  Figure  16. 
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Eg  »  Const. 


FIGURE  16.  RESULTANT  EQUIVALENT  CIRCUIT  OF  AMPLIFIER-PUMP  INTERFACE 
FROM  BATTELLE'S  INVESTIGATION  OF  THE  ELECTRONICS 


The  analytic  expression  for  R.  (6  )  is  not  obtainable  from  the  electronic  circuitry 

®  *DB 

shown  in  Figure  3,  and  to  Battelle's  knowledge,  it  has  not  been  determined  empirically. 
Battelle  performed  a  brief  investigation  into  possible  functions  R~(d  )  which  would 

®  *DB 

relate  Figures  15  and  16  to  obtain  identical  outputs  for  corresponding  inputs.  It  was 
concluded  that  R„19  )  was  a  nonlinear  function  of  9  and  would  have  to  be  deter- 

®  *DB  ®DB 

mined  empirically. 


A  reasonable  question  at  this  point  is  "If  the  simulation  of  the  electronic  amplifier 
is  not  for  the  actual  operation  of  the  electronic  amplifier,  can  values  of  power -in  truly 
be  obtained"?  The  answer  is  yes,  and  the  contribution  to  power-in  from  the  electronic 
amplifier  is  obtained  as  the  product  of  the  fixed-voltage  source  (Eg  =  0.  3  volt)  and  the 
source  current  to  the  pump  (1): 

Power  in  =  Pj^  =  0.  3  I 

This  product  does  represent  the  amplifier  power- in  contribution  to  the  system  due  to  the 
fact  that  current  from  the  fixed-voltage  source  is  the  current  to  the  pump.  (No  other 
current  flow  paths  exist  for  the  current  from  the  fixed-voltage  source. ) 

As  shown  in  Figures  15  and  16,  the  source  current  will  be  the  sum  of  the  currents 
through  the  fringing  resistance  (R^)  and  the  pump  channel.  The  pump  channel  has  a  re¬ 
sistance  (Rp)  and  an  internal  voltage  generator,  which  is  a  function  of  the  fluid  angular 
velocity.  By  appropriate  manipulation  of  terms  for  a  specific  set  of  parameters  and  a 
specific  input,  the  expression  for  1  reduces  to 


R  R  +R  R,+R^ 
p  e  p  f  r  e 

4.  56  X  lO'^  ABDjRjW 

R  R  +R  R,+R^ 
p  e  p  f  1  e 


where 


and 


1  = 


C,  = 


C,  = 
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A  constant  value,  »  0.6  watt  was  added  to  the  product  of  Egl  to  represent  the 

power  losses  in  the  electronic  inverter  and  power  supplies,  shown  in  Figure  2,  to  give 
the  total  system  power  input.  The  inverter  and  power-supply  losses  were  obtained  from 
a  G.  C.  report. (*^) 

One  cbsutge  was  necessary  in  the  system  model  for  obtaining  peak  efficiency.  This 
was  the  elimination  of  the  saturation  limiter.  The  excitation  signal  was  a  two-degree 
position  step,  which  would  normally  saturate  the  system  and  bury  the  effects  of  the 
parameter  variations.  The  two-degree  input  was  considered  nece&sary  to  obtain  a  peak 
efficiency  value  at  near-maximum  driving  conditions  for  comparison  with  G.  E.  's 
estimated  value  of  4  to  5  percent.  A  two-degree  i  iput  also  is  analogous  to  obtaining  a 
servomotor  characteristic  curve  at  rated  operating  voltage. 


Degradation  Factors.  General  Electric  has  done  a  substantial  amount  of  develop¬ 
ing  and  testiilg  <m  the  pump  and  flywheel.  However,  their  emphasis  was  not  on  obtaining 
degradalsoc  factors  and  degradation  rates,  but  rather  on  selection  of  materials,  design 
of  the  pump;^  fabrication  of  parts,  reduction  of  weight,  and  optimisation  of  certain 
parameters.  In  iieu  ^  anticipated  data  from  General  Electric,  Battelle  investi¬ 

gated  the  system  to  try  to  obtain  degradation  factors  for  the  system  parameters.  A 
considerable  axnouiU  of  time  was  expended  in  this  effort,  the  result  of  which  was  a 
suggested  test  plan  for  determining  the  degradation  factors  and  their  time  rates.  The 
necessity  for  writing  a  test  plan  points  to  the  lack  of  information  available  regarding 
degradation  factors  for  this  type  of  system. 

Battellc's  Experimental  Physics  Division  was  contacted  to  investigate  possible 
drift-type  failures  in  the  pump  suid  tubing  section,  and  to  make  suggestions  for  possible 
degradation  factors  and  their  time  rates.  The  investigation  showed  that  a  strong  possi¬ 
bility  existed  for  an  increase  in  the  friction-factor  constant  over  time  due  to  chemical 
reaction  between  the  mercury  and  pump  and  tubing  and  due  to  mass  transfer  between 
dissimilar  materials  because  of  temperature  differences  around  the  flow  loop.  Past  ex¬ 
perience  of  the  Experimental  Physics  Division  has  shown  that  this  phenomenon  takes 
place  in  all  systems  of  this  type,  with  rates  of  reactions  depending  on  the  particular 
liquid  metal  and  its  container  material.  It  was  stated  that  a  noticeable,  but  not  large, 
amount  of  material  would  go  into  solution  and  the  process  would  probably  take  several 
years.  The  reaction  takes  place  at  a  more  rapid  rate  initially  and  drops  off  as  the 
process  of  solid  dissolving  into  the  liquid  tends  toward  an  equilibrium  state. 

It  has  been  noted  that  the  rate  of  change  of  this  process  is  exponential  in  nature, 
indicating  that  the  functional  expression  for  the  friction  factor  "constant"  is  of  the  form 

fc  =  a-be'®*  ,  (30) 

where  a  is  the  friction  factor  constant  at  equilibrium  and  a-b  is  the  initial  value  of  the 
.'riction  factor  constant.  Nichol(^)  gives  the  initial  value  as 

a-be  =  a  -  b  =  0.  031,  at  t  =  0  . 
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The  values  of  a  and  a  would  have  to  be  determined  experiment;;s>lly.  For  this  work,  a 
and  cr  were  chosen  in  the  following  manner.  Nicf.ol  states  tV,^t  the  friv^tion  factor  eon 
stant  is  a  function  of  the  wetting  of  the  stainless  steel  tubing  by  the  mercury.  If  water 
were  in  the  tubing  rather  than  mercury,  the  friction  factor  would  be  0.  040,  General 
Electric  realized  that  mercury  would  not  wet  stainless  steel  as  water  wo\ild,  and  ex¬ 
perimentally  determined  the  initial  friction-factor  constant  as  0.031,  Battelle  decided 
to  choose  the  final  value  of  the  friction-factor  constant  to  be  equivalent  to  that  of  water: 

a  =  0.  040, 


Battelle  chose  5  years  as  the  time  for  attaining  equilibrium.  A  reasonable  expression 
for  the  friction-factor  "constant"  as  a  function  of  time,  and  the  one  used  in  this  report 
is 

fc  =  0,040  -  0  009  ,  (31) 


where  t  is  in  years.  The  Fanning  friction  factor  used  in  (.he  flow-loss  equation  now 
becomes 


f  = 


0.040  -  0.009  e 


(-t/1. 25) 


0.  16 


(R^) 


0.  16 


The  tolerance  value  onf^  was  assumed  to  be  ^  5  percent  and  was  assumed  to  remain 
fixed  over  time. 


Delco  Radio  Division  of  General  Motors  Corporation,  Kokomo,  Indiana,  was  con¬ 
tacted  to  supply  and  discuss  information  pertaining  to  possible  degradation  factors  in  the 
transistors  that  were  used  in  the  electronic  amplifier.  The  transistors  were  Delco 
Type  7276361  TO-36.  From  conversations  with  Delco  and  from  the  base-current 
characteristic  curve  at  a  constant  collector  current  of  5  amperes  versus  time,  the 
transistor  gain  (h^  ),  and  consequently  the  system  voltage  "constant",  K^,,  was  deter- 
e 

mined  to  vary  with  time.  From  the  above-mentioned  curve,  it  was  noted  that  the  gain 
decreased  linearly  with  time,  but  at  two  different  rates.  Initially,  the  gain  decreased  at 
a  constant  rate  until  about  500  hours  had  elapsed  and  then  decreased  at  a  lower  constant 
rate  for  the  rest  of  the  time.  The  functional  expression  for  Ky  as  a  function  of  time  is 


K  =  first  500  hours 

V  '  f\  1 


where  t  is  in  hours  and 
volts 


K„  =0.355 


degree 


K-,  =  (K,,  -  500  K  )  -  K_t  after  first  500  hours, 

V  Vq  i  c 


K  -  6.  39  X  10  volts 
1  degree-hour 


K  ••  ^ ^  ^ Q  ^  volts 
2  degre»"hour 


It  was  also  noted  that  the  tolerance  value  of  varied  over  time  in  a  similar  fashion  as 
with  the  following  functional  expression: 
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Reaulti 

The  partial  derivative*  of  the  two  performance  characteristics,  settling  time  and 
open -loop  peak  efficiency,  were  obtained  from  analog  computer  data.  The  analog  data 
were  generated  by  varying  each  parameter  in  predetermined  increments  from  its 
nominal  value  with  all  other  parameters  at  their  nominal  values,  and  exciting  the  system 
with  fixed  inputs.  Each  analog  run  gave  one  point  of  a  "performance  characteristic 
variation"  versus  "change  of  parameter"  plot.  Several  runs  gave  points ' through  which 
a  curve  could  be  drawn.  The  slope  of  this  curve  was  considered  the  partial  derivative 
of  the  performance  characteristic  with  respect  to  the  parameter  considered. 

Settling  Time.  The  time  for  the  system  to  settle  to  within  0.  05  degree  at  system 
lifetime  t  s  0  for  an  input  step  excitation  of  2.  5  degrees  was  13.  25  seconds.  System 
output  and  error  curves  versus  time  are  shown  in  Figure  17.  The  characteristic  loss 
curves  for  the  flywheel  are  shown  in  Figure  18.  Figure  18  differs  from  Nichol's  plot  in 
Figure  7  in  the  values  of  diffuser  and  back  emf  losses.  This  inconsistency  is  due  to 
Battelle's  changing  the  model  of  the  system  by  increasing  the  diffuser  efficiency  and  re¬ 
placing  the  current  source  by  a  voltage  source.  Ths  system  changes  did  net  cause  any 
noticeable  change  in  system  response,  which  can  be  noted  from  Figure  11.  Table  1  lists 
the  parameters,  their  tolerances,  partials,  and  variances,  amd  gives  the  evaluation  of 
the  variance  equation  and  a  breakdown  of  the  variance  in  percent.  The  results  of  the 
evaluation  of  the  variance  equation, shows  that  the  tolerance  value  on  the  settling-time 
characteristic  is  ,  which  is  equal  to  ±1.  08  seconds.  The  settling  time  can  be  ex¬ 
pressed  as  ^8 

t-  =  13.  25  ^  1.  08  seconds 

at  the  time  vehicle  operation  initially  starts.  A.s  can  be  seen  from  the  normalized 

^s  — 

partial  derivatives  P.  in  Table  I,  the  parameters  to  which  the  system  settling  time 

is  most  sensitive  are  tubing  diameter,  D^,  magnetic -flux  density,  B,  pump  channel 
height,  h,  ampliiier  output  impedance,  Rg,  and  the  amplifier  gain,  Ky.  The  parameter 


FIGURE  17.  SYSTEM  OUTPUT  AND  ERROR  VERSUS  TIME,  SYSTEM 
LIFETIME  0 


FIGURE  18.  CHARACTERISTIC  LOSS  CURVES,  SYSTEM  Lif-'ETIME  0 


values  with  the  larger  tolerance  values  are  amplifier  gain,  K^,  amplifier  output  im¬ 
pedance,  Rg,  friction  factor  constant,  f^,  and  voltage  saturation  level,  Vg,  The  break¬ 
down  of  the  variance  of  settling  time  shows  that  the  amplifier  gain,  K^,,  amplifier  output 
impedance,  Rg,  and  magnetic -flux  density,  B,  arc  the  parameters  that  have  the  strong- 

eot  effect  on  the  settling-time  tolerance  value.  These  results  show  that  decisions  con¬ 
cerning  the  important  parameters  must  take  into  account  both  sensitivity  (partial 
derivatives)  and  tolerance  values. 

A  worst-case  analysis  actually  showed  that  the  settling  lime  would  be  shortened  to 
7.  5  seconds  if  the  parameter  values  were  at  their  extremes;  a  care  that  would  sup¬ 
posedly  lengthen  the  settling  time.  This  apparent  discrepancy  is  considered  later  in  the 
report. 


For  an  arbitrary  upper-limit  oettling  time  of  14,  5  seconds  (and  no  lower  limit),  the 
reliability  estimate  is  computed  by  means  of  the  Moment  method  as  follows: 


_  13.  25  -  (-00) 
"l  ~  0.  360 


_  14.  50  -  13.  25  1,  25 

"u  ~  0.  360  0.  360 


3.47 


Percent  P  =50  Percent  P  =  40. 99 

"l  =  «  =  3.47 

Total  =  99.  99  percent. 

The  probability  of  a  flywheel  having  a  settling  time  of  less  than  14.  5  seconds  is  0.  9999, 
The  reliability  estimate  is  the  probability  value,  0.  9999, 

At  system  lifetime  t  =  500  hours,  the  settling  time  was  13. 85  seconds.  System 
output  and  error  curves  versus  time  arc  <^hown  in  Fi'.ure  19.  Figure  20  gives  the 
characteristic  loss  curves  for  the  flywheel  s,  s.orn.  Table  II  gives  the  same  information 
as  Table  I  but  considers  a  system  lifetime  of  500  hours.  The  tolerance  value  on  the 
settling -time  characteristic  is  ±1.  02  seconds.  The  settling  time  can  be  expressed  as 


tg  =  13.  85  ±  1. 02  seconds. 

500  h:urs 

Tubing  diameter,  D2,  and  magnetic-flux  density,  B,  are  the  two  most  important  param¬ 
eters  affecting  the  sensitivity  of  settling  time  and  are  of  equal  consequence.  The  other- 
important  parameters  affecting  the  sensitivity  of  settling  time  to  parameter  shifts  are 
the  pump  channel  height,  h,  and  the  friction  factor  constant,  f^.  The  breakdown  of  the 
variance  of  settling  time  shows  that  the  parameters  having  the  greatest  effect  on  the 
tolerance  value  for  settling  time  are  cimplifier  gain,  Ky,  magnetic  flux  density,  3, 
amplifier  output  impedance,  Rs,  and  the  friction  factor  constant,  fg.  The  rearrangement 
in  the  parameters  affecting  the  sensitivity  of  settling  time  and  the  change  in  the  break¬ 
down  of  settling -time  variance  point  to  the  fact  that  the  partial  deri\'-atives  are  not  con¬ 
stant  over  the  variation  of  the  parameters. 


A  worst-case  analysis  showed  that  the  settling  time  would  be  8.7  seconds.  The 
reliabilitv  estimate  is  as  shown  on  page  40. 
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FIGURE  19.  SYSTEM  OUTPUT  AND  ERROR  VERSUS  TIME,  SYSTEM  LIFETIME 
500  HOURS 


FIGURE  20.  CHARACTERISTIC  LOSS  CURVES,  SYSTEM  LIFETIME  500  HOURS 
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f  ; 


13.85  -  (-«) 

|\  S  ■  II  ■INI  I  Ml  M 

L  0.  340  " 


14.  5  -  13.  85  0.  65  .  ... 

n  IS  .  -  I  s  s  1.91. 

U  0.340  0.34 


P«rcent  P 


»  50 


Percent  P 


n 


=  47.  19 


U  J!  1. 91 


L  *  oe 

Total  ■  97.  19  percent 

The  reliability  estimate  is  0.9  719. 

At  system  lifetime  t  s  5  years,  the  settling  time  was  11. 75  seconds.  System  out¬ 
put  and  error  curves  versus  time  are  shown  in  Figure  21 .  The  characteristic  loss 
curves  for  the  flywheel  system  are  shown  in  Figure  22.  Table  111  gives  the  same  infor¬ 
mation  as  Tables  1  and  II  but  considers  system  lifetime  of  5  years.  The  tolerance  value 
on  the  settling-time  characteristic  is  2.  56  seconds.  The  settling  time  can  be  expressed 
as 

tg  s  1 1 . 75  2.  56  seconds. 

5  years 

The  parameters  to  which  the  settling  time  is  most  sensitive  are  tubing  diameter,  Dz, 
loop  diameter,  D^,  magnetic -flux  density,  B,  amplifier  gain,  Ky,  and  pump  channel 
height,  h,  and  amplifier  output  impedance,  Rs,  of  equal  consideration.  The  variance 
of  settling-time  breakdown  shows  that  amplifier  gain,  Ky,  and  amplifier  output  im¬ 
pedance,  Rs,  have  the  greatest  effect  on  the  variance. 

A  worst-case  analysis  showed  that  the  settling  time  would  be  15.9  seconds.  The 
reliability  estimate  is  as  follows: 


-  n.7S  -  (-«)  _ 
"l  ■  0.854  "  * 


"u  * 


14.  5  -  11.  75  2.  75 


0.  854 


0.854 


=  3.22 


Percent  P 


n. 


=  50 


Percent  P 


=  49. 99 


L  =  3.  22 


L,  =  00 

Total  a  99.  99  percent 

The  reliability  estimate  is  0.9999. 

An  obvious  difference  between  the  results  of  5-year  operation  and  both  the  initial 
starting  operation  and  500 -hour  operation  is  that  the  nominal  settling  time  at  5  years  has 
decreased.  The  decrease  is  quite  large,  being  greater  than  the  settling-time  tolerance 
values  for  both  system  lifetimes  of  t  =  0  and  t  =  500  hours.  The  three  settling -time  ex¬ 
pressions  are  repeated  here: 


a  13.  25  1  08  seconds 


a  11, 75  2.  56  seconds. 


0  hours 


5  years 


a  13.85  ±1.02  seconds 


500  hours 

The  distribution  curves  foi’  each  of  the  three  settling  times  are  shown  in  Figure  23. 
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System  Output  Error,  degrees 


FIGURE  21.  SYSTEM  OUTPUT  AND  ERROR  VERSUS  TIME,  SYSTEM  LIFETLME 
5  YEARS 


FIGURE  22.  CHARACTERISTIC  LOSS  CURVES,  SYSTEM  UFETIME  5  YEARS 
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FIGURE  *’3.  SETTLING-TIME  FREQUENCY  DISTRIBUTIONS 


The  large  decrease  in  the  nominal  settling-time  value  for  a  system  lifetime  of 
5  years  can  be  explained  by  the  fact  that  the  performance-characteristic  settling  time 
is  not  a  continuous  function  of  the  system  parameters.  The  physical  explanation  of  this 
result  can  be  given  by  considering  the  effect  of  decreasing  the  system  gain  (Ky)  and  in¬ 
creasing  system  damping  (f^.)  in  any  underdamped  closed-loop  system  and  considering 
the  general  definition  of  settling  time. 

A  general  definition  of  settling  time  is  the  time  it  takes  for  the  error  of  a  system 
excited  by  a  step  input  to  settle  writhin  a  given  magnitude.  The  acceptable  error  magni¬ 
tude  used  for  defining  settling  time  effectively  defines  two  boundaries  about  the  zero- 
error  axis.  The  time  at  which  the  error  enters  the  region  between  boundaries  and  re¬ 
mains  in  that  region  thereafter  is  the  settling  time.  If  a  system  is  underdamped,  the 
system  error  will  oscillate  about  zero  error  and  possibly  enter,  leave,  and  re-enter  the 
region  between  the  two  boundaries  several  times  before  settling  between  the  boundaries. 

As  the  system  gain  is  decreased  and  the  system  damping  is  increased,  the  nragni- 
tude  oi  tho  oscillaticn  about  the  zero-error  axis  is  decreased.  Each  cycle  of  oscillation 
about  zero-error  axis  will  approach  zero  error  from  a  positive-error  value  on  one-half 
the  cycle  and  from  a  negative-error  value  on  the  other  half  of  the  cycle.  As  the  gain  is 
decreased  and  damping  is  increased,  the  magnitude  of  the  half-cycle  peak  is  decreased 
until  they  are  enclosed  by  the  boundaries  defining  settling  time.  As  the  half  cycle  on 
which  settling  time  occurs  is  enclosed  by  a  boundary  due  to  decreasing  gain  and  increas¬ 
ing  damping,  the  settling-time  point  "jarnps"  to  the  preceding  half  cycle,  causing  the 
discontinuity  in  settling  time.  The  sitviation  described  here  would  occur  if  only  one  of 
the  two  causes,  decreasing  gain  or  increasing  damping,  took  place. 

Figure  24  demonstrates  the  discontinuity  in  the  settling -time  characteristic  by 
plotting  the  distribution  curves  and  the  line  of  nominal  settling -time  variation  versus 
time.  It  can  be  seen  from  Figure  24  that  only  one  "jump"  occurs  in  the  time  peflod 
being  considered,  and  that  it  occurs  at  approximately  a  system  lifetime  of  1/2  year.  If 
the  time  period  were  extended,  more  "jumps"  would  occur,  signifying  other  half-cycle 
peaks  had  been  enclosed  in  defining  settling -time  boundaries. 

The  results  of  using  settling  time  as  a  measure  of  performance  demonstrates  what 
one  can  expect  from  a  performance  characteristic  that  is  a  discontinuous  function  of  the 
system  parameters.  A  discontinuity  in  the  nominal  value  of  a  performance  characteris¬ 
tic  as  a  function  of  time  may  or  may  not  be  a  serious  problem  in  determining  whether  to 
use  the  performance  characteristic  as  a  measure  of  drift-type  failure.  For  the  case  of 
the  fluid  flywheel,  the  "jump"  in  the  settling -time  value  was  not  an  unfavorable  occur¬ 
rence  since  the  jump  in  settling  time  was  away  from  the  specification  limit  and  was 
faborable  to  operation.  Had  the  jump  increased  settling  time,  or  had  there  been  a  lower 
limit  on  settling  time,  then  the  jump  would  present  a  problem.  Knowledge  of  a  jump  in 
an  important  system  performance  characteristic  is  useful  information  for  design  regard¬ 
less  of  whether  the  jump  is  favorable. 

The  worst-cast  analyses  for  the  first  two  system  lifetimes  considered  were  much 
bettei  than  the  actual  settling  times  obtained  with  the  nominal  parameter  values.  The 
worst-case  result.3  are  considered  much  better  than  their  nominal  value  counterparts 
since  the  worst-case  settling  times  ace  much  lower  than  the  nominal  settling  times. 

This  contradictory  result  is  due  to  a  discontinuous  relationship  between  the  settling¬ 
time  characteristic  of  the  system  and  the  system  parameters.  This  result  could  have 
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been  surniiecd  from  the  data  in  Tables  1  through  111  by  noting  that  the  partial  derivative 
of  settling  time  with  respect  to  amplifier  gain,  Ky,  is  negative  in  sign.  The  negative 
slope  indicates  that,  for  worst-case  considerations,  tVie  Ky  should  be  decreased.  The 
question  can  now  be  asked,  "Why  did  the  settling  time  decrease  if  the  partial  has  a 
negative  slope  and  the  Ky  was  decreased?" 

The  question  is  answered  by  the  curve  shown  in  Figure  25,  which  is  a  plot  of 
settling  time  versus  the  parameter,  Ky.  The  plot  is  general  in  that  no  numerical  values 
appear  for  either  settling  time  or  amplifier  gain.  Figure  25  applies  to  any  underdamped 
system  when  the  system  gain  is  changed.  It  can  be  seen  from  the  plot  that,  other  than 
when  a  jump  in  the  curve  occurs,  the  slope  is  negative.  The  overriding  factor  in 
settling  time  as  a  function  of  system  gain  is  the  discontinuity  and  not  the  partial  deriva¬ 
tive.  The  condition  of  discontinuities  (jumps)  in  the  plot  does  not  continue  as  the  gain  is 
continually  reduced.  The  system  eventually  becomes  overdamped,  and  further  reduction 
in  gain  will  always  increase  the  settling  time. 


FIGURE  25.  PLOT  OF  SETTLING  TIME  VERSUS  SYSTEM  GAIN 


Peak  Efficiency.  At  system  lifetime  t  s  0,  the  open- loop  peak-efficiency  nominal 
value  was  4.  ^'7  percent.  The  characteristic  curves  for  the  flywheel  are  shown  in  Fig¬ 
ure  26.  Table  IV  gives  the  same  information  as  Table  I  but  considers  open- loop  peak 
efficiency  rather  than  settling  time.  The  tolerance  value  on  the  peak-efficiency  charac¬ 
teristic  is  576  percent.  The  peak  efficiency  can  be  expressed  as 

Eq  =  4.  27  ±  0.  576  percent. 

The  most  important  parameters  affecting  peak-efficiency  sensitivity  are  tubing  diam¬ 
eter,  D2,  magnetic -flux  density,  B,  and  pump  channel  height,  h;  of  equal  consequence 
are  amplifier  gain,  Ky,  and  amplifier  output  impedance.  The  breakdown  in  percent  of 
peak  efficiency  variance  shows  that  amplifier  gain,  Ky,  amplifier  output  impedance, 

RS,  and  magnetic-flux  density,  B,  have  the  greatest  effect  on  variance. 
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FIGURE  2o.  SYSTEM  CHARACTERISTIC  CURVES,  SYSTEM  LIFETIME  0 


A  worst-case  analysis  showed  that  the  peak  efficiency  would  be  1, 85  percent. 


The  reliability  estimate  is  as  follows,  where  a  lower  limit  of  3  percent  is  arbi 
trarily  assigned  for  peak  efficiency  and  no  upper  limit  is  assigned: 


4.  27  -  3  1. 27  , 

=r92 


00  -  4.  27 

•”u  ^  0. 192“ 


Percent  P  ^  »  50 

=  6.  62 


Percent  P 


-  50 


*U  =  00 


Total  ~  100 

The  reliability  estimate  is  1,0, 

At  system  lifetime  t  =  500  hours,  the  peak  efficiency  was  3.  98  percent.  Figure  27 
gives  the  characteristic  curves  for  the  flywheel.  Table  V  gives  the  same  information  as 
Table  II,  but  considers  peak  efficiency.  The  tolerance  value  for  peak- efficiency  is 
0.  708  percent.  The  peak  efficiency  can  be  expressed  as  E500  hours  =  3.  98  ±  0.  708 
percent.  The  parameters  to  which  peak  efficiency  is  most  sensitive  are  tubing  diam¬ 
eter,  D2f  rragnetic-flux  density,  B,  and  ptimp  channel  height,  h;  of  equal  consequence 
are  amplifier  gain,  Kyi  and  amplifier  output  impedance.  The  breakdown  of  the  variance 
shows  that  the  greatest  effect  on  peak  efficiency  is  caused  by  amplifier  gain,  Ky,  ampli¬ 
fier  output  impedance,  Rs,  and  magnetic -flux  density. 

A  worst-case  analysis  showed  that  the  peak  efficiency  would  be  1.  65  percent. 

The  reliability  estimate  is  as  follows: 

_  3.  98  -  3  0.  98  00  -  3,  98  _ 

"l  "  0,235  "0,  235''  ”U 0.235  " 

Percent  P  *  50  Percent  P  »  50 

"l,  =  4. 17  =  00 

Total  =100  percent. 

The  reliability  estimate  is  1,00. 

At  j'ystem  lifetime  t  =  5  years,  the  peak  efficiency  was  3.  25  percent.  The 
characteristic  curves  for  the  flywheel  are  shown  in  Figure  28.  Table  VI  gives  the  same 
information  as  Table  III,  but  considers  peak  efficiency.  The  tolerance  value  for  peak 
efficiency  is  ±1.  314  percent.  The  peak  efficiency  can  be  expressed  at  E5  years  = 

3.  25  ±  1. 314  percent.  The  parameters  to  which  the  peak  efficiency  is  .most  sensif.^re 
are  tubing  diameter,  Dg,  pump  channel  height,  h,  magnetic -flux  density,  B,  and  the 
amplifier  gain,  Ky,  The  variance  breakdown  shows  that  the  effect  of  amplifier  gair, 

Ky,  completely  swamps  the  effect  from  all  the  other  parameters. 
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A  worst-case  analysis  showed  that  the  peak  efficiency  would  be  1.  07  percent. 
The  reliability  estinnate  is  as  follows: 


0.  438  0,  438 


00-3.  25 
0.  438 


Percent  P, 


«L  =  0.  572 


a  21.64 


Percent  P_  =  50. 

”U  a  « 


Total  a  71.  64  percent. 

The  reliability  estimate  is  0.  7164. 

The  distribution  curves  for  each  of  the  three  peak  efficiencies  are  shown  in  Fig¬ 
ure  29.  A  plot  of  the  peak-efficiency  distribution  curves  and  the  line  of  nominal  peak 
efficiency  versus  time  is  shown  in  Figure  30.  It  can  be  noted  from  Figure  30  that  the 
nominal  peak  efficiency  decreases  continuously  over  time^  with  no  fluctuations  or  sharp 
breaks.  The  path  followed  by  the  line  of  nominal  peak  efficiency  indicates  a  decrease  in 
efficiency  at  a  rate  that  decreases  with  time.  The  decreasing  rate  results  from  ex¬ 
ponential  changes  in  the  friction-factor  constant,  fc  While  the  nominal  value  for  peak 
efficiency  does  not  fall  below  the  specification  limit  in  5  years,  the  spread  of  the  distri¬ 
bution  increases  so  that  approximately  30  percent  of  the  systems  built  could  be  expected 
to  be  operating  below  the  3  percent  specification. 


Comparison  of  Results.  It  is  possible  that,  when  more  than  one  system  perfor¬ 
mance  characteristic  is  used  to  define  proper  system  operation,  a  change  in  one  or  more 
parameters  will  improve  system  operation  when  one  performance  characteristic  is  con¬ 
sidered,  and  will  degrade  operation  when  one  or  more  of  the  other  performance  charac¬ 
teristics  are  considered.  The  situation  of  a  parameter  conflict  between  performance 
characteristics  occurred  in  this  investigation.  By  comparing  settling-time  and  peak- 
performance  partial  derivative  values  in  Tables  I  through  VI,  it  can  be  noted  that 
changes  in  the  loop  diameter,  Di,  pump  channel  width,  w,  electrode  length,  L,  and  the 
friction  factor  constant,  f^,  that  reduce  (increase)  settling  time  correspondingly  reduce 
(increase)  peak  efficiency.  The  fixing  of  the  design  values,  when  parameter  conflict 
exists,  would  have  to  be  made  on  some  basis  other  than  the  two  system  characteristics 
considered  in  this  investigation. 

Table  VII  gives  the  results  of  the  moment  method  for  system  lifetimes  considered, 
and  includes  the  worst-case  analysis  vsdues. 


52 


tiiiow'uj-jaMOd 

m 

in  N 

_ J _ L_ 


SUDM*  inp- 

m  Q  in 

io  CO 

odd 
j _ 1 _ I 


iU9pjed*Apuap!^;3 


53 


FIGURE  28.  SYSTEM  CHARACTERISTICS  CURVES,  SYSTEM  LIFETIME  &  YEARS 


FIGURE  30.  OPEN- LOOP  PEAK- EFFICIENCY  FREQUENCY  DISTRIBUTIONS,  TIME  RELATIONSHIPS 


TABLE  VII,  RESULTS  OF  MOMENT  METHOD 


Characteriitic 

Synem 

Lifetime 

trAi 

ffAj 

“L 

"U 

Reliability 

Wortt  Caie 

Settling  time,  tj 

0 

13, 2S  aec 

0. 1297 

0  360 

eo 

3.47 

0.9999 

7.S  fee 

SOOhr 

13. 8S  fee 

0. 11S3 

0,340 

tc 

1.91 

0.9719 

8.7  tec 

5yt 

11. 7S 

0. 729S 

0. 854 

eo 

3.22 

0.9999 

15. 9  tec 

Peak  efficiency, 

0 

4. 27<^ 

0. 0369 

0.192 

6.82 

«0 

1.000 

1.90% 

E 

SOOhr 

3.d8<^ 

0.0SS3 

0.23S 

4.17 

00 

l.OCO 

1.651fc 

Syr 

3.2S<^ 

0. 1914 

0.438 

0.572 

00 

0.7164 

1.07% 

Conclusions 


The  Moment  method  for  predicting  drift-type  failures  in  a  system  requires  infor¬ 
mation  about  the  system  that  was  only  qualitatively  known  about  the  fluid-flywheel  atti¬ 
tude  controller  investigated.  The  lack  of  information  concerning  time  degradation  of 
parameters  was  particularly  apparent.  The  qualitative  information  was  sufficient,  how¬ 
ever,  to  make  the  educated  assumptions  needed  for  the  analysis.  The  analysis  showed 
that,  for  the  system  characteristics  chosen  and  the  arbitrarily  assigned  specification 
limits,  the  system  has  a  high  reliability. 

This  investigation  demonstrates  the  ability  of  the  Moment  method  to  provide  in¬ 
sights  to  a  designer  as  to  the  behavior  of  the  system  characteristics  being  considered  as 
measures  of  performance.  It  has  been  shown  that  the  results  of  the  Moment  method 
depend  greatly  on  whether  the  performance  characteristic  being  considered  is  a  con¬ 
tinuous  or  discontinuous  function  of  the  system  parameters.  The  knowledge  of  the  func¬ 
tional  relationship  between  performance  characteristics  and  the  system  parameters  is  a 
great  aid  in  design  and  prediction  of  future  system  operation. 

The  partial  derivatives  necessary  for  the  evaluation  of  the  Propagation  of  Variance 
formula,  and  which  determine  the  performance  characteristics  sensitivity  to  the  system 
parameters,  also  provide  additional  information  when  more  than  one  performance 
characteristic  is  considered.  The  partial  derivatives  indicate  which  parameters  are 
compatible,  and  which  parameter  changes  cause  improvement  in  some  performance 
characteristics  and  degradation  in  others.  That  is,  the  partials  indicate  which  param¬ 
eters,  when  varied,  will  concurrently  improve  or  degrade  the  performance  characteris¬ 
tics  (compatible  parameters)  and  those  parameters  which,  when  varied,  will  improve 
some  performance  characteristics  and  degrade  others  (conflicting  parameters). 

The  results  of  this  investigation  show  that,  with  sxifficient  parameter  information, 
the  future  performance  can  be  predicted  and  reliability  estimates  made.  If  the  informa¬ 
tion  is  not  complete,  but  educated  assumptions  can  be  made,  important  insights  into  sys¬ 
tem  performance  can  be  obtained  by  use  of  the  Moment  method. 


TASK  II.  LIQUID  METAL  NaK-77  SERVO- ACTUATING  SUBSYSTEM 


Introduction 

Th«  objective  of  this  t««k  wak  to  develop  a  mathematical  model  of  a  high- 
tempe^ature  liquid-metal  eervoactuating  aubayatem  and,  uaing  reliability  prediction 
techniquea,  to  make  a  prediction  about  the  drift  reliability  of  the  ayatem. 

The  hardware  for  the  ayatem,  being  developed  by  another  contractor,  waa  to  con- 
aiat  of  a  centrifugal  pump,  a  aervo  control  valve,  and  an  actuator  driving  a  epring- 
loadcd  inertial  load. 

Development  work  on  the  ayatem  componenta,  however,  had  not  been  completed  at 
the  time  of  initiation  of  the  program,  and  waa  not  completed  during  the  program.  The 
reliability  analyaia  waa  limited,  therefore,  to  a  atudy  of  the  actuator  and  pump  aince  the 
aervo-valve  deaign  and  configuration  had  not  been  fully  developed.  The  development  of 
the  ayatem  hardware  concurrent  with  the  reliability  atudy  poaed  a  aecond  limitation: 
there  waa  no  time  to  obtain  degradation  data  on  the  componenta:  only  limited  performance 
datA  were  available. 

Mathematical  modela  of  the  pump  and  actuator  were,  however,  developed  and  pro¬ 
grammed  on  the  analog  computer,  and  partial  derivativea  and  variances  of  the  aignificant 
performance  characteriatica  were  calculated.  A  compariaon  of  the  aenaitivity  of  the 
actuator  performance  to  degradation  in  the  actuator  and  in  the  pump  alao  waa  included  in 
the  atudy.  Supplemental  material  including  the  mathematical  statements  for  modeling 
the  ayatem  componenta  and  containing  information  on  the  physical  characteristics  of  the 
componenta  which  may  be  conaidered  proprietary  by  other  contractora,  haa  been 
delivered  to  the  program  Technical  Monitor. 

The  next  logical  atep  in  a  reliability  analyaia  is  the  prediction  of  time  variation  of 
system  and  component  performance.  Such  an  analyaia  is,  however,  based  on  a  know¬ 
ledge  of  the  degradation  rates  of  the  critical  parameters.  Since  information  waa  not 
available  during  the  current  program,  this  phase  of  the  analyaia  waa  limited  to  identifi¬ 
cation  of  the  wear  areas  that  moat  significantly  affect  the  performance  of  the  pump  and 
actuator. 

Engineering  Activity 

Development  of  the  Mathematical  Model  of  the 
Centrifugal  Liquid-Metal  Pump 

The  liquid-metal  pump  studied  in  this  program  waa  a  ten-stage  axially  balanced 
centrifugal  pump  with  6-bladed  radial  vane  impellers,  designed  for  operation  at  35,  000 
rpm  and  rated  at  3000  psi  at  10  gpm. 

A  cross  section  of  the  pump,  showing  the  flew  and  leakage  paths,  is  given  in 
Figure  31.  The  liquid  metal  flows  through  the  first  five  stages  from  right  to  left,  and 
then  is  transferred  to  the  rear  of  the  pump  where  it  flows  through  the  last  five  stages 
from  left  to  right.  The  axial  forces  on  the  impeller  shaft  are  balanced  by  this  reversal 
in  flow  direction.  The  main  pump  shaft  is  supported  by  three  hydrostatically  pressurised 
Journal  bearings  and  one  radial  thrust  bearing.  The  presauriaation  for  the  forward 
bearing  is  taken  from  the  flow  path  between  the  fifth  and  sixth  stages  (e  and  f)  of  the  pump 
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FIGURE  31.  TEN- STAGE  CENTRIFUGAL  PUMP  SHOWING  FLOW  AND 
LEAKAGE  PATHS  AND  STAGE  IDENTIFICATION 
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and  that  for  the  tear  bearing  is  obtained  from  the  fluid  which  enters  the  sixth  stage  (f). 
The  pressure  for  the  middle  bearing  and  for  the  radial  thrust  bearing  is  taken  from  the 
inlet  to  the  tenth  stage  (j)  of  the  pump. 

Figure  shows  a  cross  section  of  the  impeller  and  diffuser.  After  fluid  passes 
through  the  pump  impeller  (1-2)* **,  it  enters  a  vane-type  diffuser  section  (3-4)  and  then 
travels  down  the  back  side  of  the  stage  housing  (4-5),  through  straightening  vanes  in  the 
housing,  which  remove  the  rotary  motion  of  the  fluid  in  preparation  for  its  entering  the 
impeller  of  the  next  stage. 

The  mathematical  model  was  designed  to  produce  a  curve  of  total  pressure  rise 
across  the  pump  as  a  function  of  net  pump  flow.  It  was  formulated  on  the  assumption  of 
steady- state  operation  and  from  the  classical  Euler  equation  for  the  total  head  rise 
developed  across  the  impeller  of  a  centrifugal  machine. 

Of  course,  a  pump  in  a  servo  control  system  would  not  operate  under  steady-state 
conditions,  unless  a  sufficiently  responsive  accumulator  was  included  in  the  circuit. 
Although  no  such  accumulator  is  indicated  in  the  existing  NaK  control  system,  a  steady- 
state  analysis  of  the  pump  was  performed  because  it  enabled  the  engineering  effort  to 
be  concentrated  on  the  modeling  of  the  various  frictional,  leakage,  and  impulse  losses  in 
the  pump,  which  are  difficult  to  predict  because  of  complex  flow  patterns  in  the  device. 

In  addition,  the  steady- state  an<-ilysis  satisfactorily  identifies  parameters  sensitive  to 
degradation  and  predicts  the  effect  of  changes  in  these  parameters  on  pump  performance. 
A  nonsteady- state  analysis  of  the  puirip  would' involve  the  prediction  of  pressure  and 
velocity  distributions  within  the  impeller  and  was  therefore  beyond  the  scope  of  this  pro¬ 
gram.  The  loss  equations  developed  for  the  steady-state  analysis  can  be  used,  however, 
in  developing  a  mathematical  model  of  the  centrifugal  pump  for  nonsteady- state  operation. 


Development  of  Equations  for  Mathematical  Model.  The  equations  developed  in  this 
analysis  are  shown  in  supplemental  material  delivered  to  the  program  Technical  Moni¬ 
tor,  and  are  based  in  part  on  the  work  of  O.  E.  Balje  published  in  the  ASME  Trans¬ 
actions,  74  (1952),  anu  entitled  "A  Contribution  to  the  Problem  of  Designing  Radial 
Turbomachines".  The  theoretical  head  developed  across  the  pump  impeller  is  ex¬ 
pressed  in  Equation  (43)ve  as  a  function  of  the  flow  through  the  impeller,  Q|,  tho  angular 
velocity  of  the  pump,  w,  the  degree  of  prerotation  of  the  fluid  as  it  enters  the  impeller, 
the  geometry  of  the  impeller,  and  a  slip  factor,  m.  It  is  seen  from  the  nature  of  the 
pump  flow  and  leakage  paths,  given  in  Figure  31,  and  from  the  flow  summary  In  the 
supplemental  material,  that  the  flow  through  the  impeller  is  net  pump  flow  plus  inter¬ 
stage  leakage. 

The  slip  factor  accounts  for  the  failure  of  the  fluid  to  leave  the  impeller  tip  tangent 
to  the  impeller  blades.  The  slip  factor  applies  an  average  correction,  greater  than  1, 
and  approaching  a  value  of  1. 0  for  an  infinite  number  of  blades,  to  the  direction  of  the 
relative  velocity  of  the  fluid  with  respect  to  the  tip  blade  angle. 

The  prerotation  factor,  'A,  must  be  considered  because,  when  the  fluid  enters  the 
impeller  blades,  it  may  have  some  rotary  notion  due  to  its  contact  with  ths  rotating 
shaft  prior  to  its  entry  into  the  impeller  or  bs  the  failure  of  the  straightening  vanes  to 
removs  the  angular  velocity  imparted  to  the  fluid  by  previous  stages.  This  rotary  motion 
reducss  the  total  amount  of  snsrgv  U<at  can  be  transfer  cd  to  the  fluid  by  She  impeller. 

*Nuniei(c*l  idar'lflcKion  of  key  paino  in  the  flow  peih  iiirou|h  the  nags. 

**  tiipiitemomkl  metoildl. 
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Impeller 


Section  of  impeller  and  diffuser 


FIGURE  32.  IMPELLER  AND  DIFFUSER  CONFIGURATION,  WITH  IDENTIFICATION 
OF  BLADE  ANGLES  AND  FLOW  LOCATIONS 
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The  total  output  head  from  each  stage  of  the  pump  is  equal  to  the  head  at  the  inlet 
to  the  pump  stage,  plus  the  theoretical  rise  in  total  head  across  the  impeller,  less  the 
various  frictional  and  impulse  losses  discussed  below. 

There  are  several  types  of  head  losses  in  the  centrifugal  pump.  The  inlet  entrance 
loss  and  the  suction-head  area  loss  [Equations  (44)  and  (45)]  are  simply  velocity-head 
losses  that  result  from  a  change  in  the  cross-sectional  area  of  the  fluid-flow  path  between 
the  pump  inlet  port  and  the  inlet  to  the  impeller  of  the  first  stage. 

For  a  given  angular  velocity,  there  is  only  one  flow  condition  at  which  the  inlet  flow 
will  enter  tangent  to  the  impeller  vanes.  At.  all  other  flow  conditions,  the  impact  of  the 
fluid  on  blade  surface  causes  shock  and  turbulance  loss.  This  loss  is  given  in  Equation 
(46).  It  is  expressed  as  a  function  of  the  flow  conditions,  the  geometry  of  the  pump,  and 
a  shock  loss  coefficient,  X,  which  is  experimentally  determined. 

Equation  (48)  gives  the  general  form  for  frictional  flow  losses  in  the  pump.  The 
flow  path  through  the  nonshrouded  impeller  (^2.2^  ^  logarithmic  spiral  form,  and  the 

friction  occurs  between  the  fluid  and  the  surfaces  of  the  stationary  impeller  housing.  The 
frictional  path  length  in  the  diffuser  and  return  sections  is  taken  as  the  mean  path  length 
through  the  diffuser  and  straightening  vanes.  The  velocity  through  these  sections  is  an 
average  velocity  equal  to  the  square  root  of  one-half  the  sum  of  the  square  of  the  magni¬ 
tudes  of  the  inlet  and  outlet  velocities  at  each  section. 

Equations  for  calculating  the  absolute  fluid  velocities  in  the  diffuser  and  the  return 
sections  are  given  in  Equations  (49),  (50),  and  (51). 

The  diffuser  mixing  loss,  as  expressed  by  Equation  (50),  is  a  turbulance  and 
momentum-exchange  loss  that  occurs  because  of  the  difference  in  the  direction  of  the 
absolute  velocities  of  adjacent  stream  lines  leaving  the  impeller  tip. 

The  diffuser  inlet  impulse  loss,  shown  in  Equation  (51),  is  a  shock  loss,  and  is 
similar  to  the  impulse  loss  that  occurs  at  the  entrance  to  the  impeller.  At  the  diffuser 
inlet,  there  is  only  one  flow  condition  under  which  the  fluid  enters  the  diffuser  section 
tangent  to  the  diffuser  vanes.  At  all  other  flow  conditions,  shock  and  turbulance  losses 
’•esult.  An  empirical  shock  loss  factor,  X',  is  used  to  match  the  analog  performance 
with  experimentally  determined  pressure  capacity  curves. 

The  interstage  leakage  flow  path  consists  of  axial  laminar  flow  around  the  pump 
shaft  and  then  a  small  length  of  outward  radial  flow  between  the  end  of  the  impeller  hub 
and  the  pump  housing.  Some  viscous  pumping  action  takes  place  in  this  latter  area; 
however,  the  contribution  of  this  action  to  the  total  leakage  is  much  smaller  than  that  due 
to  the  net  pressure  drop  across  the  total  flow  path,  and  was  therefore  neglected  in  the 
calculations.  The  magnitude  of  the  interstage  leakage  flow  is  given  by  Eciuatiun  (.54) 

Equations  (55)  through  (59)  give  the  le^akage  flow  through  the  three  journal  bearings 
and  the  thrust  bearing.  The  pressures  used  in  calculating  these  leakages  and  bearing 
flow  rates  are  obtained  by  subtracting  the  velocity  head,  (V^/2g),  from  the  total  head  at 
the  point  -.i  question.  Front  seal  leakage,  q^'p,  was  assumed  to  be  zero  in  this  study. 

Consideration  of  Pump  Efficiency.  A  second  important  aspect  of  pump  performance 
IS  efficiercy,  which  is  defined  as  the  ratio  of  the  net  increase  in  hydraulic  horsepower 
across  ;h'  pump  to  the  mechanical  power  put  into  the  pump  shaft.  Although  this  asp’"* 
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of  pump  performance  was  not  apeclfically  modeled  in  thle  study,  the  various  factors  that 
contribute  to  efficiency  losses  in  the  current  >^ump  design  are  of  significance  and  are, 
therefore,  discussed  in  this  section. 

Shaft  input  power  is  required  to: 

(1)  Impact  energy  to  the  fluid  in  the  impeller  of  each  state 

(2)  Overcome  viscous  and  frictional  losses  in  bearings,  impeller,  and  seal 

(3)  Overcome  recirculating  flow  lossea  around  the  ends  of  the  impeller  blades. 

The  fluid  receives  all  of  its  energy  while  passing  through  the  pump  impeller,  and 
the  total  shaft  power  imparted  to  the  fluid  is  equal  to  the  theoretical  power  requirements 
derived  from  the  Euler  equation  as  given  in  Equation  (81).  An  equation  for  determining 
the  power  requirements  for  nonsteady^ state  flow  operation  is  derived.  The  power 
delivered  to  the  fluid  is  not  completely  recovered  at  the  output  of  the  pump.  A  portion 
of  it  is  lost  in  overcoming  frictional  and  shock  losses  in  the  impeller  and  pump  pass¬ 
ages  and  in  throttling  losses,  which  occur  in  the  interstage  leakage  and  bearing  flow 
paths. 

The  power  required  to  overcome  viscous  drag  losses  in  the  journal  bearings  and 
interstage  clearance  spaces  and  the  viscous  losses  in  the  radial  thrust  bearing  can  be 
calculated  from  Equations  (82)  and  (83). 

The  third  type  of  power  loss  results  from  the  leakage  that  occurs  around  the  ends 
of  the  impeller  blades.  The  pressure  on  the  forward  surface  of  the  blade  is  always 
greater  than  that  on  the  back  surface,  and  the  resulting  pressure  differential  creates  a 
tangential  leakage  flow  through  the  clearance  space  between  the  edges  of  the  blade  and  tht 
sutionary  housing.  This  flow  produces  a  viscous  drag  on  the  ends  of  the  blade  and  also 
a  throttling  lose  as  the  fluid  passes  from  a  region  of  high  to  low  pressure. 

The  magnitc.de  of  this  flow  is  difficult  to  calculate  because  a  knowledge  of  the 
pressure  velocity  distributions  within  the  impeller  is  required.  Several  methods  for  " 
obtaining  approxin.ationo  of  this  distribution  are  contained  in  the  literature  and  can  be 
used  as  a  basis  for  predicting  the  pressure  drop  across  the  faces  of  the  Impeller.* 

A  simpler  method  for  obtaining  a  first  approximation  of  the  pressure  differential 
across  the  blade  surfaces  is  given  in  the  supplemental  material.  This  approach  is  based 
on  a  knowledge  of  the  general  shape  of  the  pressure  differential  along  the  blade. 

The  operating  efficiency  of  the  centrifugal  pump  is  equal  to  the  net  increase  in 
hydraulic  horsepower  across  the  pump  divided  by  the  total  input  power  requirements 
discussed  above.  The  net  hydraulic-pamp  horsepower  is  given  by  the  expression! 

(Hp  -  Hs)Qp7. 


Programming  and  Operation  of  the  Analog  Compu;,er 
Simulation  of  the  Liquid-Metal  Pump 

As  a  first  attempt  to  prograi.\  the  mathematical  model  of  the  centrifugal  pump  on 
the  analog  computer,  a  wiring  or  "patch**  diagram  wax  designed  directly  from  the  equa¬ 
tions.  The  program  was  designed  to  solve  for  the  pressure-capacity  performance  of  all 

••Msttiod  fof  Oetigiu  of  Prunp  Impellert  Lblng  a  H^h-«pee<f  Digital  Computet" .  NASA  Technical  Note  0-1569. 
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ten  stages  and  to  account  for  the  dissimilarities  among  the  stages  due  to  differences  in 
leakage  flows.  The  result  of  this  first  effort  was  a  program  approximately  four  times 
too  large  to  be  programmed  on  the  Battelle  analog  computer.  The  primary  difficulty 
was  the  very  large  number  of  variable  multipliers  needed  (about  80  channels)  to  simu¬ 
late  the  entire  pu.np. 

Subsequent  hand  calculations,  however,  showed  that  the  differences  in  the  pre- 
stagw  pressure  rises  that  are  caused  by  the  leakage  flow  variations  among  the  stages 
are  negligible  compared  to  the  very  high  static  pressures  present.  Interstage  leakage 
is,  therefore,  identical  for  all  ten  stages.  It  was  thus  possible  to  closely  simulate  the 
performance  of  the  entire  pump  by  modeling  only  the  first  four  stages  on  the  computer. 
This  required  only  20  channels  of  variable  multiplication.  In  addition,  the  first  four 
stages  are  identical  in  all  respects,  so  that  the  final  model  amounted  to  a  simulation  of 
only  one  stage. 

The  vriring  diagram  and  potentiometer  setting  sheets  are  presented  in  the  supple¬ 
mental  material.  The  potentiometer  settings  shown  are  for  the  nominal  "design  point" 
dimensions  and  performance  of  the  pump.  When  the  operation  of  the  program  on  the 
computer  was  started,  reasonable  values  of  the  shock  loss  coefficient  and  slip  factor, 

X  and  m,  were  arbitrarily  chosen. 

There  were  initially  no  integrators  in  the  program,  because  there  were  no  dynamic 
effects  modeled  and  the  calculations  involved  only  algebraic  operations.  Often,  when 
closed-loop  algebraic  circuitry  such  as  this  is  programmed  on  the  analog  computer, 
instability  occurs.  This  was  true  for  the  pump  model,  and  an  integrator  was  ttterefore 
added  to  the  circuit.  The  integrator  provided  sufficient  filtering  effect  and  was  success¬ 
ful  in  maintaining  stability  of  the  model  without  affecting  the  accuracy  or  performance  of 
the  simulation. 


Simulation  of  Experimental  Data.  In  a  contractor's  Technical  Proposal,  "A 
Study  of  Liquid  Metal,  NaK-77  for  Application  in  Flight  Control  Systems",  predicted 
pump  performance  data  for  several  impeller  speeds  were  supplied.  These  data  were 
duplicated  on  the  computer  by  manipulating  the  values  of  X  and  m.  The  value  of  m  re¬ 
quired  was  very  reasonable,  and  for  35,  000-rpm  operation  corresponded  exactly  to  the 
value  calculated  from  Equation  (43)  in  Appendix  I.  The  value  of  X  required  for  the  dif¬ 
fuser  inlet  was,  however,  somewhat  higher  than  anticipated.  Figure  33  is  a  graph  of 
computer  performance  at  35,500  rpm. 

Experimental  test  data  were  supplied  for  operation  of  this  pump  at  lower  speeds. 
In  order  to  simulate  performance  in  these  cases,  m  had  to  be  reduced  slightly  and  the 
interstage  leakage  had  to  be  increased  considerably.  From  discussions  with  other  con¬ 
tractor's  personnel,  it  was  determined  that  the  interstage  clearances  on  the  pump  used 
in  these  tests  had  been  larger  than  the  drawings  specified  and  that  the  clearances  would 
be  decreased  in  subsequent  pump  models  for  the  higher  speed  operation.  Therefore, 
the  clearances  specified  in  the  drawings  were  used  in  the  simulation  program. 

A  curve  of  predicted  net  positive  auction  head  (NPSH)  was  also  given.  During  the 
initial  runs  on  the  analog  computer,  a  very  poor  simulation  of  this  was  obtained.  Upon 
inspection  of  the  program  equations,  it  was  determined  that  the  equation  for  impen-"- 
inlet  shock  loss  as  developed  by  Balje  was  not  directly  applicable  to  the  pump 'design 
under  consideration  because  it  orovided  for  the  existence  of  an  impeller  inlet  shock 
loss  at  n^-flow  conditions.  The  equation  was  herefore  modified  to  have  the  shock  lose 
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FIGURE  33.  COMPUTED  CENTRIFUGAL-PUMP 
PERFORMANCE  AT  36»  500'  RPM 


approach  zero  at  no  flowj-nd  to  rnatth  the  curve  rasul'-.inj,  from  Balje's  equations  at 
higher  flows.  The  equations  for  the  diffuser  inlet  shoi'k  losses  were  satisfactory.  Fig¬ 
ure  34  is  a  graph  of  computed  NPSH.  This  is  lower  thun  tlat  predicted,  but  it  is  be¬ 
lieved  that  the  prediction  was  based  on  a  general  formula,  which  probably  w.^uld  not  be 
as  accurate  as  the  specific  program  used  for  this  project.  However,  the  NPSH  pre¬ 
dicted  may  include  some  idditional  local  turbulence,  anu  therefore  may  be  justifiably 
conservative.  The  partial  .’erivatives  of  particular  interest  at  this  time  are  not,  how¬ 
ever,  related  to  NPSH  and,  b*  •a'rse  a  reasonably  good  simulation  of  predicted  and 
experimental  preii  sure- capacity  •  irformance  had  been  obtained,  it  appeared  that  the 
simulation  program  could  be  usee  with  confidence  to  compute  partial  derivatives  of 
performance. 


Analysis  of  Pump  Reliabi*  '•'v 

A  study  of  the  configuration  and  matnematical  model  of  the  centrifugal  pump  reveals 
that  pump  performance,  in  terms  of  th*  pressure-capacity  curve,  is  most  dependent  on 
the  magnitude  of  t’  i  (1)  \  •  ilier  tip  radius,  (2)  impeller  blade-tip  angle,  and  {3)  the 
interstage  clearance. 

The  sensitivity  of  the  pump  performance  to  changes  in  these  paran eters..  expresaed 
as  partial  derivatives,  forms  he  basis  of  the  study  of  pump  reliability  by  the  Moment 
method  of  reliability  anol'  sis.  These  partial  derivafives,  when  expressed  in  absclute 
terms  (not  normalized),  e  used  in  the  Propagation  of  Variance  fox’mula  to  estimate  the 
expected  variation  m  performance  due  to  tolerance  variations  introduced  in  tire  manu¬ 
facturing  process.  The  normalized  partial  derivatives  fexpressedl  as  a  change  in  per¬ 
formance  per  percent  ^  .ange  in  parameter)  provide  a  measure  of  the  sensitivity  of 
performance  to  changes  or  degradation  of  the  component  parameters.  These  are  useful 
for  identifying  crit;cal  degradation  areas  and  for  redesigning  components  to  improve 
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reliability  and  performance.  When  the  time  rate  of  degradation  of  the  parameters  is 
known,  the  absolute  partial  derivatives  are  used  as  weighting  factors  in  predicting  per¬ 
formance  changes  during  operation. 

Because  quantitative  data  on  the  degradation  of  both  the  NaK  system  pump  and 
actuator  were  not  available^  the  analysis  on  these  two  components  does  not  include  a 
quantitative  estimate  of  performance  variation  with  time.  However,  Task  !II  of  this  re¬ 
port  does  contain  a  detailed  discussion  cf  the  use  of  the  Moment  method  in  predicting 
time- dependent  performance  changes  on  an  f>tectrohydraulic  servo-control  valve,  and  the 
baste  methods  used  in  this  study  could  be  applied  to  the  NaK  system  when  parameter 
degradation  data  become  available. 

Partial  derivatives  of  pump  performance  are  computed  by  varying  the  critical 
parameters  of  the  analog  model  one  at  a  time,  and  observing  the  variations  that  occur  in 
pump  performance.  The  variations  were  obtained  by  changing  the  necessary  potentiom¬ 
eter  settings  on  the  analog  computer.  For  each  parameter  variation,  a  new  head- 
capacity  curve  was  run  on  the  computer  and  plotted  on  the  X-Y  plotter.  Figure  35  shows 
the  head-capacity  curves  for  the  rknge  of  parameter  values  investigated.  Shown  are 
(1)  a  normal  curve,  (2)  curves  for  a  10  percent  increase  and  decrease  in  impeller  tip 
radius,  (3)  curves  for  a  !0  percent  increase  and  decrease  in  interstage  diametral  clear¬ 
ance,  and  (4)  a  curve  for  a  10  percent  decrease  in  impeller  blade-tip  angle. 

Partial  derivatives  were  computed  at  capacities  of  5,  15,  and  25  gpm,  and  the 
curves  of  performance  at  IS  gpm  from  which  the  partial  derivatives  were  calculated  are 
shown  in  Figures  36,  37,  and  38.  Table  VUl  presents  a  summary  of  the  partial  deriva¬ 
tives  of  pressure  at  each  of  the  flow  capacities  with  respect  to  each  of  the  parameters 
investigated.  These  are  "normalised"  partial  derivatives,  computed  on  a  percentage 
basis. 


TABLE  Vm.  SUMMARY  OF  NORMALIZED  PARTIAL  DERIVATIVES 
PREDICTED  BY  ANALOG  COMPUTER 


Partial 

Derivatives 

Value, 

pei/% 

5  gpm 

15  gpm 

25  gpm 

dP/aR2^ 

23.4 

28.7 

35.0 

7.  8 

-1.5 

-6.9 

31.  f 

40.0 

45.0 

The  nominal  values  of  the  parameters  are: 

Rj  •  0.  764  in. 

^2b  ■ 

Sba  "  In 

The  standard  deviations  for  the  three  critical  pump  parameters  are  obtained  from 
the  detailed  drawing  of  the  pump  by  assuming  that  the  specified  tolerance  range  represents 
a  totol  of  six  stondard  deviations  (^a>.  The  deviations  obtoined  are: 
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Pressure  Rise.psi 


Tip  diAmetcr  0,  OOOi^tS  ln» 

Blade  angle  0.  iU9  deg 

luteretage  0. 000033  in. 

clearance 

The  reeulting  variance  in  ontput  preeeure  at  a  nominal  flow  of  15  gpm  ie  then  obtained 
from  the  Propagation  of  Vaxiimce  formula  (Figure  1)  and  ie  computed  to  be  25.  685  pei. 
The  expected  performance  for  a3o  variation  at  15  gpm  is  therefore: 

3050  a  15.  2pai. 


Development  of  Mathematical  Model 
of  Liquid  Metal  Actuator 

Figure  39  is  a  drawing  of  the  liquid-metal  actuator  that  was  studied  in  the  program. 
It  uses  a  metallic-lip-type  rod  seal  and  metal  piston  rings.  The  rod  is  422  stainlees 
steel  with  a  LW-1  flame-plated  carbide  coating.  The  pump  housing  is  422  stainless,  and 
the  liner  and  front  bushings  are  carbide. 

The  performance  and  degradation  studi>;s  of  the  actuator  were  run  on  the  analog 
model  of  a  flow-control  servo  valve  and  linear  actuator,  which  were  developed  during 
Project  8225,  Task  No.  822502,  and  described  in  Techutcal  Documentary  Report  No. 
FDL-TDR-64-30.  The  configuration  of  the  NaK  actuator  differs  from  that  of  the  actuator 
in  chit  original  system  by  having  a  single-ended,  rather  thhu  a  double-ended,  connecting 
rod.  The  response  and  performance  of  the  NaK  system  is  therefore  dependent  on  the 
direction  of  motion  of  the  actuator.  In  order  to  study  and  compare  the  properties  of  the 
actuator  in  both  directions,  however,  the  double-ended  actuator  of  the  original  systetr 
was  utilized,  and  data  were  run  with  the  pressurized  area  of  the  do’'b'‘«.<<  ended  ac.buator 
equal  to  both  the  larger  and  the  smaller  areas  of  the  NaK  actuator.  The  larger  area  then 
represents  motion  to  the  left  in  Figure  39,  and  the  smaller  area  represents  motior.  :o 
the  right. 

The  performance  of  an  actuator  is  tied  to  the  performance  of  the  entire  system  of 
which  it  is  a  part,  and  in  the  NaK  system  the  actuator  will  ultimately  be  installed  in  a 
closed-loop  control  system  with  a  7-gpm  flow-control  valve  with  a  spring-mass  load  as 
described  in  another  contractor's  proposal.  The  predicted  performance  of  the  system 
was  described  as  "flat  to  20  radians  per  second  at  20  percent  stroke  amplitude",  A 
frequency-i'esponse  plot  of  predicted  performance  was  also  presented  in  the  proposal. 
Comparison  of  this  plot  with  other  open-loop-type  data  indicates  that  i*  is  an  open-loop 
-esponse  of  the  actuator  load  only,  without  Che  control  valve  or  fluid  in  the  circuit. 

To  evaluate  performance  and  degradation  se*isitivity,  the  modeled  system  was  set 
up  for  closed-loop  operation  with  proportional  control.  The  gain  was  set  reasonably 
below  the  marginal  stability  point,  to  providv';  good  control  without  overshooting.  Droop 
was  0.  05  percent  at  20  percent  of  maximum  stroke  from  center  position.  No  celay  or 
other  compensation  circuitry  was  uicluded,  even  though  the  system  proposed  uses  a  lag 
network  in  the  control  circuit.  In  addition,  the  par.smeters  of  the  originally  modeled 
servo  valve  were  set  to  produce  a  saturation  at  7  gpm. 


The  syetem  was  cycled  at  20  percent  of  maximum  stroke  amplitude,  and  the  re¬ 
sulting  frequency  response  plot  expressed  as  a  flow-amplitude  ratio  for  both  directions 
of  motion  is  given  in  Figure  40.  It  was  observed  that  the  flow-control  valve  remained  at 
full  stroke  during  most  of  the  cycle  and  that  the  actuator  was  therefore  limited  by  the 
maximum  flow  capacity  of  the  control  valve.  This  effect  occurs  for  large  actuator  mo¬ 
tions  when  the  response  of  the  control  valve  is  much  higher  than  that  of  the  actuator. 

The  valve  then  acts  as  an  on-off  device,  and  the  system  performance  is  nonlinear.  The 
effective  system  gain  is  also  greatly  reduced  at  the  high  error  amplitudes. 


Analysis  of  Actuator  Reliability 

The  performance  of  the  NaK  actuator  can  degrade  as  a  result  of  several  types  of 
failure.  Friction  can  build  up  due  to  misalignment  of  seals  or  to  bearing  damage,  or 
leakage  can  occur  either  through  the  end  shaft  seal  or  across  the  piston- 

If  the  actuator  is  functioning  properly,  seal-end  bearing  friction  and  leakage  of 
NaK  from  fhe  end  shaft  seal  will  be  small.  If  the  shaft  friction  or  end-seal  leakage  does 
begin  to  increase  during  operation,  the  rate  of  increase  will  generally  be  quite  rapid  and 
could  be  considered  as  a  catastrophic  failure.  Wear  in  the  piston  rings,  however,  will 
nurmally  progress  at  some  relatively  moderate  rate  and,  as  long  as  the  seals  and  bear¬ 
ings  are  satisfactory,  degradation  can  be  considered  as  arising  from  piston-ring  wear. 
If,  however,  the  performance  of  the  shaft  seal  or  bearing  begins  to  degenerate,  this 
degradation  will  completely  overshadow  any  effects  of  piston  wear.  Therefore,  across- 
piston  leakage  was  the  only  type  of  failure  of  interest  in  this  drift  reliability  analysis. 

An  across-piston  leakage  rate  of  1.  S  cis  at  500  psi  and  room  temperature,  as 
obtained  from  test  data,  was  used  aa  the  "design  point"  nominal  leakage  and  was  used 
in  determining  the  piston-ring  clearance  for  the  "design  point"  performance  in 
Figure  40. 

To  compute  the  changes  in  performance  that  would  result  from  across-piston  leak¬ 
age,  frequency  response  was  computed  for  the  servoactuator  system  for  each  of  three 
off-design-point  conditions: 

(1)  Low  leakage,  0.  15  cis  at  500  psi 

(2)  High  leakage,  3  cis  at  500  psi 

(3)  High  leakage,  6  cis  at  500  psi. 

Normal  leakage  was  assumed  to  be  1.5  cis  at  500  psi,  which  was  supplied  from  room- 
temperature  tests.  The  leakage  may  be  considerably  higher  at  high  temperature,  thus 
providing  a  different  design- point  performance.  However,  partial  derivatives  taken  for 
the  room-temperature  performance  should  be  useful. 

Figure  41  is  a  graph  cf  the  servoactuator-system  frequency  response  for  the  large 
area  and  for  the  three  off-design- point  leakage  rates  previously  discussed.  By  com¬ 
paring  these  and  data  for  the  smaller  area  with  the  "normal"  responses  of  Figu.  e  40, 
the  normalised  partial-derivative  data  of  Table  IX  were  compiled.  The  partial  deriva¬ 
tives  were  computed  with  respect  to  the  diametral  clearance,  a,  the  effective  value  of 
which  was  computed  to  be  0.  004  in.  at  room  temperature. 


74 


FIGURE  40.  SYSTEM  FREQUENCY  RESPONSE 


,r. .. WH  -■ 


TABliE  IX.  SUMMARY  OF  NORMALIZED  DERIVATIVE  DATA 
OF  ACTUATOR  PERFORMANCE, 


20  Radians/ 

29.  3  Radians/ 

52.  1  Radians/ 

Sec 

Sec 

Sec 

Large  area 

0.0  db/% 

-0. 0081  db/% 

-0.014  db/% 

Small  area 

0.0  db/% 

-0.0085  db/% 

-0. 023  db/% 

Analysis  of  Syatem  Reliability 

After  the  critical  partial  derivatives  cf  the  performance  of  each  o'  the  elements  in 
a  system  have  been  obtained,  it  is  then  necessary  to  determine  which  of  the  elements  are 
the  most  critical  as  far  as  degradation  of  total  system  performance  is  concerned.  In 
this  study  the  relative  sensitivity  of  system  performance  (actuator  frequency  response) 
to  degradation  in  the  pump,  as  compared  with  that  in  the  actuator,  was  investigated. 

This  was  done  by  first  determining  the  absolute  partial  derivatives  of  actuator  fre> 
quency  response  with  respect  to  system  supply  pressure,  dFR/dP.  The  resulting  par- 
tials  for  the  areas  indicated  at  a  nominal  across-piston  leakage  of  1.5  cis  are  as  follows: 


Large  Area 
Small  Area 


20  Radians/ 

Sec 

-0.  00072  db/psi 

0.  0 


29.3  Radians/ 
_ Sec 

-0.00168  db/ps} 

-0.00112 


52.  1  Radians/ 
_ Sec 

-0.00168  db/psi 
-0.00276  db/psi 


The  normalised  partial  derivatives  of  actuator  frequency  response  with  respect  to  the 
critical  pump  parameters  can  then  be  calculated  from  the  following  equations; 


dFR 

dFR 

ap 

aR2% 

"  dP 

bR2% 

&FR 

3fr 

hP 

d/S% 

dP 

a^% 

3fr 

3fr 

ap 

^®ba^ 

"  dP 

^a 

The  -esulting  normalised  partial  derivatives  of  frequency  response  with  respect  to  across- 
piston  leakage,  impeller-tip  radius,  blade-tip  angle,  and  interstage  clearance  at  5-gpm 
valve  glow,  across-piston  leakage  of  1.  5  cis,  and  frequency  ot  52.  1  radians  cer  second 
arc: 
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arR/ac^^^ 

apR/aRi% 

aFR/aab^% 


Large  Area 
-0.  022  db/% 
-0.0593  db/% 
-0.0530  db/% 
-0.00} 3 1  db/% 


Small  Area 
-0.  0562  db/% 
-0.  0647  db/% 
-0.  0871  db/% 
-0.  00214  db/% 


Xt  icv  seen  from  the  above  results  that  degradation  of  the  pump  parameters  has  a 
strong  effect  on  system  performance  as  well  as  the  degradation  in  the  actuator  piston 
rings  on  a  percentage  basis.  This  once  again  holds  true  only  as  long  as  the  front  shaft 
seal  and  bearing  on  the  actuator  are  operable. 


For  an  estimate  of  the  numerical  ef^'ect  of  piston  leakage  on  system  frequency  re- 
aponse,  calculations  are  made  for  an  assumed  piston  equivalent  clearance  tolerance. 
The  system  variance  in  actuator  fj  equency  response  is  obtained  from  the  equation: 


With  a  nominal  piston-ring  clearance  tolerance  of  0.  0004  in.  ,  the  standard  devia¬ 
tion  of  the  across-pistun  clearance  is  assumed  to  be  0.  00013  in.  and  the  system  pressure 
standard  deviation  is  assumed  to  be  25.63  psi.  The  resulting  deviation  of  the  system 
frequency  response  in  terms  of  actuator  and  pump  manufacturing  tolerances  is  then 
0.  0436  db  and  0,  0728  db,  and  the  expected  tolerance  in  system  performance  is  -7.  Z  ± 

C.  13  db  and  -5.  4  i;  0.  22  db  for  the  large  and  small  areas,  respectively,  at  52.  1  radians/ 
sec. 


The  ralctilations  show  also  that  the  assumed  tolerance  in  the  actuator  would  con- 
tribute  approximately  2.  6  percent  of  the  variance  when  the  large  area  is  active  and  6.  2 
percent  when  the  small  area  is  active. 


Conclusions 


One  obvious  but  important  result  of  Task  11  of  this  program  was  a  demor  atrrtion  of 
the  ability  of  the  Moment  method  reliability  analysis  technique  in  conjunction  with  mathe- 
:  latical  and  analog  models  to  predict  the  variation  that  can  be  expected  in  the  performance 
of  a  Lervo  control  system  and  in  the  individual  system  components,  and  to  aid  in  the 
eva’  ation  of  the  sensitivity  of  th^se  performance  characteristics  to  gradual  degradations 
in  the  physical  parameters  of  the  system. 

The  analysis  has  shown  that  the  variance  in  system  performance,  which  in  this 
case  is  taken  to  be  actuator  frequency  response,  is  due  to  tolerances  in  the  effective 
clearance  around  the  actuator  piston  rings  and  tolerances  in  pump  parameters.  System 
performance  is  sensitive  to  degradationo  in  the  tip  radius  of  the  pump  impeller  and  in  the 
angle  of  the  blade  tip;  these  parameters  produce  about  twice  the  change  in  system  per¬ 
formance  as  an  equal  per  .entage  increaoe  in  the  piston-ring  clearance.  Leakage  betweer 
stages  of  the  centrifugal  pump  has  only  a  sJ'ght  influence  on  the  pump  and  system 
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performance.  It  v/ould,  however,  have  more  significance  if  pump  efficiency  were  being 
considered  as  the  pump  performance  characteristic. 

While  there  were  no  data  available  on  the  extent  of  parameter  degradation,  a  study 
of  the  pump  and  actuator  configurations  leads  to  the  conclusion  that  changes  in  blade->tip 
radius  and  angle  due  to  liquid-metal  erosion  would  be  very  slight.  Also,  since  inter¬ 
stage  leakage  does  not  significantly  affect  frequency  response,  it  could  be  expected  that 
wear  in  the  piston  rings  of  the  linear  actuator  would  be  the  most  critical  in  the  system. 
This  assumes  that  the  design  of  the  actuator  rod  seal  and  bearing  is  fully  developed,  so 
that  their  performance  is  satisfactory  over  the  required  life  of  the  system.  If  the  design 
is  not  fully  developed,  then  this  area  may  be  the  most  critical  because  degradation  in  the 
bearing  and  seal  can  lead  to  a  "catastrophic"  type  of  failure. 

The  above  reliability  observations  provide  useful  design  guidance  during  the  devel¬ 
opment  stages  of  new  equipment.  Decisions  to  reduce  or  increase  manufacturing  toler¬ 
ances  on  parts  that  affect  the  various  system  parameters  can  be  made  on  the  basis  of  the 
degree  of  contribution  of  parameter  tolerances  to  the  total  variance  in  system  perfor¬ 
mance.  The  normalized  partial  derivatives  of  system  performance  indicate  the  sensi¬ 
tivity  of  performance  to  changes  in  system  parameters;  when  these  paxtials  are  con¬ 
sidered  along  with  the  relative  degradation  rates  of  tlie  parameters,  critical  wear  areas 
can  be  pinpointed.  For  improvements  in  reliability,  these  critical  areas  should  be 
studied  first  to  see  where  changes  in  design  configuration  or  in  material  selection  could 
produce  longer  component  life. 
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TASK  in.  REDESIGN  OF  AN  ELECTROHYDRAUUC  SERVO  VALVE 


This  acction  of  tb«  raport  cover*  the  application  of  the  Moment  method  of  reliability 
prediction  to  a  two-etage  electrohydraulic  servo  valve.  Task  HI  objectives  were  in  two 
parts: 

(1)  The  objective  of  the  first  part  was  to  use  the  Moment  method  and  the  mathe¬ 
matical  model  previously  developed  for  a  two-stage  servo  valve  to  improve 
the  design  of  the  valve  from  the  standpoint  of  reliability. 

(2)  The  second  part  consisted  of  analytically  predicting  the  drift  reliability  of  the 
improved  servo  valve  relative  to  that  of  the  original  design  and  conducting 
life  tests  to  verify  the  predicted  relative  reliability. 


Engineering  Activity 


The  engineering  activity  for  Task  m  consisted  primarily  of  three  subta«ks.  First, 
the  Moment  method  of  statistical  analysis  was  applied  to  a  two-stage  servo  valve. 

Second,  two  modified  and  two  unmodified  servo  valves  were  life  tested.  The  thir'c’.  part 
consisted  of  analysis  of  the  data  gathered  during  the  life  testing  in  order  to  show  the 
degree  of  correlation  between  the  predicted  and  actual  improvements  in  reliability. 


Valve  Redesign 

This  section  describes  the  procedure  and  results  of  the  effort  to  identify  parameter 
changes  that  would  improve  the  reliability  of  an  electrohydraulic  servo  valve.  The 
methods  used  for  this  subtask  were  those  associated  with  the  Moment  method  of  relia¬ 
bility  prediction.  The  Moment  method  is  a  statistical  technique  used  for  theoretically 
predicting  certain  attributes  of  the  distribution  of  performance  characteristics  of  a  device 
or  system.  Use  of  the  Moment  method  requires  a  mathematical  model  of  the  device  or 
system  under  consideration,  and  also  requires  data  on  the  actual  or  assumed  distributions 
of  the  parameters  governing  the  performance  characteristics  for  which  the  distribution 
attributes  are  to  be  derived. 

The  first  step  in  the  procedure  was  the  recognition  and  description  of  possible 
mechanisms  by  which  degradation  can  occur  in  the  servo  valve  and  the  identification  of 
those  areas  and  the  specific  parameters  in  the  valve  that  are  subject  to  the  effects  of 
degradation.  It  was  assumed,  in  this  study^  that  the  valve  design  is  at  anadv  mced  stage 
of  development  where  gradual,  drifHtype  failures  —  not  sudden,  catastrophic  failures  — 
are  of  concern. 

Mechanics  of  Degradation.  A  list  of  the  degradation  modes  and  the  affected 
parameters  considered  in  the  analysis  is  given  in  Table  X. 


TABLE  X.  MODES  OF  DEGRADATION  AND  AFFECTED  PARAMETERS 


Modes  of  Degradation 


A.  Wear 

B.  Dimensional  Instability 

C.  Parameter  Instability 

D.  Component  Shift 
£.  Human  Error 


II.  Identification  of  Affected  Parameters 

A.  Wear 


1 .  Metering  Edge  Rounding  and  Port  Coefficient 

2.  Spool-Bushing  Clearance 

3.  Flapper  Wear 

4.  Nozzle  Coefficient 

5.  Orifice  Coefficient 

6.  Feedback-Spring  Ball 

7.  Spool  Drag  Coefficient 

B.  Dimensional  Instability 

1 .  Spool  and  Bushing  Diameters 

2.  Nozzle  Diameter 

3.  Orifice  Diameter 

4.  Armature  Thickness  and  Pole-Piece  Height 

5.  Flapper  Thickness 

6.  Spool  Length 

7.  Nozzle  Location 

8.  Armature-Assembly  Location 

9.  Bushing  Location 

10.  Magnet  and  Pole-Piece  Location 

C.  Electromagnetic  Instability 

1 .  Magnet  and  Pole- Piece  Permeability 

2.  Coil  Resistance 

3.  Magnet  Strength 

4.  Coil  Inductance 


D.  Component  Shift 

1.  Nozzle  Location 

2.  Ar.nature-Assembly  Location 

3 .  Bushing  Location 

4.  Magnet  and  Pole- Piece  Location 

E.  Human  Error 


1.  Affects  the  Degree  of  Severity  of  A,  B,  C,  and  O,  Above 
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W«>r  is  dcfinsd  *s  ths  chang«  in  ths  shape  of  a  valve  component  due  to  erosion  by 
ddfitAaSSaAts  in  tiie  hydraulic  fluid,  and  as  a  change  in  the  surface  condition  of  a  com* 
poaeat  due  to  abrasion,  scoring,  or  impingement  by  solid  contaminant  particles. 

Dtmeneiooal  iastabiltty  is  defined  as  the  change  in  physical  dimensions  that  occurs 
as  a  fiiactioa  ol  time,  stress,  and  temperature  and  is  usually  associated  with  creep  or 
grain  growth. 

Electrwnt^gnetic  instability  is  defined  as  the  change  that  takes  place  in  the  electri¬ 
cal  and  magnetic  characteristics  of  ?be  valve  over  a  period  of  time  as  a  result  of  the 
thermal  and  magnetic  environment  in  which  the  valve  is  located. 

Component  shift  is  de'  .ne.i  as  the  changes  that  occur  in  the  relative  positions  of 
valve  components  due  to  diffeT:.'ences  in  thermal-expansion  rates  and  to  acceleration 
forces  produced  by  mechanical  shock  or  vibration. 

Human  error  is  an  indirect  mechanism  of  degradation,  as  it  describes  the  increases 
that  can  occur  in  the  severity  of  the  above  four  modes  of  degradation  due  to  human 
errors  introduced  in  the  manufacturing  or  assembly  process  that  go  undetected  in  the 
final  chsekeut: 


Selection  of  Most  Significant  Mechaniams  of  Degradation.  The  next  step  in  the  pro¬ 
gram  was  a  study  of  the  modes  of  degradation  and  the  identification  of  those  modes  with 
the  corresponding  parameters  that  significantly  contribute  to  the  degradation  of  valve 
performance.  This  evaluation  was  based  on  Battelle's  knowledge  and  experience  from 
past  programs  and  on  a  study  of  applicable  literature. 

Th }  study  revealed  that  wear  is  the  primary  cause  of  valve-performance  degrada¬ 
tion  and  hat  this  manifests  itself  in: 

(1)  An  increase  in  the  discherge  flow  coefficient  of  the  orifices  and  nozzles 

(2)  Deformation  and  pitting  of  the  flaoper  surface 

(3)  Rounding  of  the  spool  metrving  edge 

(4)  An  increase  in  the  ^adial  spool  clearance 

(5)  An  increase  in  the  rculomb  friction  force  on  the  spool  due  to  scoring. 

The  rounding  of  the  spool  metering  edge  was  modeled  as  an  increase  in  valve  underlap 
and  in  spool  radial  clearance. 

Experimental  and  analytical  work  to  date  indicates  that  serious  problems  are  not  to 
be  expected  in  the  valve  because  of  electromagntiic  and  dimensional  instability.  Valve 
materials,  especially  those  machined  to  close  tolerance?,  are  chosen  for  their  high 
stability.  The  normal  environment  in  which  these  valves  are  used  is  generally  not 
severe  enough  to  came  excessive  degradation  of  tbe  electromagnetic  or  mechanical 
properties. 

It  was  also  concluded  that  there  are  no  serious  effects  from  thermal  cycling  and 
mechanical  and  thermal  shuck,  because  the  valve  is  stabilized  for  mechanical  vibrations 
and  thermal  cycling  during  the  final  stages  of  manufacture.  In  addition,  the  valves  in 
normal  service  arc  t.os  s.  ibjectea  to  c  shock  or  vibration. 
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The  problem  of  human  error  ia,  of  courae,  alwaya  preaent,  but  it  was  concluded 
that  the  occurrence  of  thia  in  a  highly  developed  and  carefully  inspected  valve¬ 
manufacturing  process  is  small,  and  that  it  would  be  difficult  and  beyond  the  scope  of 
this  program  to  attempt  to  define  the  rate  of  occurrence  and  the  effects  of  human  error 
in  any  quantitative  way. 

The  objective  of  Task  in  included  the  investigation  of  effects  of  design  changes  in 
the  critical  parameters  on  the  reliability  of  the  servo  valve.  Because  very  little  informa¬ 
tion  is  available  on  the  time  rate  of  degradation  of  these  parameters,  the  technique  is 
limited  to  predicting  relative  rather  than  absolute  reliability.  It  is  noted  that  this  re¬ 
striction  on  the  application  of  the  technique  is  set,  not  by  limitations  in  the  mathematical 
and  analog  models,  but  rather  by  limitations  in  the  ability  to  predict  rates  of  degradation. 
A  parameter  as  it  appears  in  the  model  can  be  manipulated  if  there  is  some  knowledge 
of  how  it  should  be  manipulated,  i.  e. ,  if  the  rate  of  change  of  the  parameter  with,  say, 
flow,  is  available,  this  information  can  be  corrected  to  a  time  basis  and  used  in  the 
reliability  analysis. 


Consideration  of  Possible  Design  Changes  for  Improved  Reliability.  Improvement 
in  reliability  can  generally  be  achieved  in  three  ways.  The  first  is  to  reduce  the  sensi¬ 
tivity  of  the  valve-performance  characteristics  to  degradation  of  the  parameters  and  to 
variations  introduced  in  manufacture,  i.  e. ,  to  reduce  the  partial  derivatives.  The  sec¬ 
ond  is  to  improve  the  nominal  valve  performance  to  create  a  greater  tolerance  for 
degradation,  and  the  third  is  to  make  design  and  material  changes  that  would  reduce  the 
rate  of  parameter  degradation.  Table  XI  shows  the  critical  valve  parameters  and  the 
valve-performance  characteristics  that  are  affected  by  changes  in  these  parameters. 
While  most  of  the  characteristics  are  described  in  conventional  terms,  an  explanation  of 
some  of  the  terms  is  helpful. 

The  blocked- spool  pressure  gain  is  a  measure  of  the  ability  of  the  valve  to  over¬ 
come  frictional  forces  on  the  spool  when  the  sp^  1  i  ,  blocked.  I’i  is  defined  as  the  net 
axial  force  on  the  spool,  created  by  the  differential  control  pressure,  per  milliampere 
of  differential  current  input  per  inch  of  spool  circumference,  i,  e. , 

dg  lb 

Blocked-spool  pressure  gain  e  — ; -  —7-  ■"  '  , 

1  4  ma  in.  ' 


(See  list  of  symbols). 

It  is  assumed  that  the  frictional  force  on  the  spoc  l  Is  ci  eated  hy  particles  that  nav*^ 
been  forced  into  the  spool  clearance  space  and  that  the  magnitude  of  this  force  is  propor¬ 
tional  to  the  number  of  particles  entrapped  around  tue  spool.  For  any  clearance,  the 
number  of  entrapped  particles  is  a  ^unction  of  tae  spool  cii  Cv.mference. 

The  threshold  current  is  the  amo'  nt  of  input  current  required  to  initiate  spool 
motion,  and  is  related  to  the  blocked-spool  pressure  gain. 

Null  bias  is  the  amount  of  input  current  required  to  center  the  valve  spool.  A 
change  in  bias,  or  a  shift  in  the  null,  develops  if  the  orific  e  or  nozzle  on  one  side  of  the 
spool  wears  at  a  different  rate  from  that  on  the  other  sidr . 
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After  initial  study  of  the  existing  valve  design,  it  was  concluded  Aat  reductions  in 
sensitivity  to  degradation  and  improvements  in  the  nominal  /alue  of  performance 
characteristics  might  best  be  achieved  by: 

(1)  Changing  the  relative  si^e  of  the  nozzle  and  orifice  diameters 

,2)  Decreasing  the  slot  width  and  spool  diameter,  and  decreasing  the 
spring  rate  of  the  feedback  wire  to  increase  the  spool  stroke. 

The  relative  size  of  the  nozzle  and  orifice  diameters  is  expressed  as  a  ratio  of 
nozzle  diameter  to  orifice  diameter, 

d 

D 

Figures  42  through  45  show  the  effect  of  a  design  change  in  the  ~  ratio  on  the  sig¬ 
nificant  valve  performance  characteristics  and  on  the  normalized  parnal  derivates  of 
quiescent  leakage,  static  flow  gain,  and  control  pressure  gain  with  respect  to  changes  in 
the  nozzle,  orifice,  ard  discharge  coefficients.  Figure  4fi  shows  the  change  in  the  input 
current  required  to  bring  the  spool  to  the  null  position  for  a  given  percentage  change  in 
the  orifice  and  nozzle  coefficients  on  one  side  of  tlie  valve.  The  nozzle  discharge  area 
is  equal  to  TrDh;  to  facilitate  manipulating  this  q>iantity  on  the  analog  computer  model,  the 
parameter  h  was  chosen  as  the  variable  to  represent  discharge  area.  Thus,  the  partlals 
are  expressed  relative  to  the  parameter  h.  The  data  for  these  curves  were  obtained  from 
the  analog  model  of  the  servo  valve  d.iveloped  during  Project  8225,  Task  822502,  and 
described  in  FDL-TDR-64-30. 


Ihe  existing  design  condition  is  represented  by  a  ratio  of  2.  02.  New  values 
of  this  ratio  wc^re  obtained  by  increasing  tbi  area  of  the  orifice  and  decreasing  the 
discharge  area  of  the  nozzle  (the  circumferential  orifice  between  the  nozzle  tip  and  the 
flapper)  by  the  same  factor  and  vice  versa. 

In  Figure  42,  frequency  response  is  represented  as  the  decrease  in  flow  (or  spool 
travel)  between  5  cycles  per  second  and  150  cycles  per  second  expressed  in  dscibsls. 

An  additional  dynamic  characteristic,  maximum  armature  deflection,  is  shown  on  the 
curves.  This  characteristic  is  the  observed  maximum  armature  deueetion  at  130  cycles 
per  second  and  z4-.na  differential-curront  input,  and  represents  the  tendency  of  the 
flapper  to  hit  the  nozzle  tip  under  dynamic  conditions.  The  possible  design  alternatives 
were  determined,  in  part,  on  the  basis  that  the  armature  deflsction  in  the  new  design 
would  be  no  £  eater  than  that  in  the  original  design.  This  assumes,  therefore,  that  the 
ratb  of  degradation  of  the  nozzle  discharge  orifice  in  the  new  design,  due  to  impingement 
against  the  nozzle  tip,  is  no  greater  than  in  the  original  design. 

The  rcflults  of  design  changes  in  the  second  stage  of  the  valve  are  shown  in  Figures 
47  through  52.  All  the  curves  were  obtained  with  the  4  ratio  at  the  existing  value  of 
2.  02.  Figure  47  shows  the  effect  on  the  performance  characteristics  of  a  design  change 
in  the  spool  diameter.  Figures  48  and  49  show  the  effect  of  a  decrease  in  slot  width  and 
feedback-wire  spring  rate  (to  keep  flow  gaip  constant)  on  the  performance  character istics 
at  a  given  spool  diameter.  Figures  50  and  51  show  the  effect  of  a  change  in  slot  width  cn 
the  normalized  partial  derivative  of  quiescent  leakage  with  respect  to  the  valve  uiiderlap 
and  spool  radial  clearance.  Computer  runs  made  at  two  different  spool  diameters  show 
that  these  partlals  are  not  affected  by  changes  in  spool  diameter.  Figure  52  shows  the 
effect  of  changes  in  spool  diameter  and  feedback-wire  spring  rate  on  the  blocked-spoul 
pressure  gain. 
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FIGURE  42.  EFFECT  OF  D/d  RATIO  ON  PERFORMANCE  CHARACTERISTICS  OF  ELECTROHYDRAULIC 
SERVO  VALVE 
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FIGURE  44.  EFFECT  OF  D/d  RATIO  ON  NORMALIZED  PARTIAL  DERIVATIVES  OF  PERFORMANCE 
WITH  RESPECT  TO  NOZZLE  DISCHARGE  COEFFICIENT 
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FIGURE  48.  EFFECT  OF  SLOT  WIDTH  ON  PERFORMANCE  CHARACTERISTICS  FOR  CONSTANT 
FLOW  GAIN  (SPOOL  DIAM,  0.  7  d^;  D/d  RATIO,  2,  02) 
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FIGURE  50.  EFFECT  OF  SLOT  WIDTH  ON  NORMALIZED  QUIESCENT  LEAKAGE  PARTIALS 
WITH  RESPECT  TO  YsLAP  AN*^  CrS  (SPOOL  DIAM,  0.  7  d,) 
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FIGURE  52.  EFFECT  OF  SPOOL  DIAMETLR  AND  FEEDBACK-WIRE  STIFFNESS  ON  BLOCKED-SPOOL 
PRESSURE  GAIN  FOR  CONSTANT  FLOW  GAIN 
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Since  the  valve  is  rated  according  to  its  flow  gain,  all  design  changes  in  the  second 
iitage  were  made  with  the  limitation  that  the  static  flow  gain  of  the  valve  remain  the  same. 
If  the  slot  width  is  reduced,  the  spool  must  travel  further  to  provide  the  same  output 
flow,  and  this  is  obtained  by  a  proportionate  reduction  in  the  stiffness  of  the  feedback 
wire.  Thus,  a  direct  relationship  exists  between  changes  in  slot  width  and  required 
changes  in  feedback-wire  stiffness. 

From  the  curves,  it  can  be  seen  that  quiescent  leakage  and  the  partial  derivative  of 
this  leakage  with  respect  to  underlap  and  spool  clearance  can  be  reduced  if  tlie  slot  width 
(and  correspondingly,  the  feedbt ;k-wire  spring  rate)  is  reduced.  From  the  investigation 
of  the  first  stage,  it  is  seen  that  a  reduction  in  the  ^  ratio  can  produce  a  net  reduction 
in  the  normalized  partial  of  flow  gain  with  respect  to  nozzle  and  orifice  degradation  (the 
decrease  in  the  partial  derivative  with  respect  to  a  change  in  nozale  coefficient  is  greater 
than  the  increase  in  the  derivative  with  respect  to  change  ^  orifice  coefficient).  The 
other  effects  produced  by  a  decrease  in  slot  width  and  in  ratio  must,  however,  be 
correlated  to  arrive  at  possible  design  alternatives. 

Table  XII  summarizes  the  effects  produced  by  the  various  possible  design  changes. 


Determination  of  Characteristics  of  Improved  Design.  A  decrease  in  the  ~  ratio 
and  in  the  slot  width  must  be  accompanied  by  a  decrease  in  the  spool  diameter  in  order  to 
keep  the  maximum  armature  deflection  within  the  specified  limitation  and  the  frequency 
response  as  flat  as  in  the  existing  valve.  A  decrease  in  spool  diameter  Is  accompanied, 
However,  by  a  decrease  in  blocked- spool  pressure  gain,  and  this  must  be  balanced  by 
the  increase  that  results  from  the  decrease  in  die  ^  ratio. 

D 

Figure  53  shows  the  method  by  which  allowable  combinations  of  ^  ratio,  spool 
diameter,  and  slot  width  for  constant  flow  gain  were  determined.  It  is  based  on  Figure 
52  but  with  slot  width  rather  than  feedback-wire  spring  rate  as  the  abscissa.  Loci  are 
drawn  that  locate  the  minimum  allowable  slot  width  at  any  spool  «li?.meter  and  ratio 
to  keep  the  maximum  armature  deflection  below  80  percent  at  150  cycles  per  second  (this 
is  the  value  achieved  in  the  present  design).  These  Icci  are  based  on  the  data  in  Figures 
42,  48,  and  49,  and  on  the  assumption  that  the  minimum  allowable  slot  width  is  that  value 
below  which  the  decrease  m  maximum  armature  deflection  di'  to  the  decrease  in  spool 
diameter  no  longer  compensates  for  the  increase  in  armature  deflection  due  to  the  change 
in  sloe  width  and  ratio.  Thus,  for  any  ratio,  the  aJlo'.vable  slot  width  must  lie  to 
the  right  of  the  corresponding  locus  line. 


The  consideration  oi  allowable  blocked-spool  pressure  gain  establishes  a  second 
limitation  on  the  design  alternatives.  The  decrease  in  the  blocked  spool  pressure  gain 
resulting  from  a  decrease  in  spool  diameter  must  be  balanced,  essentially,  by  the  in¬ 
crease  that  results  from  a  decrease  in  the  -9  ratio.  Tnus,  a  new  set  of  loci  can  be 
shown  to  locate  the  minimum  allowable  blocked-spool  pressure  gain  for  any  combination 
of  sp'^ol  diameter,  slot  width,  and  ^  ratio.  In  this  case,  the  allowable  combinations  of 
slot  width  and  spool  diameter  for  any  ratio  must  lie  above  the  corresponding  locus  line. 
These  two  limitations  now  ettiblish  the  possible  new  design  combinations  since,  for  any 
ratio,  the  combination  of  spool  diameter,  slot  width,  and  the  corresponding  feedback- 
wire  spring  rate  must  lie  to  the  right  and  above  the  intersection  of  the  two  appropriate 
loci.  The  set  of  intersections  establishes  an  envelope  of  poaeible  designs,  and  the  most 
desirable  point  cn  this  envelope  is  that  which  results  in  a  minimum  slot  width.  This 
point  on  the  envelope  is  represented  by  the  following  parameters,  which  describe  a  new 
valve  design; 
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TABLE  Xn.  SUMMARY  OF  EFFECTS  OF  POSSIBLE  DESIGN  CHANCES 
ON  VALVE  CHARACTERISTICS 


Design  Change^*^ 

V«lv«  Characteristic 

Decrease 

in 

D 

d 

Decrease 

in 

Slot  Width 

Decrease 

in 

Spool 

Diameter 

Static  flow  gain 

+ 

Constant 

Constant 

Quiescent  leakage 

+ 

+ 

0 

Frequency  response  (db) 

- 

- 

+ 

Maximum  armature  deflection 

- 

- 

+ 

Blocked-spool  pressure  gafn 
Quiescent  leakage  partial 

+ 

+ 

* 

with  respect  to  Cq  and  C^^ 
Quiescent  leakage  partial 
with  respect  to  YeL,Ap 

+ 

Not  applicable- 

Not  applicable 

^RS 

Not  applicable 

+ 

0 

Null  shift  current 

Flow  gain  partial  with 

0 

Not  applicable 

No*'  applicable 

respect  to  Cq  and  Cj^ 

+ 

Not  applicable 

Not  applicable 

(«)  ♦  lndic«K«  desirable  change. 

-  Indicawi  undetiiable  change. 

0  Indlcacei  Insignincant  change. 
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a  1.57 


_D 
d 

D  s  0.  0135  in. 

d  T  0.  00862  in. 

s  0.  ii250  in. 
u  =  0.  150  in. 

Kfba  =  11*5  Ib/in. 

Table  XIll  gives  a  summary  comparison  of  the  performs'. ce  characteristics  and 
normalized  partial  derivatives  of  the  new  design  with  the  existing  design.  The  Per  Cent 
Change  column  shows  the  improvement  in  the  sum  of  the  normalized  partial  derivatives, 
assuming  an  equal  per  cent  change  in  each  parameter. 

For  comparable  flow  gains  and  blocked-spool  pressure  gain,  the  nominal  quiescent 
leakage  and  valve  attenuation  have  been  decreased,  and  the  sensitivity  of  the  flow  gain, 
blocked- spool  pressure  gain,  quiescent  leakage,  and  frequency  response  to  parameter 
degradation  have  been  reduced. 

Measurements  of  differential-load  pressure  gain  were  not  made  for  all  design 
conditions,  because  the  high  ratio  of  this  gain  makes  it  difficult  tc  get  an  accurate  mea¬ 
surement  on  the  analog  computer.  The  strut  pressure  gain  is,  however,  inversely 
proportional  to  quiescent  leakage  at  small  values  of  input  current,  and  it  can  be  assumed 
that  the  changes  in  pressure  gain  are  related  directly  to  changes  in  quiescent  leakage. 

The  reduction  in  the  ratio  to  1.57  was  accomplished  by  increasing  the  inlet- 
orifice  diameter  and  decreasing  the  nozzle  diameter.  This  has  the  additional  benefit  of 
reducing  the  tendency  of  the  orifice  to  be  blocked  by  contaminant  particles.  Since  the 
orifice  is  the  smaller  of  the  two  diameters  and  has,  therefore,  a  greater  susceptibility 
to  being  blocked,  an  increase  in  this  diameter  produces  an  increase  in  reliability. 

It  is  clear  from  the  results  of  this  design  investigation  that,  while  some  improve¬ 
ment  in  valve  reliability  can  be  achieved,  the  parameters  selected  by  the  manufacturer 
in  the  present  design  configuration  are  very  close  to  their  optinium  values. 


Predicting  Relative  Reliability 

The  next  subtask  of  project  activity  was  to  estimate  quantitatively  the  improvement 
in  valve  reliability  that  could  be  anticipated  from  the  new  valve  design.  The  analysis 
was  based  on  the  performance  characteristic  of  quiescent  leakage  at  null.  The  relative 
reliability  was  expressed  as  the  ratio  of  the  expected  number  of  cycles  of  operation  of 
the  new  valve  design  for  the  quiescent  leakage  to  reach  a  specified  limit  to  the  number  of 
cycles  of  the  original  valve  design  to  reach  the  same  limit.  Quiescent  leakage  was  chosen 
as  the  focus  of  study  because  past  laboratory  experience  had  shown  that  this  valve  char¬ 
acteristic  was  most  susceptible  to  degradation,  and  that  the  change  in  performance  over 
the  period  of  the  life  test  could  be  readily  obsc»‘ved. 
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TABLE  Xin.  COMPARISON  OF  ORIGINAL  AND  NEW  DESIGN 


Value 

Original 

E/esig-v 

New 

Design 

Change,  % 

Qy  at  4ma 

18. 64  cis 

18. 59  cis 

6  at  47na 

5.06% 

7.52% 

Y3  at  4ma 

28. 74% 

29.82% 

APj^  at  4ma 

125. 2  psi 

16'j.  8  psi 

Pj^j  at  4ma 

1844  psi 

2240  psi 

at  null 

1. 240  cis 

1 .  13 j  cis 

-8.6 

APj^  (ds)/(4i) 

17.2  Ib/ma  in. 

16.  ?  Tb/ma  in. 

FR  at  150  cps 

-5. 98  db 

-4.38  db 

-26.7 

<1>  at  IS?  cps 

-64. 2  deg 

-54.  5  deg 

6  at  150  cps 

80% 

80% 

^^Null^^^l 

0,  024  ma/% 

-0.0241  ma/%  3 

0 

^^Null/^^Ol 

-0.021  ma/% 

-0.  0195  ma/%  / 

aQv/3c 

-0, 001  cis/% 

-0.0083  cis/%  1 

-28 

-0.014  ci3/% 

-0.0025  cis/%  / 

a  {APj^dg/4i) 

aco 

-0.  013  Ib/ma  in. 

-0.  058  Ib/ma  in.”j 

»  -32 

a{APNd8/4i) 

3Cn 

-0.  093  Ib/ma  in. 

-0.  016  Ib/ma  in.  > 

3Q£/Sc 

+0. 0015  cis/% 

+0.0012  ci8/%  1 

+&. 0016  cis/% 

+0.  002  ci8/% 

-7 

ao^/aYsLAP 

+C. 0013  ci8/% 

0.00125  cis/% 

ao^/acRs 

+0.  0227.  ci8/% 

0.0203  cis/%  J 

aFR/aCo  at  150 

-0.033  db/% 

-0.  022  db/%  3 

-38 

aFR/aCj^  at  150 

-0. 056  db/% 

-0.033  db/%  i 

i  < 


4 


I 


r>C»tr  -  -  . . . ,  _  .  , 

. . . . .  —  . . . — I  . . . umi — 7^ - f..,— ...>,>1.1  .  i  »  y. 


The  prediction  of  reliability  involves: 

(1)  Calculating  the  variance  that  cun  be  expected  in  the  performance 
characteristics  under  study  as  a  function  of  the  manufacturing 
tolerances  of  the  critical  parameters 

(2)  Deriving  a  mathematical  expression  for  the  change  in  the  mean 
periarmance  characteristics  resulting  from  degradation  in  the  critical 
valve  parameters  during  operation 

Combining  the  variance  and  average  degradation  equation  to  determine 
the  range  cf  cycles  of  operation  during  which  the  performance 
specifications  \?ill  be  exceeded. 

The  parameters  that  have  the  most  cignificar.t  influence  cn  the  quiescent  leakage 
and  are  subject  to  degradation  during  operation  are; 

(1)  Noszie  flow  coefficient  (Cj^) 

(2)  Orifice  flow  coefficient  (Cq) 

(3)  Clearance  between  spool  a;  d  bushing  (Cnc-) 

XV 

(4)  Magnitude  of  spool  underlap 

Tne  quiescent  leakage  at  null  is  expressed  by  the  equation; 

Qg  s  Qi  i.  ♦  ^^^SLLl-DP  ‘^BLLl-DP^  '  **DP^  ^  ^^^PSL^^S2^  ^^Ll  " 

-p 

\Kp/  V  'Ll  DP/  SLAP  ’ 

where 

Qj  *TrDCN(2/p)‘'^h(PNi  - 
Qg  =•  nD  C^{Z/o)^^^  h  (Pj^2  ‘ 

^SLLl-DP  *  Cps^/12fxg 

^BLLl-DP  *  ’'^^B 
CpsL  * 

Ysg.  1.47 

1/K  p  Xsma:^^ 

P  p  X  104 

In  the  original  design,  Q g  ■  1 . 240  cis  and  in  the  new  design  it  is  reduced  to 
1. 198  cis. 
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Calculation  of  Leakage  Variance.  The  variance  that  can  be  expected  tor  the  above 
leakage  values  is  calculated  from  the  partial,  derivatives  developed  from  the  analog 
model  and  from  the  variance  in  the  valve  pararr  eters  by  means  oi  the  Propagation  of 
Variancf?  formula  as  given  in  Figi-.re  1. 

It  is  assumed  that  the  noazle  and  orifice  flew  coefficients  have  zero  variance  and 
that  the  leakage  from  the  end  of  the  spool  through  the  spool  clearance  is  small,  so  that 
the  variance  in  the  nozzle  flow  at  null  is  due  only  to  the  manufac' uring  tolerances  in  the 
nozzle  and  orifice  diameters.  The  laminar  spool  end  leakage,  calculated 

from  the  equation: 


®END  * 


ZTTds 

IZy-i 


^RS 


Values  for  QeND  were  calculated  to  be  0.  00372  cio  and  0.  00223  cis  foi  the  original  and 
new  designs,  respectively,  *  1783  psi  in  the  original  design  and  2160  psi  in  the  new 

design,  as  obtained  from  the  analog  model. 


The  partial  derivitaves  with  respect  to  flow  area  can  be  obtained  from  the  partials 
with  respect  to  flow  coefficient  by  noting  that  the  normalized  values  of  these  partials  are 
equal.  This  results  from  the  fact  that  the  flow  coefficient  is  simply  an  area  correction 
factor. 


We  can  thus  write  that: 


and  that: 


iSi. 


, 

and  -rp; — 

dCjsj 


the  analog  model. 


do 

dCp  dCj^  dAf^ 

Cq  Aq  Cm  An 


dQjt 

-  .^9 

dOi  ^ 

aof 

cn 

dQl 

c)Aq 

aco 

"  An 

are  the  absolute  values  of  the  partial  derivatives  as  obtained  from 


It  can  further  be  assumed  that  the  dimensional  tolerances  in  the  nozzle  and  orifice 
diameters  represent  the  range  of  ±3  standard  deviations  or  99.9  percent  of  all  parts 
manufactured.  The  nozzle-diameter  tolerances  of  ±0.0005  in.  and  the  orifice  tolerance 
of  ±0.0001  in.  then  give  the  standard  deviations  in  the  respective  flow  areas,  as  listed 
in  Table  XIV. 

TABLE  XIV.  STANDARD  DEVIATIONS  IN  NOZZLE 
AND  ORIFICE  AREAS 


Original 

Design 

New  Design 

Nozzle  area 

0. 042  X 

lO-**  in.  ^ 

0.  036  X  10"4  in.  2 

Orifice  area 

0. 0042  X 

c 

1 

o 

0.0045  X  10“^  in.  ^ 

.03 


The  spool  clearance  and  undcrlap  arc  the  same  in  both  designs,  and  the  standard 
deviations  in  these  values,  due  to  manufacturing  tolerances,  are  as  given  in  FOL<-TDR- 
64>50,  page  52: 

-6 

Spool  Clearance  Deviation  ■  1.67  x  10  in. 

Spool  Underlap  Deviation  a  0.29  x  10~^  percent. 

These  variances  taken  with  the  computed  partial  derivatives  then  result  in  quiescent 
leakage  standard  deviations  of  0.0676  cis  for  the  original  design  and  0.0618  cis  for  the 
new  design. 

If  a  tolerance  range  of  ±3  standard  deviations  is  assumed,  the  expected  quiescent 
leakage  from  the  servo  valves  is: 

Original  Design:  1.24  ±  0.20  cis  or  0.322  ±  0,053  gpm. 

New  Design;  1.  198  ±  0.  19  cis  or  0.309  ±  0.048  gpm. 

The  quiescent  leakage  will  be  affected  by  two  types  of  wear  degradation.  The  first 
is  due  to  erosion  of  metal  surfaces  by  the  impact  of  solid  particles  contained  in  the  fluid, 
and  the  second  is  abrasive  wear  due  to  the  sliding  action  between  contacting  metal 
surfaces.  The  leakage  that  is  affected  by  erosive  degradation  will  be  identified  as 
and  is  equal  to  the  noazle  flow  (Q  j  +  Q^)  less  the  laminar  flow  out  the  ends  of  the  spool, 
plus  the  flow  through  the  lap  space.  The  abrasive  wear  will  affect  the  radial  spool 
clearance,  and  the  leakage  through  this  clearance  will  be  identified  as  The  total 

quiescent  leakage  is  then  Q/  ■  Q^a  *  Qfe  +  bushing  leakage. 

The  sUndard  deviations  for  the  degradable  leakage  components  is  as  follows: 

Original  Design  New  Design 

0.  0054  cis  0.  00626  cis 

0.  0675  cis  0.  0616  cis 


Qfe 

Cifa 


Leakage  Degradation  Equation.  The  expression  for  the  change  in  leakage  as  a 
function  of  time  or  of  cycles  of  operation  is  based  on  the  differential  equations  of  the 


form: 


dAj 

dN 


m 

V  ^  £PJ 

Li  dPj  dV 

J*1 


dV  , 

-  for  erosive  wear 

dN 


(32a) 


and 


dAi  \  dAt  *^^3  ,  . 

—  -  at  >  ■  ■  V  for  abrasive  wear, 

dN  L  dPr  dN 

J=1  ^ 


(32b) 


where 


dAj 

dN 


is  the  change  in  performance  characteristic,  A^,  per  cycle  of  operation. 


— ^  is  the  partial  derivative  of  the  performance  characteristic  with  respect 
dPj  to  the  parameter  Pj,  as  obtained  from  the  mathematical  model. 

dPj 

is  the  change  in  the  value  of  the  parameter  Pj  per  unit  volume  of 
material  removed  and  is  obtained  from  the  geometrical  configuration 
of  the  valve  component  in  question. 


is  the  total  volume  removed  V,  per  cycle.  For  erosive  wear,  this 
differential  is  the  product  of  the  volume  of  material  removed  per 
particle  (v),  number  of  particles  per  unit  volume  of  fluid  and 

the  volume  of  fluid  passing  through  the  region  psr  cycle- 

The  volume  of  material  removed  per  particle  v  can  be  expressed  as 

V  ■  u,  where  Kj  is  a  constant  that  depends  on: 

(1)  The  mass  of  the  particle,  or  the  average  particle  mass 

{Z)  The  angle  of  incidence  of  the  particle  on  the  surface  being  eroded 

(3)  The  plastic  flow  stress  of  the  material, 

and  u  is  the  average  fluid  velocity  in  the  region  of  erosion. 

The  number  of  particles  per  unit  volt.me  of  fluid  is  determined  from  a  fluid 
contamination  count. 

The  fluid  flow  per  cycle  is  a  function  of  the  flow  at  null  (q)  and  is  expressed  as 
K3q,  where  Kj  is  a  constant. 


Then, 


a  K^u  X  X  K^q. 


Since  V  «  0  when  N  ■  0  ,  Vs  KuqN, 


where  K  a  K^K^K^. 


.  is  the  rate  of  change  of  parameter,  Pj,  per  cycle,  N.  For  abrasive 

wear  in  the  servo  valve,  this  parameter  is  the  radial  clearance. 

* 

Degradation  is  assumed  to  be  a  function  of  the  distance  traveled 
by  the  spool  per  cycle  and  is  expressed  as: 

■  ■  My,  where  M  is  a  constant  that  depends  on: 

(1)  The  material  characteristics 

(2)  The  magnitude  of  the  input  signal  in  ma 

(3)  An  arbitrary  maximum  spool  displacement,  and 


y  is  a  measure  of  the  amplitude  of  spool  motion  and  is  expressed  as  a  percentage  of  the 
arbitrary  maximum  spool  amplitude  per  milliampere  of  input  current.  Equations  (32a) 
and  (32b)  can  now  be  rewritten  for  quiescent  leakage  as: 


E/oiion  of  SurUcet  by  Solid  Particiei".  Wei.  March.  I960. 


i 


dOie 

dN 


!2i 


dCp 


dV 


O 


dCj^  dV 


*^n“n*5n  * 


dQ/ 


N 


dYsLAP 


^^SLiAP 


K y^y^y 


(33) 


where  the  cubscripte  O,  N,  and  Y  refer  to  the  parameter  values  of  the  orifice  nozr,le  and 
lap  space,  respectively. 


(34) 


dCp 

The  values  of  “jr?*, 


dV, 


iSiL 


dV 


and 


dY 


SLAP 


N 


dVy 


are  derived  from  the  geometry  of  the 


valve  elements.  If  one  now  considers  the  factors  affecting  the  magnitudes  of  Kq, 
and  Ky*  above,  it  could  be  assumed  that  these  constants  are  to  some  first  order  of 
approximation  almost  equal,  or  that  they  at  least  a?‘e  of  the  same  order  of  magnitude. 
If  it  is  initially  assumed  that  tliey  are  equal,  it  is  possible  to  further  simplify  Equation 
(33):  in  doing  so,  it  can  be  shown  that  the  assumption  that  Kq, 
same  order  of  magnitude  is  sufficient  for  the  simplification. 


and  Ky  be  of  the 


If  we  let  Kq  ••  Kjyj  »»  Ky  ••  K  and  divide  by  K,  we  can  write  Equation  (33)  as: 


Then,  since  uq  ■ 


uQqo  + 


hQt  dC 


N 


ac 


N 


dV 


Wn  + 


dY 


SLAP 


N 


aY 


SLAP 


dV, 


Uyqy 


KN 


we  can  rewrite  this  as: 


dCo 

'^^SLAP 

h 

dKN 

. 

dVo 

KN 

dVj,, 

^^SLAP 

dV 

Y 

KN 

dQf, 


(35) 


versus  KN 


By  assuming  that  the  partial  derivatives  remain  constant,  a  curve  of 

can  be  plotted.  The  area  under  this  curve  between  any  two  values  of  KN  represents  the 
change  in  between  these  values.  It  can  be  shown  that  the  terms  relating  to  wear  in 
the  orifice  and  nozzle  are  insignificant  when  compared  with  the  terms  involving  YslAP’ 
Thus,  as  long  as  Kq,  Kj^,  and  Ky  are  of  the  same  order  of  magnitude.  Equation  (33)  can 
be  simplified  by  considering  only  the  effect  of  changes  in  spool  underlap.  The  same 
conclusion  can  be  drawn  for  the  new  design. 


Equation  (35)  then  becomes: 


dC£e  ‘^^SLAP 

^Kn  ^Ygi^p  dVy  KN 

From  the  supplemental  material; 

tlYsLAP  1  6840 

— — -  ■  6840  L-r^  ■  - 

'^^Y  ^  ^Y  yuyqyKN 

substituting  this  into  Equation  (36)  and  integrating  gives; 


(36) 
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(37) 


Q£g  8  6840  - - yuyqyKN  +  ,  where! 

^’’^SLAP 

Q^'  is  the  initial  value  of  Q^^- 

Now  by  integrating  Equation  (34),  the  abrasive  degradation  equation  is  obtained: 

bOi 

Qia  -  y  M  N  +  Q^;  .  (38) 

Table  15  lists  the  calculated  values  of  the  parameters  of  Equations  (37)  and  (38)  as 
applied  to  the  original  and  the  new  valve  designs. 


TABLE  XV.  VALUES  OF  QUIESCENT  LEAKAGE 
DEGRADATION  PARAMETERS 


Parameter 

Original  Design 

New  Design 

dQj^/^YgLAp 

1. 33  cis/%  max  stroke 

1.  25  cie/%  max  stroke 

dOf/aCg^ 

0.  0404  X  10^  cis/in. 

0.0369  X  10^  cis/in. 

Uy 

4520  in. /sec 

4520  in. /sec 

qy 

0.  0162  cis 

0.  0156  cis 

Q£e 

0.  5341 

0. 5222 

Q£a 

0.6839  .  . 

0.  6441 

y 

11 

B 

00 

7.46%/ma(^) 

(a)  Baled  on  dau  from  analog  model. 


From  the  values  in  Table  15  and  the  consideration  of  a  ±3  standard  deviation 
tolerance,  the  final  leakage  degradation  equations  are: 

Original  Design 


8  0.  1555  X  10^  Vlw 
Q£a  »  0.  29  X  106  M  N 


+  0,  5341  dk  0.  0162  cis 
+  0.  684  i  0.  2025  cis 


New  Design 

Q£g  a  0.  1488  X  10^  Vk  N  +  0.  5222  ±  0.  0188  cis* 
Ufa  *  0-  275  X  10^  M  N  +  0.  644  ±  0.  1848  cis 


(39) 

(40) 


(41) 

(42) 


The  bushing  leakage  is  calculated  from  the 
has  a  value  of  0.  022  cis  for  both  designs. 


equation  j_j2p(PLi  -  Pop)  and 


D*r  over  t  lymbol  indicates  ne  new  design. 
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Relative  ReliabiUty.  If  the  limiting  value  of  quiescent  leakage  is  taken  &s  2  cis, 
and  since  +  bushing  leakage,  then,  using  Equations  (39)  through  (42), 

expressions  for  numbers  of  cycles  to  failure  for  the  two  designs  are: 

0.  1555  X  10^  \/irTrj  +  0.29  X  10^  M  Nj  ■  0.760*  0.219/for  the  original  design 

0.  1488  X  10^  Vk  Nj  +  0,27S  x  10^  M  Nj  ■  0.812*  0.204/for  the  new  design, 

/ 

where  N|  and  Nj  represent  the  number  of  cycles  to  failure  for  the  originp.l  a.id  new  de¬ 
signs  respectively. 

_  Depending  on  the  values  of  K  and  M,  we  would  expect,  on  the  average,  that  the  ratio 
Ni/Ni,  the  relative  reliability  for  the  two  designs,  would  be  between  1.  125  and  1.25. 

This  indicates  a  net  improvement  in  valve  reliability.  If  K  is  much  larger  than  M  so  that 
the  abrasive  degradation  can  be  neglected,  the  higher  ratio  could  be  expected. 

It  should  be  noted,  however,  that  the  above  ratios  represent  only  mean  values  and 
that,  in  a  direct  comparison  of  a  small  sample  of  valves,  the  relative  reliabilities  will 
not  necessarily  fall  on  this  average.  In  fact,  if  t.  sample  of,  say,  two  valves  of  each 
design  is  taken,  one  should  consider  that  the  valve  performance  could  faU  anywhere 
within  a  range  of  *3  standard  deviations  from  the  average.  Thus,  at  the  extreme,  if  the 
leakage  of  the  new  design  were  at  the  high  limit  and  that  of  the  original  design  at  the  low 
limit,  cne  could  expect  relative  reliability  ratio  to  be  about  0.  65,  which  would  indicate 
poorer  reliability.  If,  on  the  other  hand,  the  new  design  were  at  the  low  limit  and  the 
original  design  at  the  high  limit,  a  96  percent  improvement  in  reliability  could  be 
expected. 


Experimental  Verification  of  Reliability  Prediction 


Construction  of  Modified  S?rvo  Valves.  In  order  to  examine  experimentally  the 
reliability  resulting  from  an  application  of  the  Moment  method  analysis  to  a  two-stage 
servo  valve,  a  sample  of  four  valves  was  ordered  from  a  supplier.  Two  of  the  valves 
incorporated  the  design  changes  derived  'rom  the  Moment  method  analysis.  The  other 
twc  were  of  conventional  design  for  comparison. 

Although  the  two  modified  valves  formed  a  small  sample,  it  was  hoped  that  they 
would  provide  useful  information  concerning  the  effectiveness  of  the  application  of  che 
Moment  method  to  the  reliability  optimization  of  an  actual  device.  As  stated  earlier,  the 
servo  valve  design  was  found  to  be  already  close  to  "optimum.  However,  some  signifi¬ 
cant  changes  were  evolved  that  would  improve  u»e  predicted  reliability.  The  sample 
valves  were  gradually  degraded  by  life  testing  under  continuous  cycling.  A  record  of  the 
progress  of  degradation  was  kept  by  periodically  removing  the  valves  frt  m  the  life¬ 
testing  stand  and  placing  them  on  an  instrumented  test  stand  to  undergo  a  s-ries  of  tests 
covering  the  major  parameters  contributing  to  the  predicted  reliability  improvements. 


Description  of  Testing  Program.  Initial  tests  (before  lift  testing)  were  made  in 
order  to  establish  a  standard  against  which  degradation  would  be  measured.  The  initial 
tests  applied  to  each  valve  aie  listed  below: 
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(1)  Record  Null 

(2)  Quiescent  Flov/  /ersus  Differential  Current 

(3)  Differential  Load  Pressure  Versus  Differential  Current 

(4)  Short-Circuited  Load  Pressure  Versus  Differential  Current 

(5)  Total  Flow  Versus  Differential  Current 

(6)  Load  Flow  Versus  Differential  Current 

(7)  Frequency  Response  With  Blocked  Load. 

In  later  tests,  this  series  was  shortened  to  cover  the  characteristics  most  directly 
applicable  to  the  results  of  the  Moment  method  analysis: 

(1)  Record  Null 

(^)  Quiescent  Flow  Versus  Differential  Current 

(3)  Differential  Load  Pressure  Versus  Differential  Current 

(4)  Total  Flow  Versus  Differential  Current 

(5)  Frequency  Response  With  Blocked  Load. 

After  initial  testing,  the  valves  were  placed  on  a  life-testing  stand  and  operated  by 
continuously  reversing  the  direction  of  flow.  At  intervals  during  the  life  testing,  the 
valves  were  removed  and  the  individual  tests  were  repeated  to  record  the  progress  of 
degradation  in  the  operating  characteristics  of  the  valves.  The  test  series  was  made  a 
tc'^al  of  seven  times  at  0;  140,  000;  300,  OOC;  900,  000;  1,  600,  000;  3,  000,  OOO;  and  8,  600,  000 
cycles. 

The  final  step  in  the  testing  program  was  the  plotting  of  trend  curves  to  show  the 
rate  of  degradation  of  the  operating  characteristics  as  a  function  of  totel  cycles  or  total 
hours.  The  characteristics  represented  by  the  trend  curves  were  null  bias,  peak 
quiescent  flow,  minimum  quiescent  flow,  difference  between  maximum  and  minimum 
quiescent  flow,  pressure  gain,  and  frequency  response  at  100  cycles  per  second. 


Description  of  Test  Equipment.  The  equipment  used  in  testing  the  valves  consisted 
primarily  of  two  independent  systems:  an  instrumented  hydraulic  test  bench  and  a  lif  - 
testing  stand. 


Instrumented  Test  Bench.  The  test  bench  was  originally  designed  for  the  moni¬ 
toring  of  servo- valve  characteristics.  It  includes  a  pump,  reservoir,  va.iable  load,  and 
instrumentation.  A  schematic  diagram  of  the  test  bench  is  shown  in  Figure  54.  Oil  is 
delivered  under  pre-^sure  to  the  manifold  block  through  a  passage  leading  to  the  valve 
supply  ports.  It  is  then  routed  ly  the  valve  (depending  on  spool  position)  to  either  Load 
Line  A  or  Load  Line  B.  ihe  v.ii  >  ..le-load  device  is  not  shown  on  the  diagram.  It 
consists  of  a  main  valve  and  a  bypasr  valve  for  fine  adjustments.  The  load  may  be  ad¬ 
justed  anywhere  from  fully  blocked  to  fully  open,  i.  e,  very  low  hydraulic  resistance.  A 
closed  loai  valve  is  referrt^d  to  as  '  blocked  load",  and  a  fully  open  load  valve  is  called 
"short-circuited  load".  Except  in  the  null  or  centered  spool  position,  one  load  line  is 
always  connected  .  o  supply  and  the  other  to  drain.  The  main  flow  leakage  and  nozzle  flow 
are  mixed  inside  the  manifold  at  the  manifold  drain  port. 

The  test  bench  instrumentation  consists  of  a  servo  crcplifier,  flowmeters,  pressure 
gages,  and  thermometers.  A  signal  generator  is  used  as  an  accessory  during  frequency 
response  tests.  The  flowmeters  re  ar*'an^ed  in  two  parallel  banks,  each  containing  two 
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Hydraulic  Filters 


FIGURE  54.  SJ'.RVO- VALVE  KYDRAULIC-TEST-STAND  SCHEMATIC 
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meters  in  series.  One  bank  cover  flows  up  to  1  gpm,  and  the  other  registers  flows  from 
1  to  20  gpm.  The  flowmeters  are  of  the  variable- orifice  type,  using  a  metal  float  in  a 
tapered  glass  tube.  The  meters  may  he,  read  easily  to  three  significant  figures.  Pres¬ 
sure  gages  are  provided  at  the  supply  port,  in  both  load  linen,  and  in  the  drain  line.  The 
difference  between  the  load-line  pressure  readings  is  governed  by  the  setting  of  the  load 
valve  and  spool  position.  Oil  is  supplied  to  the  test  bench  circuit  by  a  fixed  displacement, 
pump  delivering  10  gpm  at  pressures  up  to  5000  psi.  Tests  are  normally  run  with  a 
supply  pressure  of  3000  psi,  which  is  held  constant  under  varying  conditions  by  adjusting 
a  bypass  valve.  Flow  from  the  pump  is  smr othed  by  an  accumulator  in  the  supply  line. 

Objectives  and  Procedures  for  Tests.  The  following  describes  the  objectives  and 
summarizes  the  method  used  for  the  five  tests  that  formed  the  series  made  at  each  of  the 
seven  testing  points  during  the  life  tests. 

Record  Null.  The  purpose  of  this  test  is  to  determine  the  amount  of  differential 
current  required  to  place  the  spool  so  that  the  pressures  in  Load  Lines  A  and  B  are 
equal  with  the  load  valve  completely  shut  off.  The  amount  of  current  required  to  achieve 
this  "hydraulic  null"  should  be  very  small,  and  is  a  measure  of  the  difference  between 
electric  null  (zero  differential  current)  and  an  actual  pressure  balance  in  the  valve. 
Ideally,  these  conditions  would  coincide,  though  in  a  real  valve  it  is  only  reasonable  to 
'xpect  a  small  offset  between  them.  A  plot  of  the  "null  bias"  or  current  required  to 
balance  the  valve,  versus  time,  can  be  valuable  in  revealing  changes  in  the  behavior  of 
the  valve  in  the  critical  area  around  null.  Scoring  of  the  spool,  bending  of  the  feedback 
wire,  wear  of  the  feedback-wire  tip,  or  damage  to  the  flapper  can  result  in  changes  in 
the  value  of  the  null  bias  current  with  tinne. 

The  record  null,  or  null  bias  test,  is  run  by  bringing  the  spool  slowly  from  full 
displacement  to  a  pressure-balanced  condition.  This  is  done  from  both  directions  to 
cancel  out  the  effect  of  spool  friction.  In  the  tests,  the  average  of  the  null  diiferential 
curr<-nts  from  both  directions  was  taken  as  the  null  bias,  expressed  in  milliamperes. 

Quiescent  Flow  Versus  Differential  Current.  This  test  is  performed  with  the  load 
valve  shut  off.  The  "quiescent  flow"  measured  is  the  sum  of  the  flow  through  the  nozzles 
and  the  spool  leakage  flow.  The  quies'ent  flow  characteristienliy  reaches  a  peak  near 
null  and  decreases  to  a  lelatively  constant  value  in  the  0.  1  to  0.  2-gpm  range  on  either 
side  of  null.  Since  the  spool  leak  age  is  very  small  away  from  null,  and  since  the  total 
nozzle  flow  remains  essentially  constant  throughout  the  normal  spool  stroke,  the 
quiescent  leakage  curve  may  be  viewed  as  a  portrait  of  variable  spool  leakage  super¬ 
imposed  on  constant  nozzle  flow,  as  a  function  of  spool  position  (spool  position  is  repre¬ 
sented  by  differential  current  in  ther  linear  region  of  valve  operation).  The  maximum 
value  of  quiescent  leakage  is  a  sensitive  indicator  of  spool  metering  edge  wear,  and  was 
a  characterintic  of  primary  importance  in  the  modification  of  the  valve  for  better 
reliability . 

The  quiescent  flow  test  is  run  by  stepping  the  differential  current  away  from  null  in 
both  directions  ("positive"  and  "negative"  differential  current),  and  recording  the  total 
flow  from  the  valvf  at  each  step.  No  load  flow  is  permitted.  A  tyoical  quiescent  flow  test 
curve  is  shown  in  I'igure  55, 
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Differential  L.oad  Pressure  Versus  Differential  Current.  The  differential  load 
pressure  test  is  a  blocked-load  procedure  used  to  evaluate  the  "resolution"  or  propor¬ 
tionality  of  response  of  the  valve  near  hydraulic  null.  The  test  is  performed  by  closing 
tile  load  valve  and  noting  the  difference  between  load-line  pressure  A  and  B  as  the  spool 
is  shifted  slightly  to  either  side  from  the  pressure-balanced  condition.  The  slope  of  the 
differential-pressure  curve  if  known  as  "pressure  gain",  and  is  an  index  of  nozzle,  spool, 
and  meter ing-port  wear.  The  existence  of  a  finite  pressure  gain  is  due  primarily  to  the 
slight  underlap  built  into  the  valve  to  improve  its  proportionality.  Figure  56  shows  a 
typical  differential-load-pressure  curve. 


Total  Flow  Versus  Differential  Current.  In  this  test  the  load  valve  is  opened  and 
adjusted  so  that,  for  any  value  of  differential  current,  a  specified  constant  pressure  drop 
occurs  from  Load  Line  A  to  Load  Line  B,  or  vice  versa.  Supply  pressure  is  also  ad¬ 
justed  (as  in  all  tests)  to  maintain  3000  psi.  Differential  current,  i.e.  ,  spool  displace¬ 
ment,  is  increased  from  null  in  steps  under  the  above  conditions,  and  total  flow  is 
recorded  at  each  step.  The  resulting  graph  takes  the  characteristic  form  of  two  lines 
sloping  upward  in  a  vee  from  the  origin,  which  represents  zero  differential  current.  The 
right  leg  of  the  vee  shows  an  increase  in  flow  with  increasing  spool  motion  in  one  direc¬ 
tion  through  the  simulated  load  reuiatance,,  and  the  left  leg  represents  a  similar  flow 
increase  in  the  opposite  direction.  At  the  vertex  of  the  vee,  only  quiescent  flow  occurs. 
One  of  the  objectives  of  the  cotal-flow-versus-differential  current  test  is  to  illustrate 
graphically  the  extent  to  which  flow  increases  linearly  with  increasing  differential 
current.  The  slope  of  the  total-flow  curve  is  called  the  "flow  gain".  The  flow  gain  gen¬ 
erally  remains  quite  constant  for  differential  currents  from  zero  up  to  about  6  ma.  It 
should  be  emphasized  that,  in  order  to  obtain  this  linearity,  both  supply  pressure  and 
load-pressure  drop  must  be  held  at  fixed  values  by  varying  the  supply  bypass  and  load- 
valve  settings  after  each  change  in  differential  current.  This  is,  of  course,  a  static  or 
steady- state  test,  as  are  all  those  in  this  series  except  frequency  response.  Figure  57 
shows  a  represent.- tive  total-flow  curve. 


Frequency-Response  Test.  For  this  program,  a  form  of  frequency-response  test 
was  devised  to  conform  to  the  available  instrumentation  on  the  static-test  stand.  The 
test  was  run  with  the  load  valve  closed,  and  the  spool  was  oscillated  about  hydraulic  null 
with  a  sinusoidal  torque-motor  signal  of  7  millivolts  peak-to-peak  magnitude,  and  fre¬ 
quencies  ranging  from  5  to  150  cps.  A  signal  generator  was  connected  to  the  regular 
test-stand  d-c  amplifier  to  provide  the  sinusoidal  input.  The  output,  or  dependent  varia¬ 
ble,  was  the  mean  flow  as  read  visually  from  the  flowmeter.  The  response  of  the  flow¬ 
meter  was  smooth  above  5  cps,  and  provided  an  accurate  and  repeatable  reading.  The 
individual  test  curves  were  plotted  directly  as  recorded,  giving  flow  in  gpm  versus  fre¬ 
quency  in  cps.  For  the  trend  curves  representing  changes  from  test  to  test,  the  ratio  of 
flow  at  100  cps  to  maximum  flow  was  recorded,  in  decibels,  as  a  function  of  total  cycles 
on  the  valves.  A  representative  frequency-response  curve  is  shown  in  Figure  58. 


Application  of  the  Moment  Method.  Before  reviewing  the  results  of  the  life  tests  of 
the  two  modified  and  two  unmodified  nervo  valves,  it  will  be  helpful  to  review  some  of 
the  general  features  of  the  Moment  method  of  statistical  inference. 

The  Moment  method  cannot,  by  itself,  take  time  into  account  as  a  variable.  It 
yields  only  the  first  and  second  moments  (mean  and  variance)  of  a 
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P/^-Pq=  2000  psi 
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FIGURE  57.  VALVE  FLOW  VERSUS  DIFFERENTIAL  CURRENT 


F^«3CXX)psi 


FIGURE  58.  BLOCKED- LOAD  FLOW  VERSUS  FREQUENCY 


performance-characteristic  distribution  at  a  particular  point  in  time.  The  full  shape  of 
the  distribution  is  not  known  unless  some  specific  form  such  as  the  normal  is  assumed. 
The  success  of  the  method  depends  on  accurate  knowledge  of  the  form  of  the  para.neter 
distributions  involved,  and  also  on  the  validity  and  depth  of  the  mathematical  model  used 
to  represent  the  device  or  system  under  analysis.  In  order  to  extend  the  Moment  method 
to  take  into  account  the  passage  of  time,  either  experimental  or  theoretical  estimates  of 
time  rates  of  parameter  degradation  must  be  secured.  One  of  the  objectives  of  the  pres¬ 
ent  life  tests  was  to  determine  empirically  the  long-term  rates  of  change  of  servo-valve 
parameters,  such  as  quiescent  leakage  and  frequency  response.  In  addition,  the  com- 
parative*performances  of  the  modified  and  unmodified  valves  were  expected  to  allow 
inferences  about  the  extent  to  which  predicted  and  realized  reliability  improvements  were 
matched  in  practice.  In  formulating  the  design  changes  incorporated  in  the  two  modified 
valves,  ways  to  improve  reliability  were  found  by  noting  where  in  the  design  an  applica¬ 
tion  of  one  or  more  of  three  cardinal  criteria  could  be  applied.  The  first  criterion  was 
the  reduction  of  the  values  of  the  partial  derivatives  representing  the  rate  of  change  of  a 
given  performance  characteristic  with  respect  to  a  physical  parameter  of  the  valve.  For 
instance,  if  the  maximum  quiescent  leakage  can  be  made  less  sensitive  to  spool  metering 
edge  location  (without  affecting  other  performance  characteristics),  then  the  design 
change  producing  this  reduction  of  sensitivity  may  be  said  to  increase  the  relative  relia- 
bil’  ’’C  valve  with  respect  to  statistical  variations  in  the  longitudinal  spacing  of  the 

^  Iges. 

.>econd  criterion  for  reliability  improvement  was  the  tightening  of  the  nominal 
specifications  of  parameters,  so  that  the  probability  of  their  being  a  given  amount  off 
ideal  to  begin  with  would  be  lessened.  The  variations  in  physical  parameters  may  be 
thought  of  as  operating  on  performance  characteristics  by  means  of  a  variable  leverage 
factor  —  the  leverage  factor  being  the  partial  derivative  of  the  characteristic  with  respect 
to  the  parameter.  A  reduction  in  either  the  original  parameter  variation  or  the  partial 
derivative  relating  it  to  a  performance  characteristic  can  effect  an  improvement  in 
reliability. 

The  third  primary  means  for  reliability  improvement  was  the  reduction  of  the  time 
rate  of  change  of  physical  parameters.  This  area  of  improvement  is  usually  concerned 
with  the  wear  and  aging  of  parts.  It  may  be  considered  an  extension  of  the  notion  of 
initial  parameter  variation.  Not  only  is  the  original  spread  or  distribution  of  a  tolerance 
dimension  or  property  of  concern,  but  also  the  further  broadening  of  the  range  of  possible 
values  due  to  hours  of  operation  or  hours  of  shelf  time. 


Test  Results.  As  the  bench  tests  progressed,  it  was  possible  to  estimate  the  rela¬ 
tive  repeatability  of  the  tests  by  noting  the  scatter  evident  in  the  trend  curves,  and  the 
geometric  similarities  between  curves  taken  in  adjacent  series.  The  repeatability  of  all 
tests  appeared  to  be  good,  with  the  exception  of  the  differential-load-pressure  test. 

During  this  test,  apparent  instabilities  in  valve  operation  often  occurred,  and  control  of 
the  d-c  servo  amplifier  to  the  required  degree  of  precision  was  difficult.  An  example  of 
instability  in  the  differential-load-pressure  test  is  shown  in  Figure  59.  The  vertical 
section  of  the  curve  represents  a  region  in  which  continuous  control  of  the  valve  was  not 
possible  because  of  a  sudden  shift  of  the  spool.  The  shift  occurred  repeatedly  and  seemed 
independent  of  the  late  at  which  the  differential  current  was  increased. 

With  the  exception  of  pres  sure- gain  curves,  all  the  bench  tests  performed  seemed 
sufficiently  repeatable  to  justify  confidence  in  using  successive  sets  of  curves  to  trace 
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Lood  Differentiol,  psi 


FIGURE  59.  DIFFERENTIAL  LOAD  PRESSURE  VERSUS  DIFFERENTIAL  CURRENT 


parameter  changes  with  time  by  means  of  trend  curves.  The  trend  curves  were  made 
by  plotting  some  significant  feature  of  each  graph  for  a  given  valve  against  total  cycles. 
For  example,  the  tre.id  curve  for  quiescent  flow  shows  the  peak  value  of  this  performance 
characteristic  as  found  at  each  testing  period.  The  testing  periods  occurred  at  0; 

150, 000;  300, 000;  900, 000;  1, 500, 000;  3,  000, 000;  and  8,  600, OCO  cycles,  A  set  of  trend 
curve's  for  a  given  test  requirv^s  four  graphs,  one  for  each  valve.  Six  sets  of  trend 
curves  were  plotted.  They  illustrate  peak  quiescent  flow,  quiescent  spool  leakage, 
quiescent  nozzle  flow,  frequency  response  at  100  cps,  pressure  gain  ir.  the  Interval  from 
1000  to  2000-psi  blocked- load  differential,  and  null  bias.  The  set  of  quiescent-flow- 
trend  curves  was  expected  to  be  most  significant  in  efforts  to  relate  actual  to  predicted 
reliability. 


Quiescent  Flow,  The  trend  analysis  of  quiescent  flow  empio/ed  three  types  of 
curves.  The  first  was  a  plot  of  peak  quiescent  flow  for  each  test  series.  The  second 
was  a  graph  of  the  difference  between  the  peak  flow  at  null  and  the  nearly  constant  flow 
beyond  4-  or  5-ma  differential  current  on  either  side  of  null.  Following  the  reasoning 
given  in  the  description  of  the  quiescent  flow  test,  this  difference  may  be  considered  to 
represent  the  portion  of  the  maximum  quiescent  flow  accounted  for  by  leakage  past  the 
underlapped-spool  metering  edges.  It  should  be  the  most  sensitive  indicator  of  metering 
edge  wear  available,  aside  from  direct  measurement.  The  third  quiescent-flow  trend 
curve  showed  the  flow  at  8-me,  differential  current  for  each  test  series.  Since  the  8-ma 
point  is  quite  far  from  null,  it  is  on  the  constant  portion  of  the  quiescent-flow  test  curve, 
and  may  be  taken  to  represent  nozzle  flow.  The  slope  of  nozzle-flow  trend  curve  is  an 
indicator  of  nozzle  wear,  and  was  used  to  check  this  parameter. 

Figure  60  illustrates  a  typical  peak  quiescent- flow  trend  curve.  It  shows  a  period 
of  rapid  initial  increase  extending  through  the  first  2  or  3  million  cycles,  and  a  gradual 
leveling  off  into  a  much  lower  long-term  rate  of  increase.  The  rapidly  increasing  sec¬ 
tion  probably  represents  an  initial  rounding  of  the  spool  metering  edge.  Because  the 
rounded  edge  presents  more  and  more  area  to  the  flow  after  each  increment  of  erosion, 
it  is  to  be  expected  that  the  rate  of  wear  will  decrease  and  produce  a  peak  leakage  that 
changes  very  little  as  the  life  testing  progresses.  Comparison  of  the  curves  for  modified 
and  unmodified  valves  showed  that  no  definite  conclusions  may  be  drawn  concerning 
relative  reliabilities,  at  least  on  the  basis  of  the  small  sample  treated  in  this  project. 
There  were  some  marked  differences  in  the  degradation  of  the  valves,  but  some  of  these 
differences  were  split  among  the  valves  so  that  they  must  be  considered  intrasample  and 
not  intersample  phenomena.  I  ur  instance,  both  unmodified  valves  showed  an  increase 
in  total  quiescent  flow  of  about  0.  1  gpm  over  the  entire  test  program,  while  one  modified 
valve  increased  by  only  0.  05  gpm,  and  the  other  by  nearly  0.  2  gpm.  The  final  slopes  of 
the  quiescent-flow  trend  curves  for  three  valves  were  similar,  and  a  fourth,  a  modified 
valve,  was  somewhat  steeper. 

The  second  set  of  trend  curves  associated  with  quiescent  flow  was  that  of  spool 
leakage.  This  quantity  is  the  differenc-  between  nozzle  flow  and  peak  quiescent  flow. 
Since  it  turned  out  that  the  nozzles  degraded  only  a  negligible  amount,  the  spool-leakage 
curve  was  nearly  equal  to  the  total  quiescent  curve  minus  a  constant  value.  Hence,  the 
conclusions  drawn  from  the  total  quiescent- flow  curves  also  apply  to  the  spool- leakage 
curves.  Typical  spool-leakage  and  nozzle-flow  trend  curves  are  shown  in  Figure  61  and 
62.  Spool  leakage  is  plotted  against  total  hours  rather  than  total  cycles,  because  it  was 
felt  that  spool  metering  edge  wear  (the  primary  cause  of  spool  leakage  at  null)  would  be 
depende:it  on  total  hours  of  flow  rather  accumulated  cycles.  Hours  and  cycles  were  not 
uniformly  related  throughout  the  life  test,  because  the  frequency  of  operation  was 
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increased  during  the  latter  portion.  None  of  the  valves  exceeded  the  quiescent-flow 
specification  limit  during  the  *eat  program. 

Conclusions  concerning  the  quiescent-flow  trend  curves  are  limited  by  the  rather 
low  final  steady  rate  of  spool  -leakage  degradation,  and  the  extremely  small  rate  of  nozzle 
degradation.  It  was  hoped  that  renioval  of  the  2-micron  filter  from  the  life-testing  cir¬ 
cuit  would  increase  contamination  considerably  beyond  that  experienced  in  previous  tests. 
An  increase  by  a  factor  of  two  did  occur,  but  the  degradation  rate  was  still  slower  than 
originally  hoped  for.  The  possibility  of  adding  contamination  or  switching  to  a  contami¬ 
nation  stand  was  considered,  but  it  was  decided  that  this  type  of  forced  general  wear  is 
not  realistic  when  compared  with  the  actual  conditions  of  use  for  the  valves,  and  in  fact 
might  not  proceed  with  the  tame  relative  balance  as  the  slow  wear  associated  with 
systems  as  used  in  nortrkal  applications.  In  addition,  time  and  cost  considerations  made 
it  Impractical  to  use  a  contamination- stand  prucadure  in  which  only  one  or  two  valves 
could  be  run  at  a  time.  The  cumulatlve-contamination-count  curve  for  the  present  tests 
is  shown  in  Figure  63,  The  contamination  count  was  made  on  a  Coultier  counter. 

Although  contaminatior  counts  obtained  by  this  method  are  not  as  accurate  as  obtained 
by  more  laborious  methods,  the  results  appear  to  be  consistent  with  the  filter  system 
that  was  used. 


Frequency  Response.  The  frequency-response  trend  curves  showed  that,  generally 
speaking,  the  response  of  the  modified  valvee  at  100  cps  was  better  than  the  response  of 
the  unmodified  valvee.  At  the  conclusion  of  the  life-testing  period,  the  two  unmodified 
valves  were  down  7  and  3  db,  while  the  modified  were  down  1  and  2  db.  Initially,  both 
unmodified  valvet  were  down  approximately  4  db,  and  both  modified  were  down  about 
2  db.  During  the  tests,  then,  the  one  unmodified  lost  3  db  of  response  and  one  gained 
1  db;  for  the  modified  valves,  one  remalred  the  same  and  one  also  gained  1  db.  With 
samples  as  small  as  two  units,  it  cannot  necessarily  be  concluded  that  the  modified  valves 
are  inherently  better  in  frequency  response,  although  such  a  conclusion  was  not  contra¬ 
dicted  by  the  actual  valves  tested.  A  typical  frequency-response  trend  curve  is  shown  in 
Figure  64.  All  valves  demonstrated  an  Initial  dip  in  response,  but  for  two  of  the  valves 
(one  modified  and  one  unmodified}  the  dip  was  considerably  more  pronounced  than  that 
shuwn. 


Pressure  Gain.  The  pressure-gf.ln  trend  curves  were,  unfortunately,  Indicative 
of  the  instabilities  and  other  difficulties  sneountered  during  the  differential-load-pressure 
versus  different.al-current  tests.  Only  one  valve  (an  unmodified)  showed  the  initial 
increase  and  subsequent  decrease  that  would  be  expected  in  the  light  of  experience  gained 
in  a  previous  test  program.  Two  of  the  vslvspi,  one  modifie  .d  one  unmodified,  showed 
a  long-term  Increase  in  pressure  gain,  which  is  difficult  to  sxplsin  in  view  of  the  slow 
but  normal  degr-^dation  evidenesd  by  ths  othsr  tssts.  It  was  fslt  that  ths  pressure-gain 
trend  curves  must  be  regarded  ss  inconclusive  data  for  this  test  ssriss. 


FIGURE  64.  FREQUENCY  RESPONSE  (BLOCKED  LOAD)  AT  100  CPS  VERSUS  TOTAL  CYCLES 


Visual  Inspection  of  Valves 


In  order  to  bring  to  light  any  specific  features  that  might  account  for  some  of  the 
tost  results,  all  of  the  valves  were  disassembled  and  inspected  after  the  conclusion  of 
the  tests.  Table  16  shows  a  summary  of  the  visual  inspection.  The  two  most  significant 
findings  were  the  definite  flats  worn  on  the  ball  ends  of  the  feedback  wire  of  the  unmodi¬ 
fied  valves.  The  wear  on  the  modified-valve  feedback  wires  was,  by  contrast,  not 
detectable  with  a  micrometer.  Figure  65  is  an  enlarged  photograph  of  the  ball  tips  of  a 
modified  and  an  unmodified  valve  at  the  conclusion  of  the  tests.  The  magnitude  of  the 
feedback-wire  wear  on  the  unmodified  valves  was  0.001  in.,  or  approximately  3  percent 
of  the  tip  ball  diameter.  The  cause  of  this  difference  in  wear  is  difficult  to  ascertain 
with  certainty.  The  feedback-wire  spring  rate  of  the  modified  design  is  lower,  but  this 
is  balanced  by  a  longer  stroke  required  by  the  narrower  ports.  In  fact,  the  spring  rate 
was  reduced  in  order  to  produce  the  longer  stroke.  It  is  felt  that  the  differences  in  wear 
rate  are  m.ost  probably  the  result  of  different  hardnesses  on  the  ball  surfaces.  According 
to  the  valve  manufacturer,  this  is  the  most  likely  explanation.  If  the  feedback-wire  wear 
reduction  is  an  inherent  feature  of  the  modified  design,  it  may  be  significant  in  improving 
the  performance  of  the  valve  near  null  in  the  later  stages  of  valve  life. 


FIGURE  65.  BALL  ENDS  OF  FEEDBACK  V/IRE  ON  MODIFIED  AND 
UNMODIFIED  VALVE  DESIGNS  AFTER  8,  600,  000 
CYCLES  OF  OPERATION 
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TABLE  XVI. 


SERVO-VALVE  VISUAL  INSPECTION 


Valve  No, 

Spool 

Bushing 

Feedback 

Wire 

Other 

65295 

(Unmodified) 

Very  slight 
scoring,  one 
outside  land 

ok 

Worn  O.UOl  in.  on 
ball  (Ball  diame¬ 
ter  0.  03 1  in. ) 

Damaged  O-ring 
on  one  bushing 
plug 

65296 

(Unmodified) 

Moderate  scoring, 
one  outside 
land 

ok 

Worn  0.001  in.  on 
ball 

Damaged  O-ring 
on  one  bushing 
plug 

65297 

(Modified) 

Very  slight 
scoring,  cen¬ 
ter  and  out¬ 
side  lands 

ok 

No  wear  on  ball 

65298 

(Modified) 

Heavy  scoring  on 
middle  and  at 
one  end.  One 
short,  very  deep 
score  at  end  of 
outside  meter¬ 
ing  edge 

ok 

No  wear  on  ball 

•mm* 

One  possible  effect  of  the  feedback- wire  wear  experienced  in  the  unmodified  valves 
is  a  tendency  for  the  valve  to  be  unstable  near  null  with  a  blocked  load.  The  test  best 
able  to  reveal  such  a  condition  is  the  differential-pressure  versus  differential-current 
test.  A  check  of  the  original  data  sheets  showed  that  only  one  instance  of  "jumping" 
instability  occurred  with  the  modified  valves,  as  opposed  to  three  for  the  unmodified' 
valves.  The  discovery  of  a  possibly  important  but  unanticipated  difference  between  the 
designs,  such  as  the  difference  in  feedback-wire  wear,  was  one  of  the  goals  of  the  testing 
program.  Such  findings  help  to  compensate  for  the  inconclusiveness  arising  from  the 
low  statistical  value  of  the  small  samples  involved. 

The  other  anomalies  noted  in  the  visual  inspection  table  —  a  rather  heavy  spool 
scoring  in  one  valve,  and  damaged  O-rings  on  two  valves  —did  not  appear  to  cause 
noticeable  defects  in  operation. 


Conclusions 


The  first  part  of  die  analytical  study  employed  an  analog  computer  simulation  of  a 
two-stage  servo  valve  to  obtain  partial  derivatives  of  various  performance  characteristics 
with  respect  to  valve  parameters  that  are  susceptible  to  degradation.  The  partial  deriva¬ 
tives  were  calculated  and  plotted  for  different  values  of  key  design  parameters.  These 
plots  were  evaluated  and  their  interrelationships  noted  in  order  to  identify  design  changes 


that  would  reduce  the  over-all  sensitivity  of  the  valve  to  parameter  degradation  and  to 
the  effects  of  parameter  variations  that  are  introduced  in  the  manufacturing  process. 

The  study  showed  that  the  existing  design  was  already  close  to  optimum  as  far  as 
the  critical  areas  of  degradation  were  concerned,  and  it  was  not  expected  that  design 
changes  would  bring  about  significant  improvement  in  the  gross  reliability  of  the  valve. 
The  critical  areas  were  identified  in  previous  programs  as  those  parts  of  the  spool  and 
bushing  that  influence  the  quiescent  leakage  at  null. 

I 

The  design  changes  that  were  ultimately  derived  did  reduce  slightly  the  sensitivity 
of  quiescent  leakage  to  degradation  effects,  and  also  reduced  the  initial  value  of  quiescent 
leakage  and  improved  the  frequency  response  of  the  valve.  In  addition,  these  modifica¬ 
tions  reduced  significantly  the  sensitivity  of  several  performance  characteristics  to 
degradations  in  the  first-stage  orifice  and  nozzle  flow  characteristics.  Previous  experi¬ 
ments,  however,  have  shown  that  spool  and  orifice  wear  are  not  significant  under  normal 
operating  conditions. 

Even  though  the  improvement  in  anticipated  valve  reliability  resulting  from  the 
derived  design  modifications  was  relatively  minor,  the  results  of  the  study  are  signifi¬ 
cant  because  they  serve  to  demonstrate  the  successful  application  of  the  Moment  method 
of  reliability  analysis  to  the  redesign  and  improvement  of  an  existing  control- system 
component. 

The  prediction  of  the  expected  reliability  improvement  in  the  valve  was  based  on  the 
performance  characteristic  of  quiescent  leakage  at  null. 

In  the  reliability  improvement  derivation  for  quiescent  leakage,  it  was  shown  that 
in  the  statistical  sense,  i.  e. ,  for  the  averaged  histories  of  a  large  number  of  valves,  a 
definite  improvement  in  time  to  failure  with  regard  to  allowable  leakage  would  occur. 
However,  the  study  also  showed  that,  given  the  normal  parameter  variations  expected 
among  valves,  any  single  modified  valve  might  or  might  not  show  improvement  over  the 
existing  design,  depending  on  its  particular  combination  of  initial  parameter  values. 

Thus,  it  was  not  possible  to  predict  that  a  small  sampling  of  old  and  new  servo  valves 
would  show  the  relative  improvements  in  average  characteristics  that  were  predicted  for 
large  samples.  Hence,  it  was  concluded  that  a  testing  program  on  a  small  sample  of 
valves  might  be  more  useful  in  showing  up  unexpected  differences  between  the  designs 
than  in  verifying  quantitatively  the  predicted  improvements  in  drift-failure  reliability. 

When  the  experimental  life-te sting  program  was  carried  out  on  a  sample  of  two 
modified  and  two  unmodified  valves,  the  above  conclusions  were  generally  confirmed. 

The  spread  in  initial  characteristics  of  the  valves  was  relatively  quite  great,  both  in  Che 
inter  sample  and  in  the  intrasample  sense.  Although  it  was  not  possible  to  draw  firm 
conclusions  regarding  relative  reliability,  an  Interesting  unexpected  difference  did  show 
up.  This  difference  was  a  noticeable  improvement  in  the  wear  of  the  ball  on  the  end  of  the 
feedback  wire  in  the  new  design.  This  improvement  would  lead  to  significantly  better 
performance  of  the  "aged"  valve  in  the  critical  area  very  close  to  null.  It  is  believed 
that  this  wear  reduction  is  the  result  of  increased  hardness  of  the  ball  tip  in  the  new 
design. 
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•  S-ABSTRACT 

JI^A  previous  reliability  research  effortiiSy  Battell^generated  an  analytical 
technique  called  the  Moment  method  for  predicting  drift-type  failures  in  com¬ 
ponents  and  systems.  This  technique  is  based  on  the  Propagation  of  Variance 
formula,  which  combines  data  derived  from  the  system,  or  system  modcd,  and 
parameter  tolerance  data,  to  formulate  an  estimate  of  the  probability  that  the 
system  will  perform  within  specifications  at  the  initial  or  ^me  future  time  of 
operation.- . . p. 

This  report  discusse^the  application  of  the  reliability-prediction  technique 
to  a  momentum- exchange- type  attitude  control  and  stabilisation  subsystem  for  a 
satellite,  a  liquid-metal  servo  actuating  subsystem,  and  an  electrohydraulic  servo 
valve. 

mathematical  model  was  completed  for  the  mom sntum-exc)j>nge- type  atti¬ 
tude  control  (fluid  flywheel),  significant  system  parameters  identified,  the  attalysis 
of  the  expected  system  degradation  over  time  was  made,  and  a  prediction  of  ef¬ 
ficiency  and  reliability  is  presentsd(J|i  the  report. 
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